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AN  EVALUATION  OF  NONSINUSOIDAL  RADAR  TECHNIQUES 


1.0 

INTRODUCTION 

During  the  past  decade,  a  seemingly  insignificant,  but  highly 
emotional,  controversy  has  existed  on  the  use  of  nonsinusoidal 
"carrier-free"  waveforms  for  applications  in  radar  and  communication. 
The  debate  apparently  started  on  the  utility  of  Walsh  functions  as 
compared  to  sinusoids  [1,2];  it  later  evolved  to  the  system  Implications 
of  the  utility  of  the  large  relative-bandwidth  carrier-free  waveforms  as 
contrasted  to  the  narrow  relative-bandwidth  of  conventional 
carrier-based  systems  [33,53,54].  Indeed,  the  debate  continued  to  such 
a  point  that  proponents  of  the  nonsinusoidal  approach  only  published 
their  work  in  the  IEEE  Transactions  on  Electromagnetic  Compatibility; 
this  action  further  fueled  the  skepticism  of  opponents  since  the 
objectivity  of  such  a  small  peer  group  could  be  easily  questioned.  As  a 
result,  there  had  been  no  objective  comparison  made  between  carrier-free 
and  conventional  radar  techniques  for  military  applications.  Such  an 
assessment  was  needed  by  the  Department  of  Defense  (DoD)  in  order  to 
understand  the  potential  of  nonsi r usoidal  radar  to  provide  unique 
capabilities  not  achievable  with  conventional  techniques.  For  example, 
it  had  been  argued  that  nonsinusoidal  radar  could  be  used  as  an 
anti-stealth  system;  i.e.,  that  it  could  be  used  to  detect  and  track 
low-observable  platforms  [46].  A  validation  of  that  claim  alone  was  of 
sufficient  importance  to  merit  an  objective  evaluation;  however,  there 
were  undoubtedly  other  important  applications  for  which  nonsinusoidal 
radar  may  be  well  suited.  The  heated  controversy  which  existed  over  a 
decade,  combined  with  the  need  by  the  DoD  to  understand  the  military 
potential  of  such  a  radar,  provided  the  motivation  for  the  study  whose 
results  are  presented  herein. 


1.1  REVIEW  OF  THE  NONSINUSOIDAL  RADAR  CONTROVERSY 


This  section  presents  a  brief  overview  of  the  "common"  arguments 
for  and  against  the  use  of  nonsinusoidal  waves  in  radar  applications. 
It  provides  an  important  historical  perspective  of  the  controversy  and 
thus  also  provides  insight  into  why  the  issues  were  not  easily  resolved, 

1.1.1  THE  CASE  FOR  NONSINUSOIDAL  WAVES 

The  beginning  arguments  for  nonsinusoidal-based  systems  apparently 
goes  back  to  the  discovery  of  a  set  of  orthogonal  functions  —  Walsh 
functions  --  that  could  be  used  as  a  basis  set  for  waveform  analysis, 
coding,  etc.,  in  lieu  of  the  traditional  sinusoids.  Since  the  functions 
were  binary  valued,  it  was  felt  that  the  computation  of  the  required 
coefficients  could  be  done  much  more  efficiently  than  in  the  case  of 
sinusoids.  Even  the  pros  and  cons  of  Walsh  function  representation  were 
presented  in  the  midst  of  heated  debate  [1,2],  and  the  IEEE  Transactions 
of  Electromagnetic  Compatibility  became  the  forum  for  such  papers. 

Soon  after  that,  H.F.  Harmuth  developed  the  notion  that  all  of 
electrical  engineering  was  so  biased  by  the  use  of  sinusoids,  starting 
with  the  practical  hardware  limitations  faced  by  Mar:oni ,  that  alternate 
system  concepts,  based  upon  nonsinusoidal  waves,  never  received  adequate 
treatment.  This  theme  is  repeated  throughout  his  paoers  [32-51]  and  is 
adequately  summarized  in  his  recent  book  [52]. 

Current  views  in  favor  of  nonsinusoidal-based  sy items  are  no  longer 
exclusively  tied  to  the  merits  of  Walsh  functions  alone,  but  to  the 
notion  of  the  large  relative  bandwidths  which  can  be  achieved  with 
"carrier-free"  waveforms.  Here,  relative  bandwidth  n  is  defined  by 


where  f^(f^)  1s  highest  (lowest)  significant  frequency  component  in 
the  signal  of  interest.*  Typical  audio  signals  and  pulse  video 
waveforms  have  n  =  0.75  to  1,  while  <.  0.1  is  typical  for  radar.  Thus, 
it  is  argued  that  to  obtain  the  large  bandwidths  needed  for  high 
resolution,  conventional  radars  must  operate  in  the  millimeter  wave 
region  (e.g.,  35  or  95  GHz)  with  the  attendant  problems  of  high 
absorptive  losses,  high  noise  temperatures,  and  strong  dependence  on 
weather  effects  [33].  For  example,  a  1  nanosecond  pulse  system  would 
have  significant  spectral  content  up  to  1  GHz;  a  conventional  system 
with  the  same  bandwidth  would  require  either  a  35  or  95  GHz  carrier. 

With  these  concepts  in  mind,  the  proponents  for  nonsinusoidal  radar 
present  an  impressive  list  of  potential  military  applications;  it,  is 
provided  in  Table  1-1.  Note  that  all  the  justifications  are  based,  in 
part,  on  the  wide  relative  bandwidth  of  the  technique  and  on  the  fact 
that  many  limiting  propagation  effects  can  be  ameliorated  by  operating 
in  a  lower  frequency  band.  A  few  quantitative  analyses  have  been 
presented  [32-51  ],  but  the  results  are  technically  incomplete  as 
compared  to  similar  system  studies  using  conventional  techniques. 

It  is  interesting  to  note  here  that  numerous  nonsinusoidal  radars 
have  been  built  commercially  for  geophysical  probing  applications.  In 
fact,  one  such  radar  was  used  to  detect  the  tail  section  of  the  Air 
Florida  Flight  90  which  crashed  Into  the  Potomac  several  years  ago. 
Positive  results  as  these  tend  to  make  proponents  all  the  more  steadfast 
in  their  views.  In  fact,  the  public  claims  for  utility  are  so  often 

*In  the  context  of  the  present  study,  nonsinusoidal  radar  (NSR)  and 
wideband  radar  (WBR)  were  used  interchangeably  for  radar  systems  whose 
waveforms  were  of  a  large  relative  bandwidth.  It  was  judged  that  the 
fundamental  issue  was  that  the  waveform  had  a  large  relative  bandwidth 
and  not  that  it  had  to  be  considered  carrier-free.  In  fact,  Appendix  A 
discusses  that  the  concept  of  a  carrier  frequency  is  completely  arbi¬ 
trary  when  using  complex  notation  to  describe  a  waveform  with  large 
relative  bandwidths. 
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TABLE  1-1.  POTENTIAL  MILITARY  APPLICATIONS  OF  NONSINUSOIDAL  RADAR 


MILITARY  APPLICATION 
Anti-stealth  Radar 


All-weather,  High-Resolution 
Llne-of-SIght  (LOS)  Radar 

Antl-balllstlc  Missile  (ABM) 
Defense  Radar 

Over-the-Horlzon  (OTH)  Radar 

Look- Down- Radar  for  detecting 
low-cross  section  platforms 

Synthetic  Aperture  Radar 


REASON 

Wide  bandwidth  easily  penetrate 
coating  materials. 

Hard  to  defeat  with  active 
stealth  techniques. 

Excites  (unknown)  body 
resonances. 

Improved  resolution  with  wide 
bandwidth. 

Better  performance  In  rain,  fog, 
snow,  since  spectrum  below 
5  GHz. 

Improved  resolution  to  separate 
decoys  from  warheads. 

Microsecond  pulses  provide 
better  resolution  with 
reduced  dispersion. 

Better  range  resolution  to 
reduce  clutter  effects. 

High  bandwidth  enhances  radar 
cross-section. 

Improved  range  resolution  allows 
for  enhanced  Image  resolution. 


overstated  that  their  credibility  Is  In  question.  (It  Is  important  to 
observe  that  radar  remote  sensing  is  a  "soft  science"  and,  as  such,  is 
not  conclusive  support  for  complete  success  in  military  applications.) 
In  short,  nonsinusoidal  techniques  for  radar  have  been  advocated  by  a 
number  of  individuals.  The  advantages  have  probably  been  overstated 
based  upon  heuristic  reasoning  and  have  not  been  quantitatively  proven. 

1.1.2  THE  CASE  AGAINST  NONSINUSOIDAL  WAVES 

The  proponents  of  a  conventional,  sinusoidal-based  radar  have 
argued  against  the  nonsinusoidal  radar  based  primarily  upon  hardware  and 
complexity  issues.  Stated  simply,  they  say  that  since  the 
implementation  is  so  difficult  and  since  there  is  no  solid  prediction  of 
an  attendant  gain  in  performance,  then  there  is  no  motivation  to  proceed 
further.  For  the  record,  Table  1-2  presents  a  partial  list  of  the 
hardware  and  system  related  difficulties.  In  all  fairness,  the 
conclusions  are  also  heurlstical ly  correct  but  quantitatively  soft. 
Without  specific  system  applications  and  performance  parameters,  it  is 
difficult  to  predict  which  areas  of  difficulty  are  "fundamental"  or 
which  are  just  perceived  based  upon  current  design  approaches.  It  is 
exactly  that  mismatch  of  glowing  applications  and  sour  implementation 
issues  that  motivated  this  study. 

1.2  PROGRAM  OBJECTIVE  AND  APPROACH 

With  the  nonsinusoidal  controversy  as  a  background,  the  objective 
of  the  subject  study  was  to  provide  an  unbiased  evaluation  of 
nonsinusoidal  radar  (NSR)  techniques  for  military  applications  and  to 
identify  future  activities  needed  to  realize  the  potential  gain  of  NSR. 
The  study  objective  was  accomplished  by  a  combination  of  (1)  interviews 
with  several  of  the  key  NSR  proponents,  including  H.F.  Harmuth,  (2) 
original  radar  system  modeling  and  performance  prediction  calculations. 


TABLE  1-2.  DIFFICULTIES  OF  NONSINUSQIDAL- BASED  RADAR 


RESTRICTIVE  ELEMENT 


COMMENTS 


Antenna 


Waveform  Generator/ 
Power  Ampl Ifier 

Propagation  Media 


Timing/Synchronization 


Receiver  Processing 


Component  (Front-end) 
Dynamic  Range 


Range  Resolution 


Poor  radiative  efficiency. 

Poor  directivity. 

Phase  center  varies  In 
time/space. 

Dispersive  nature  shapes 
pulse  spectrum. 

Requires  extreme  power. 

Difficult  switching  requirments. 

Dispersive/absorptive  media 
destroys  specific  pulse 
structure. 

Extremely  tight  tolerances 
needed  for  resolution  and 
reduction  of  "co-channel" 
interference. 

Analog  tapped  delay  line  with 
very  tight  tolerances  and 
thousands  of  taps. 

Must  be  maintained  without 
destroying  temporal 
character. 

Too  much.  May  unnecessarily 
resolve  individual  scattering 
centers  of  a  target. 


and  (3)  critical  review  and  critiques  of  numerous  published  literature 
on  NSR  and  other  related  topics.  Great  care  was  taken  to  ensure  that 
fundamental  Issues  were  not  confused  with  current  hardware  limitations 
so  that  "pure"  conclusions  could  be  deduced  from  the  analysis. 

1.3  FINAL  REPORT  ORGANIZATION  AND  RESULTS 

This  report  is  organized  as  follows:  Section  2.0  provides  a 
summary  of  radar  system  requirements  and  conventional  radar  limitations 
for  military  applications.  Emphasis  is  placed  on  target  detection,  and 
tracking  In  a  heavy  clutter  environment.  The  study  scope  was  limited  to 
this  application  because,  for  certain  military  scenarios,  it  is  still 
very  stressing  for  conventional  radars  and  there  appeared  to  ne  few 
additional  alternatives  remaining.  Moreover,  it  was  judged  that  NSR 
techniques  could  be  promising  in  that  case. 

Section  3.0  presents  an  end-to-end  analysis  and  radar  system  model 
for  detection  of  stationary  targets  in  clutter.  Using  crude,  but  useful 
models  for  clutter  and  target  frequency  dependence,  optimum  radar 
waveforms  are  designed  to  maximize  the  target  detectability.  In  some 
instances,  the  formalism  dictates  that  optimum  waveforms  should  have 
large  relative  bandwldths.  The  Importance  of  the  result  is  not  the 
specific  signal  design,  but  rather  the  "existance  proof"  that  signals 
with  wide  relative  bandwidth  may  be  optimum  for  detection  of  targets  in 
clutter.  This  is  particularly  true  when  the  target  response  and/or  the 
intervening  propagation  media  have  significant  frequency  dependencies. 

Section  4.0  presents  critiques  of  many  relevant  papers  published  in 
the  area  of  WBR  techniques.  Among  many  results,  it  demonstrates  that 
the  use  of  singularity-expansion  methods  (SEM)  to  describe  target 
scattering  provide  a  means  to  do  target  classification.  (SEM  techniques 
require  WBR  waveforms  and  are  virtually  unknown  to  the  NSR  proponents.) 
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Finally,  Section  5.0  presents  a  summary  of  conlusions  and  recommen¬ 
dations  for  future  activities.  It  Is  concluded  that  the  formalism 
developed  In  Section  3.0  should  be  applied  In  detail  to  specific 
military  problems  to  Identify  one  or  two  most  promising  for  which 
additional  system/component  hardware  development  work  should  be 
directed. 
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2.0 

RADAR  SYSTEM  CONSIDERATIONS 


2.1  GENERAL  COMMENTS 

The  potential  uses  of  radar  for  military  applications  are  quite 
numerous  so  that  to  provide  an  extensive  description  Is  beyond  the  scope 
of  the  present  study.  However,  It  is  useful  to  present  some  general 
considerations  to  focus  the  deliberation  for  NSR  techniques.  In  the 
weapon  system  context,  the  required  functions  are  (1)  surveillance  of  an 
area  of  interest  for  potential  military  targets,  (2)  tracking  of 
potential  targets  of  interest,  (3)  classification/  identification  of 
tracked  target,  (4)  fire  control  for  the  weapon  delivery,  and  (5)  damage 
assessment  after  weapon  delivery.  Of  course,  these  categories  are 
"soft"  depending  upon  the  specific  mission,  e.g.,  the  tactical  problem 
of  destroying  armor  on  a  road,  or  the  strategic  problem  of  destroying 
missile  launch  sites,  or  defending  specific  area  against  an  air  attack. 

The  type  of  radar  required  to  do  a  specific  function  depends 
strongly  on  the  location/p'latform  constraints  of  the  radar  and  expected 
targets/backgrounds.  Radars  are  routinely  considered  for  basing  on 
satellites,  aircraft  (both  large  AWACS  and  fighter-type  aircraft),  ships 
and  ground  installations.  The  platforms  typically  dictate  limitations 
on  prime  power  and  available  antenna  aperture  sizes  and  the  areas  to  be 
searched  (relative  to  the  platform  motion/location).  In  addition,  the 
radar  basing  will  dictate  the  type  of  clutter  background  against  which 
the  target  will  be  observed. 

Depending  on  the  mission,  the  target  type  will  also  dictate  the 
available  conventional  (or  narrowband)  radar  techniques  which  are 
applicable  for  a  specific  function.  For  example,  an  airborne  radar 
searching  for  a  moving  tank  column  on  a  road  can  use  MTI  techniques, 
whereas  detection  of  stationary  tanks  would  require  either  CFAR  or 
imaging  (spatial)  discrimination  techniques. 
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It  was  decided  that  the  ultimate  end-to-end  system  analysis 
(presented  In  Section  3.0)  would  need  tc  be  limited  to  a  specific 
scenario  to  yield  quantitatively  useful  results.  Thus,  the  detection  of 
a  target  In  clutter  was  selected  to  be  the  most  interesting  problem.  It 
Is  representative  for  a  number  of  scenarios,  Including  an  autonomous 
missile  directed  against  tanks,  an  AWACS  searching  for  penetrating  enemy 
aircraft,  as  well  as  a  surface-to-air  missile  (SAM)  battery  defending  a 
specific  area. 

In  the  next  section,  we  review  the  use  of  conventional  (narrowband) 
radar  techniques  for  such  a  problem.  This  review  Is  necessary  to 
underscore  the  fundamental  limitations  of  conventional  system  and  to 
illustrate  the  sophisticated  signal  processing  used  in  such  systems. 

2.2  REVIEW  OF  NARROWBAND  RADAR  TECHNIQUES  FOR  CLUTTER- 

LIMITED  APPLICATIONS 

2.2.1  STATIONARY  TARGET  DETECTION 

To  motivate  the  stationary  target  detection  problem,  consider  an 
anti-armor  missile  system.  A  key  feature  of  such  a  system  is  the 
ability  to  reliably  detect,  descriminate  and  track  both  stationary  and 
moving  enemy  vehicles  such  as  tanks,  self-propelled  artillery,  APC's, 
ammunition  and  fuel  trucks,  etc.  These  systems  must  operate  in  an 
autonomous  mode  once  they  are  launched.  There  Is  no  data  link  from 
missile-to-mlssile  or  from  missile-to-laimch  aircraft,  and  the  cn-board 
missile  processor  performs  all  navigation,  search,  detection, 
classification  and  homing  functions  after  launch. 

A  typical  ai r-to-surface  concept  has  a  range  of  approximately  10  km 
and  can  be  launched  at  a  low  altitude  so  that  the  target  area  does  not 
need  to  be  visible  from  the  aircraft  at  the  time  of  launch  -  thus 
increasing  aircraft  survivability.  The  launch  aircraft  is  assumed  to  be 
vectored  to  the  launch  point  by  a  surveillance  system,  such  as  J-STARS. 


The  mapped  swath  width  and  search  length  are  chosen  to  permit  reliable 
operation  when  both  aircraft  vectoring  errors  and  missile  guidance 
errors  are  present. 

Since  most  of  the  primary  and  secondary  roads  In  Europe  and  Asia 
are  tree-lined,  it  is  seen  that  a  primary  requirement  of  such  a  missile 
system  is  to  detect  targets  against  a  harsh,  non-homogeneous  clutter 
background.  This  must  be  accomplished  when  armor  is  both  moving  and 
stationary,  for  example,  at  rest-stops  and  staging  areas  where  tree 
cover  will  be  used.  Also,  certain  vehicles  will  usually  be  stationary, 
for  example,  parked  anti-air  units  that  protect  advancing  armor.  The 
reliable  detection  of  stationary  and  moving  targets  embedded  in  clutter 
is  thus  the  primary  problem  that  must  be  solved. 

A  constant  false  alarm  rate  (CFAR)  is  usually  a  requirement  for  any 
radar  that  is  used  to  detect  stationary  or  slowly  moving  ground  targets 
embedded  in  a  clutter  background.  This  means  that  the  false  alarm  rate 
should  Ideally  be  independent  of  the  clutter  distribution  and  clutter 
statistics.  The  usual  way  that  this  is  accomplished  is  by  comparing  the 
power  returned  from  a  cell  (or  cells)  being  tested  for  the  presence  of  a 
target  with  the  power  returned  from  a  group  of  cells  surrounding  the 
test  cell(s).  If  the  ratio  exceeds  a  threshold,  then  a  potential  target 
is  declared.  The  cells  surrounding  the  potential  target  are  used  to 
statistically  characterize  the  local  clutter  background  and  the 
detection  threshold  is  chosen  (as  a  function  of  the  local  clutter 
background)  to  achieve  a  given  probably  of  false  alarm. 

The  CFAR  property  is  achieved  by  adaptively  changing  the  detection 
threshold  as  a  function  of  the  observed  clutter  statistics.  An  example 
of  this  type  of  CFAR  "window"  process  is  given  in  Figure  2-1  where  the 
clutter  returns  from  the  24  cells  surrounding  the  test  cell  are  used  to 
control  the  detection  threshold  for  the  center  target  cell.  As  first 
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Figure  2—1.  Sketch  of  Area  CFAR  Window  (5  Cells  by  5  Ceils) 
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demonstrated  by  Finn  [3]  and  Finn  and  Johnson  [4],  an  adaptive  threshold 
of  this  type  Is  a  key  requirement,  since  It  Is  possible  to  obtain  false 
alarm  rates  that  are  orders  of  magnitude  different  if  the  form  of  the 
clutter  distribution  changes  and  the  threshold  is  kept  the  same. 

Usually  the  threshold  adjustment  is  made  on  the  basis  of  estimates 
of  the  first  two  moments  of  the  local  clutter  distribution.  In  the 
detection  literature  these  are  usually  called  one  or  two  parameter  CFAR 
tests.  A  general  two-parameter  CFAR  procedure  was  first  given  by 
Goldstein  [5],  His  processor  will  maintain  a  constant  false  alarm  rate 
in  either  log-normal  or  Weibull  clutter,  however,  a  given  threshold 
setting  will  give  a  larger  Pp^  in  log-normal  clutter  as  compared  to 
Weibull  clutter  even  though  the  clutter  parameters  are  the  same.  Thus, 
his  detection  threshold  must  be  changed  as  the  clutter  distribution 
changes.  The  Goldstein  processor  does  not  Include  a  method  for 
distinguishing  clutter  types  and  adjusting  the  detection  threshold 
accordingly.  More  robust  CFAR  processors  that  are  extensions  of  the 
Goldstein  procedure  have  been  analyzed. 

The  primary  limitation  on  the  detection  performance  of  such  a 
system  is  the  effective  signal-to-clutter  that  can  be  obtained.  Present 
and  planned  tactical  CFAR  systems  usually  operate  with  a 
target-to-cl utter  ratio  of  5  to  15  dB.  As  will  be  shown,  even  a  small 
increase  (3  to  5  dB)  in  signal-to-clutter  ratio  will  have  a  large  payoff 
in  terms  of  detection  performance. 

It  should  finally  be  noted  that  any  potential  targets  that  pass  the 
detection  threshold  of  the  above  type  of  CFAR  processor  would  be 
subjected  to  additional  discrimination  tests  in  order  to  further  reduce 
the  probability  of  false  alarm  and  produce  a  final  selected  target. 
Such  tests  could  include  ranking  by  target-to-cl utter  ratio,  target 
shape  tests,  length-to-width  ratio  tests,  analysis  of  detection  patterns 
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and  edge  detections,  use  of  polarization  discriminants,  scene 
segmentation  algorithms,  taking  another  look  at  a  finer  resolution,  and 
finally,  scintillation  vs.  frequency  tests. 

All  of  these  tests  (except  perhaps  for  the  last  one)  will  also  be 
helped  by  any  increase  in  effective  target- to-cl utter  ratio  that  can  be 
achieved.  The  majority  of  the  research  being  conducted  at  the  present 
time  on  signal  processing  for  tactical  missile  seekers  is  in  the  area  of 
target  discrimination  tests  like  those  listed  above.  A  large  number  of 
both  analytically  derived  and  heuristic  algorithms  are  being  tested 
using  both  simulated  and  live  test  data. 

As  shown  in  an  article  by  Novak  and  Vote  [6],  many  complex  targets 
(assuming  a  narrowband  waveform)  may  be  characterized  by  a  two-part 
statistical  model.  At  any  given  aspect  angle,  the  target  radar  cross 
section  (RCS)  exhibits  a  Rayleigh  scintillation  in  amplitude  while  the 
average  RCS  over  many  aspect  angles  may  be  characterized  by  a  log-normal 
density  function.  If  a  frequency  diverse  (and/or  polarization  diverse) 
pulse  train  is  used  to  average  out  the  Rayleigh  scintillation  component, 
then  only  the  log-normal  variation  of  average  RCS  with  aspect  angle 
remains.  It  is  this  component  which  dominates  Pp-Ppp  performance  since 
the  aspect  angle  at  which  the  target  will  be  viewed  is  almost  always 
unknown.  Thus  PD‘PFA  performance  predictions  must  be  averaged  over 
aspect  angle.  (This  point  will  be  discussed  further  in  Section  3.7.) 

Thus,  the  first  way  in  which  a  large  relative  bandwidth  waveform 
might  assist  the  detection  problem  is  in  the  elimination  of  the  Rayleigh 
component  of  target  fluctuation.  Another  way,  which  is  probably  of 
greater  importance,  is  that  such  a  waveform  will  result  in  an  effective 
increase  in  target  cross-section  and  an  effective  decrease  in  measured 
clutter  cross-section  (assuming  the  resolution  cell  size  remains  the 
same).  The  decrease  in  the  reflection  coefficient  of  clutter  as 
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wavelength  decreases  is  a  well  documented  experimental  fact  as  is  the 
increase  in  target  cross-section  as  the  length  of  various  structures  on 
the  target  become  comparable  with  the  radiating  wavelength.  This  effect 
will  also  be  quantified  in  Section  3.0. 


The  overall  effect  of  a  wideband,  non-sinusoidal  waveform  should 
therefore  be  an  increase  in  the  effective  signal-to-clutter  or  target- 
to-cl utter  ratio.  It  is  important  to  note  that  even  a  small  increase  (3 
to  5  dB)  in  signal-to-clutter  ratio  will  have  a  large  payoff  in  terms  of 
detection  performance  (i.e.,  present  performance  is  below  the  "knee"  of 
the  detection  curve).  This  can  be  seen  on  Figure  2-2  taken  from  the 
Goldstein  paper  which  gives  curves  of  Pp  vs.  target-to-  clutter  ratio, 
TCR,  for  three  levels  of  Pp^.  Log-normal  clutter  is  assumed  with  a 
power  (cross-section)  mean-to-median  ratio  equal  to  R=2.07.  A  Rayleigh 
target  is  postulated  in  the  cell  under  test.  A  very  large  number  of 
auxiliary  clutter  cells  are  assumed  so  that  there  is  no  CFAR  loss.  The 
figure  dramatically  illustrates  how  difficult  it  is  to  detect  a  target 
in  a  severe  clutter  environment  with  small  values  of  Pp^.  This  is  a 
result  which  is  well  known  to  engineers  who  must  design  and  develop  CFAR 
processors  to  be  used  in  a  tactical  environment.  It  is  seen,  for 


example,  that  for  PpA= 
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,  an  increase  of  target-to-cutter  ratio  from 

This  increase  in  TCR 
will  also  improve  the  performance  of  non-parametric  processors  (i.e., 
those  which  utilize  order  statistics). 


15  to  19.5  dB  will  increase  Pp  from  0.5  to  0.8. 


For  the  usual  CFAR  processor,  a  detection  loss  can  be  caused  by 
non-homogeneous  clutter  in  the  auxiliary  clutter  cells  used  to  determine 
the  clutter  statistics  (and  threshold).  In  general,  the  threshold  does 
not  raise  fast  enough  in  the  vicinity  of  a  clutter  discontinuity  so  that 
false  alarms  are  obtained  along  the  clutter  edge.  In  studies  at 
MIT-Lincoln  Laboratory  [6,7],  it  has  been  found  that  the  effects  of 
clutter  discontinuities  can  be  decreased  by  the  use  of  lead/lag  CFAR 
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windows.  These  lead/lag  techniques  sense  the  presence  of  a  clutter 
discontinuity  and  adapt  the  detection  threshold  accordingly.  Such  a 
technique  should  be  even  more  effective  when  combined  with  the  smoothing 
and  reflectivity  reduction  that  might  be  achieved  by  a  wideband 
waveform. 

2.2.2  MOVING  TARGET  DETECTION:  MTI 

Target  motion  information  is  obtained  from  two  general  classes  of 
coherent  radars  -  MTI  and  pulse  Doppler.  Although  the  distinction 
between  these  two  types  is  somewhat  hazy,  the  usual  distinction  is  that 
an  MTI  radar  normally  operates  with  a  pulse  repetition  frequency  (PRF) 
that  is  low  enough  to  avoid  range  ambiguities  (accepting  possible 
velocity  ambiguities)  while  a  pulse  Doppler  radar  operates  with  a  very 
high  PRF  that  avoids  all  target  velocity  ambiguities  (and  nearly  always 
has  range  ambiguities).  The  discussion  of  MTI  systems  will  be  in  the 
general  context  of  an  airborne  system.  A  hypothetical  ground-based 
pulse  Doppler  system  will  then  be  discussed  in  sligthly  less  detail  in 
Section  2.2.3. 

The  Doppler  spectrum  seen  by  a  moving  airborne  MTI  radar  is 
illustrated  in  Figure  2-3.  The  altitude  line  at  the  carrier  frequency 
fQ  can  be  avoided  by  gating  of  the  receiver.  Assuming  the  antenna 
sidelobes  are  uniform,  the  sidelobe  clutter  will  extend  +  2V/X  on  either 
side  of  f  .  The  mainbeam  clutter  spike  is  the  strong  return  from  the 
mainbeam  striking  the  ground.  The  width  of  the  mainbeam  clutter 
spectrum  depends  on  a  number  of  factors  including  antenna  beamwidth, 
aircraft  velocity,  antenna  look  angle  with  respect  to  the  velocity 
vector,  type  and  rate  of  antenna  scanning,  rain  clutter  width  and  wind 
speed  across  the  ground  (causing  an  increase  in  the  width  of  the 
internal  clutter  spectrum).  Because  the  data  is  sampled  at  the  PRF, 
there  are  replicas  of  the  mainbeam  region  separated  from  f  by  integer 
multiples  of  the  PRF. 


bank  of  digital  filters.  The  magnitude  of  the  clutter  return  Is  usually 
much  larger  than  the  magnitude  of  the  target,  return.  Thus  some  type  of 
filtering  must  be  used  to  Increase  the  effective  processing  gain.  This 
is  done  by  means  of  a  clutter  canceller  after  the  signal  Is 
phase-detected.  The  theory  behind  suco  a  canceller  is  that  the 
successive  returns  from  clutter  will  have  the  same  phase  while  the 
successive  returns  from  a  target  will  have  &  phase  shift  caused  by  the 
component  of  target  velocity  in  the  direction  of  the  radar,  Tfus,  for 
example,  a  two-phas;e  canceller  y/hico  subtracts  successive  pukes  will 
have  only  a  very  small  (residual)  output  when  clutter  only  is  present 
and  a  much  larger  output  when  a  moving  target  is  present. 

The  canceller  output  is  then  passed  through  a  Doppler  filter  bank. 
A  system  diagram  is  shown  in  Figure  2-4  where  the  MTI  filter  is  a 
clutter  canceller  and  Doppler  filter.  The  returns  from  N  pulse 
transmissions  are  processed  by  the  FFT  to  provide  N  filter  outputs 
covering  the  Doppler  frequencies  between  zero  and  the  PRF.  The 
signal-to-noise  ratio  in  each  filter  depends  on  the  filter  response  and 
the  clutter  spectrum.  Circuitry  will  offset  the  IF  or  RF  frequency  to 
center  the  clutter  spectrum  at  zero  frequency  so  that  it  can  be  removed 
by  the  pulse  canceller. 

The  output  of  each  of  the  Doppler  filters  is  compared  with  a 
threshold.  If  the  threshold  is  exceeded,  it  is  assumed  that  a  target  is 
present,  and  a  target  will  be  indicated  on  the  display  for  that  range 
and  azimuth  position.  If  rain  clutter  is  present  in  some  of  the 
filters,  these  filters  can  be  desensitized. 

Since  the  radar  will  detect  only  targets  out  of  mainbean  clutter, 
the  radar  must  have  a  long1  antenna  and  operate  at  a  high  frequency  so 
that  the  Doppler  shift  of  slowly  moving  targets  will  be  outside  ma inbeam 
clutter.  (Note  however  that  the  size  of  the  aircraft  carrying  the  radar 


Figure  2  —  4.  A  Coherent  MTI  Implementation 
(with  square— law  detector) 
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will  impose  Important  physical  constraints  on  the  antenna,  prime  power 
available,  etc.)  With  a  high  operating  frequency,  rain  clutter  becomes 
significant. 

The  PRF  Is  selected  so  that  the  unambiguous  range  will  exceed  the 
maximum  surveillance  range.  At  the  high  operating  frequency  of  the 
radar,  blind  speeds  will  then  occur  within  the  ground  target  velocity 
spread  expected.  The  PRF  can  be  staggered  to  raise  the  lowest  blind 
speed.  The  radar  could  transmit  two  bursts  of  pulses  with  different 
PRFs  at  one  beam  position  or  transmit  pulse  bursts  with  different  PRFs 
on  adjacent  beam  positions. 

If  we  consider  only  mainbeam  clutter,  the  signal-to-clutter  ratio 
before  MTI  processing  will  be  the  ratio  of  the  target  cross  section 
to  the  clutter  cross  section  of  a  resolution  cell  given  by 

1  -  gt _  (2-2) 

C  $  oQ  p  sec  6  Rq 

Typical  values  for  these  parameters  are 

at  ~  20m  (target  cross  section  of  a  tank) 

<P  =  0.01  rad  (azimuth  beamwidth  =  0.6°) 

a  =  0.03  (-15  dB,  normalized  clutter  cross  section) 

U  5 

Rq  =  10  m  (range  from  aircraft  to  ground) 
p  =  10m  (range  resolution) 
sec  <5  »  1  (  6  is  depression  angle) 

for  which  S/C  -  0.07  or  -12  dB.  Obviously,  the  S/C  will  be  even  smaller 
if  targets  with  smaller  c?t  are  of  interest. 

Since  a  signal-to-cluter  ratio  of  over  30/1  is  typically  needed  for 
reliable  detection,  the  MTI  processing  must  improve  the 
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With  an  airborne  MTI  radar,  one  attempts  to  center  the  mainbeam 
clutter  spectrum  in  the  notch  at  zero  frequency  in  a  Doppler  filter  bank 
by  shifting  an  IF  or  RF  frequency.  Therefore,  such  a  system  attempts  to 
filter  out  the  mainbeam  clutter  by  attenuating  it  in  the  stop  band  of  a 
filter.  If  the  mainbeam  clutter  spectrum  is  broad  enough,  it  might  also 
be  necessary  to  ignore  the  output  from  the  "lowest  order"  filters  In  the 
filter  bank. 

A  major  limitation  of  a  conventional  MTI  radar  is  that  the  minimum 
detectable  velocity  is  essentially  governed  by  the  antenna  beamwidth. 
Only  targets  outside  of  the  mainbeam  clutter  region  will  be  detected. 
The  equation  for  the  ideal  width  of  the  clutter  spectrum,  Af^  is 

Af  -  2V4>  sin  e  cos  5  ,0  .. 
*  d  "  X  (2-1' 


where 

V  =  aircraft  velocity, 

4>  =  antenna  beamwidth, 

0  *  squint  angle  from  velocity  vector, 

<5  =  depression  angle, 

X  =  wavelength. 

A  narrow  azimuth  beamwidth  minimizes  the  clutter  spread  caused  by 
platform  motion.  In  elevation,  the  beam  must  be  wide  enough  to  cover 
the  surveillance  area,  and  the  pattern  is  shaped  to  provide  coverage 
over  a  range  of  depression  angles. 

Such  a  radar  typically  employs  a  linear  FM  waveform  or  some  other 
form  of  pulse  compression  to  obtain  an  adequate  signal-to-noise  ratio. 
At  each  beam  position,  the  radar  transmits  a  burst  of  pulses,  and  the 
returns  are  digitized  and  processed  using  digital  techniques.  The 
processor  usually  consists  of  a  digital  pulse  canceller  followed  by  a 
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slgnal-to-clutter  ratio  by  about  27  dB  to  detect  the  target.  A  figure 
of  merit,  I,  for  an  MTI  system  Is  the  output  target- to-cl utter  ratio 
divided  by  the  target-to-  clutter  ratio  at  the  Input  averaged  over  all 
target  velocities. 

Depending  on  the  type  of  filter  employed  (two-pulse  canceller, 
three-pulse  canceller,  digital  canceller  which  employs  feed-forward  and 
feed-back),  an  overall  Improvement  of  25  to  35  dB  In  the  target-to- 
clutter-ratlo  can  be  achieved.  Higher  ratios  are  very  difficult  to 
achieve  In  a  conventional  operational  system.  Factors  limiting  improved 
performance  include:  Doppler  filter  straddling  loss,  increase  in 
spectral  width  caused  by  platform  motion  and  antenna  scan  motion, 
multipath  loss,  target  fluctuation  loss,  timing  jitter,  pulse  width 
variation ,  ampl Itude  instability,  second  time  around  clutter, 
limiting  In  the  receiver,  and  finally,  noise  quantization. 

This  means  that  such  an  MTI  radar  has  another  important  limitation 
(In  addition  to  a  minimum  target  velocity)  in  that  enough  processing 
gain  may  not  be  available  to  see  low  cross-section  targets  embedded  in 
harsh  ground  clutter.  These  limitations  are  discussed  briefly  in  the 
following  paragraphs. 

The  characteristics  of  the  Doppler  filter  bank  and  the  clutter 
spectrum  affect  the  clutter  improvement  factor  I.  Figure  2-5a  shows  the 
frequency  response  of  the  Doppler  filter  bank  as  well  as  the  clutter 
spectrum.  Ideally,  the  control  system  would  adjust  the  reference 
frequencies  to  center  the  clutter  spectrum  at  zero  frequency;  however, 
there  will  normally  be  an  error  and  the  clutter  spectrum  will  be  offset 
from  zero.  A  similar  situation  exists  when  a  Doppler  filter  bank 
processes  the  data  (see  Figure  2- 5b ) .  As  with  the  single  filter,  the 
clutter  is  moved  into  the  notch  at  zero  frequency.  Although  weighting 
will  be  applied  to  reduce  filter  sidelobes,  the  sidelobe  level  will 
limit  the  achievable  value  of  I  to  something  less  than  40  dB. 


Amplitude 


Figure  2 -5a.  The  Effect  of  Doppler  Offset  Error  E 


Figure  2-5b.  Offset  Error  with  a  Digital  Filter  Bank 
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The  clutter  power  In  the  antenna  sldelobes  also  limits  I  because 
this  power  will  be  spread  over  Doppler  frequencies  between  ±2V/X.  To 
achieve  theoretical  performance,  there  must  be  adequate  short-term 
stability  In  the  reference  frequencies,  the  timing  circuitry,  and  the 
signal  processing  circuitry;  hence,  system  stability  Is  another  factor 
limiting  I.  The  width  of  the  clutter  spectrum  In  the  return  depends  on 
the  antenna  beamwldth  and  the  aircraft  velocity;  the  narrower  the 
beamwldth  and  the  slower  the  velocity,  the  narrower  the  spectrum.  For 
the  monostatic  radar  to  do  the  job  required,  the  moving  targets  of 
interest  must  lie  outside  the  clutter  in  the  mainbeam  of  the  antenna. 

When  the  targets  of  Interest  are  inside  the  mainbeam  clutter, 
displaced  phase  center  antenna  techniques  (DPCA)  are  necessary  [8], 
These  techniques  require  that  the  phase  center  of  the  antenna  move 
toward  the  rear  of  the  aircraft  to  compensate  for  the  forward  motion  of 
the  aircraft.  Subtracting  two  successive  pulses  greatly  reduces  the 
clutter  and  permits  the  detection  of  targets  with  lower  velocities.  One 
DPCA  technique  uses  an  antenna  with  two  different  phase  centers  which 
can  be  switched  In  alternately  from  pulse  to  pulse.  The  conditions  for 
cancellation  are  that  the  first  pulse  transmitted  from  the  first  phase 
center  and  the  second  pulse  transmitted  from  the  second  phase  center 
line  up  in  the  boresight  direction  and  that  the  beams  are  matched. 
Another  DPCA  technique  employs  a  two-lobe  antenna  pattern  similar  to  a 
monopulse  antenna  with  a  sum  pattern  and  a  difference  pattern,  and  the 
difference  pattern  is  used  to  compensate  for  platform  motions.  Finally, 
it  is  possible  to  combine  SAR  and  DPCA  principles  to  achieve  detection 
of  even  slower  targets  [9'J.  Significant  clutter  reduction  (>80  d B )  may 
be  possible.  It  is  also  important  to  note  that  such  systems  are  very 
complex,  both  in  terms  of  hardware  and  signal  processing. 
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Assuming  the  clutter  is  Gaussian,  we  can  add  the  clutter  power  and 

noise  power  together  to  get  the  total  background  power  which  is  Gaussian 

like  noise.  Standard  curves  for  the  detection  of  power  targets  In  noise 

then  apply.  Let  S/B  be  the  signal-to-background  power  ratio.  Then 

J_  _  J_  ,  1  (?  , 

s/b  "  s/n  Ts/cTT  u 

where  S/N  and  S/C  are  the  signal-to-noise  and  signal-to-cl utter  ratios 
respectively.  (Usually  the  system  is  clutter  limited  and  S/N  can  be 
neglected.)  After  S/B  is  calculated,  the  probability  of  detecting  a 
target  can  be  determined  from  standard  detection  curves  or  from 
equations. 


A  Swirling  case  I  target  model  (scan-to-scan  Rayleigh  fluctuation 
in  amplitude)  is  usually  assumed.  Figure  2-6  has  plots  of  the  proba¬ 
bility  of  detection  versus  the  signal-to-noise  ratio  (signal-to- 
background  ratio  for  the  case  at  hand).  There  are  curves  for  several 
false  alarm  numbers.  The  false  alarm  number  n'  is  the  number  of 
independent  opportunities  for  a  false  alarm  in  the  false  alarm  time. 
The  false  alarm  time  is  the  time  during  which  the  probability  is  0.5 
that  there  will  be  no  false  alarm.  The  false  alarm  probability  Pp^  is 
approximately 

0.693 


PFA  " 


n ' 


(2-4) 


Examples  of  false  alarm  calculations  are  given  in  the  section  discussing 
pulse  Doppler  radars.  After  the  signal-to-background  ratio  has  been 
computed,  the  probability  of  detection  can  be  read  off  of  Figure  2-6  or 
calculated. 


For  each  data  cell,  a  threshold  will  be  established  and  a  detection 
will  occur  whenever  the  threshold  is  exceeded.  The  threshold  settings 
will  differ  with  the  Doppler  filter  because  there  will  be  more  clutter 
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In  the  Doppler  cells  closest  to  the  malnlobe  clutter.  The  probability 
of  detecting  a  Swerling  case  I  target  at  a  given  slant  range  will  depend 
on  the  radial  velocity.  The  filter  response  varies  with  Doppler,  and 
targets  near  zero  velocity  or  the  first  blind  speed  will  not  be 
detectable.  Since  S/N  and  S/C  are  functions  of  the  slant  range,  the 
probability  of  detection  will  vary  with  slant  range  for  a  constant 
radial  velocity. 

It  is  also  possible  to  construct  a  non-coherent  MTI  system  which 
essentially  detects  the  beat  frequency  of  the  moving  target  and  the 
clutter.  Such  a  radar  depends  on  a  strong  clutter  return  being  present 
in  the  target  cell  and  is  not  very  effective  in  detecting  moving  targets 
against  a  low  clutter  background.  The  detection  of  moving  vehicles  on 
roads  might  be  very  difficult  while  the  detection  of  airborne  targets  is 
generally  impossible  with  a  non-coherent  MTI  system.  Another  limitation 
of  such  non-coherent  systems  is  that  improvement  factors  are  3  to  5  dB 
lower  than  is  the  case  with  coherent  MTI  systems.  This  further  reduces 
the  possibility  of  seeing  low  cross-section  moving  targets. 

2.2.3  MOVING  TARGET  DETECTION:  PULSE  DOPPLER 

A  pulse  Doppler  radar  does  not  use  a  clutter  canceller  and, 
instead,  uses  the  fact  that  a  high  PRF  will  separate  the  clutter  and 
moving  target  frequency  spectrums  and  prevent  the  target  Doppler  from 
being  folded  back  into  the  clutter  spectrum.  A  clutter  rejection  filter 
is  used  to  separate  the  clutter  and  target  returns.  The  output  from 
each  range  gate  is  then  passed  through  a  bank  of  Doppler  filters  (which, 
for  modern  systems,  is  synthesized  by  an  FFT  operation).  Each  filter 
output  is  detected.,  non-coherently  Integrated,  and  finally  thresholded 
to  determine  the  presence  of  a  target.  A  very  understandable  discussion 
of  the  performance  of  pulse  Doppler  radars  is  given  by  Meltzer  and 
Thaler  [10].  A  schemetic  diagram  of  such  a  system  is  given  in  Figure 
2-7. 


Figure  2-7.  Block  Diagram  of  a  Pulse  Doppler  Receiver 


To  illustrate  the  uses  and  limitations  of  this  concept,  a  ground- 
based  radar  will  be  hypothesized.  A  typical  next  generation  ground- 
based  SAM  radar  could  have  the  following  characteristics: 


pRP  _  ^15,000  pps ,  seach  mode 
100,000  pps,  track  mode 


center  frequency,  f0 
wavelength,  X 
azimuthal  beamwidth,  BW 
antenna  scan  rate,  SR 


3  GHz, 

0. 1m, 

10° 

180°/sec. 


The  radar  is  designed  to  detect  aircraft  and  helicopters  flying  at 
tree-top  altitude  and  against  a  foliage  background.  Another  requirement 
is  that  such  a  system  be  capable  of  detecting  sensors  on  helicopter  hubs 
when  the  main  body  of  the  helicopter  Is  beneath  the  tree  tops.  It  is 
desired  that  the  system  have  a  range  of  about  35  km. 


One  minor  limitation  of  this  type  of  pulse  Doppler  system  is  that 
it  is  range-ambiguous.  This  can  be  corrected,  however,  by  the  use  of 
multiple  PRF's  that  are  staggered  from  one  burst  of  pulses  to  the  next. 
For  the  hypothesized  system,  the  unambiguous  range,  Ramb»  has  the  values 

,4. 


R 


amb 


C 

2TPRF7 


10  m  (Search  Mode) 
1500m  (Tracking  Mode) 


The  first  blind  speed,  V^,  for  this  radar  occurs  at 

Vb  =  =  750  m/s 


and  all  radial  velocities  smaller  than  this  are  velocity  unambiguous. 
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The  number  of  pulses,  N,  which  illuminate  a  target  during  an 
antenna  scan  is  approximately  given  by 

N  =  |p|  x  PRF  =  830. 


It  is  postulated  that  a  three  PRF  technique  will  be  used  to  resolve 
range  ambiguities  where  the  PRF's  are  related  by  the  ratios  of  closely 
spaced,  relatively  prime,  integers.  This  means  that  approximately  N/3 
of  the  pulses  per  target  scan  will  be  at  each  PRF.  For  this  reason,  it 
will  be  assumed  that  the  number  of  pulses  at  each  PRF,  denoted  by  Np, 
which  illuminate  the  target  during  a  target  scan,  has  the  value 

N  =  256 
P 

This  is  the  number  of  pulses  that  are  coherently  integrated  by  the 
Doppler  filter. 


The  overall  probability  of  detection,  Pp,  per  scan  depends  on  how 
the  probabilities  for  each  PRF  are  combined.  For  example,  if  the  target 
must  be  detected  by  each  of  the  three  PRF  groups,  then  the  overall 
probability  of  detection  per  scan,  Pp,  is  given  by 


Thus,  if  Pp-0.99,  it  follows  that  Pp=0.97. 


The  probability  of  false  alarm,  PpA,  is  the  false  alarm  per 
threshold  test.  For  purposes  of  illustration,  a  value  of  PpA  =  5  x  10“6 
will  be  chosen.  Thu;;,  if  there  are  M  Doppler  filters  tested  per  range 
gate  (resolution  cell)  and  the  filter  outputs  are  independent,  then  it 
follows  that  the  false  alarm  probability,  PpA,  per  range  gate  tested 
(which  is  defined  as  the  probability  that  at  least  one  of  the  filter 
clutter  noise  outputs  for  one  of  the  PRF  groups  exceeds  a  threshold)  is 


W  "  1 
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given  by 

1  -('  -  pfa)3M"  3MPfa 

The  approximation  holds  since  PpA  Is  usually  a  very  small  number. 

For  example,  suppose  the  radar  is  scanning  over  a  complete  360° 
circular  sector  which  corresponds  to  approximately  36  independent  range 
strobes  since  the  azimuth  beamwidth  during  search  is  10°.  A  range 
resolution  of  about  1,500m  In  the  search  mode  Is  assumed.  If  the  range 
from  5  to  35  km  is  being  searched  for  targets,  then  there  are 
approximately  20  cells  or  range  gates  per  strobe  or  about  720 
Independent  cells  per  360°  sector.  Finally,  this  corresponds  to 
approximately 

EFA  =  760  ^FA  '  760  (3M)  PFA 


expected  false  alarms  per  rotation  (i.e.,  every  2  sec). 

As  an  example,  assume  that  velocity  will  be  coarsely  quantized  in 
the  search  mode  of  operation  and  a  value  of 


M  =  number  of  Doppler  filters  per  range  gate  *  3, 


-6 

will  be  chosen  for  illustration.  If  a  value  of  Pp^=5xl0  is  chosen  as 
the  operating  point,  then  the  expected  false  alarm  rate  is 


EpA  =  0.03  expected  false  alarms  every  two  seconds. 

This  corresponds  to  about  1  expected  false  alarm  every  minute.  The 
verification  of  a  detection  before  going  into  a  tracking  lock-on  can  be 
accomplished  by  requiring  a  detection  at  the  "same"  place  on  two 
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consecutive  scans.  This  win  lower  the  overall  detection  probability 
discussed  above  (0.97  per  second)  to  a  value  of  0.94,  which  Is  the 
probability  of  getting  a  detection  on  two  consecutive  scans. 


The  false  alarm  number  associated  with  PFA  15  given  by 


n'  =  =*  1.95  x  105 


As  before,  the  slgnal-to-clutter  ratio  before  MTI  processing  is  given  by 

S  =  °t 
C  <j)a^  R  p  sec  6 

2 

where  a*  *  5m  (helicopter  hub  cross  section) 

<j>  =0.17  rad  (azimuth  beamwidth  =  10°) 
a  =  0.03  (-15  dB,  normalized  clutter  cross  section) 

U  yl 

R  *  10\  (range  to  target) 
p  =  1,500m  (range  resolution  in  search  mode) 

<S  =  0°  (depression  angle) 


This  results  in  a  slgnal-to-clutter  ratio  equal  to 

S/C  =  6.5  x  10"5  or  -42  dB. 

With  a  50  dB  improvement  factor,  this  is  equivalent  to  an  8  dB  signal- 
to-clutter  ratio  after  Doppler  processing;  while  a  60  dB  improvement 
factor  results  in  an  18  dB  signal-to-clutter  ratio  after  processing. 

Using  the  graph  of  Pp  vs  S/N  given  in  Figure  2-6  with  Pp^  =  5x10"®, 
the  result  is 

PD  =  .25,  for  I  =  50  dB, 
and 

Pp  =  .85,  for  I  =  60  dB. 


It  would  be  noted  that  these  are  the  probability  of  detection  for  a 
single  scan  (every  2  sec).  Thus,  for  example,  the  probability  of 
obtaining  a  detection  within  5  scans  (10  sec)  is 

Pq  (5  scans)  =  0.76,  I  =  50  dB, 

PD  (5  scans)  a  0.99,  I  *  60  dB. 

For  future  pulse  Doppler  radars,  the  improvement  factor  is  in  the 
range  of  45  to  55  dB  and  50  dB  is  a  reasonable  operational  number. 
Experimental  systems  are  under  development  to  achieve  an  improvement 
factor  as  large  as  60  dB.  Thus,  the  hypothetical  system  is  pushing  the 
state-of-the-art.  The  limitations  on  the  achievable  value  of  I  are: 
waveform  stability  and  purity  requirements,  sidelobe  clutter,  various 
types  of  spurious  modulation  that  can  appear  on  the  received  signal, 
power  supply  and  line  ripple,  pul se- to-pul se  random  modulation  which 
causes  clutter  spreading,  limiting  in  the  receiver,  and  finally,  noise 
quanti zation. 

The  biggest  limitation  is  the  stability  requirement.  The  reference 
frequencies,  timing  signals  and  signal  processing  circuitry  must  have 
adequate  short-term  stability.  The  most  severe  stability  requirement 
relates  to  the  generation  of  spurious  modulation  sidebands  on  the 
mainbeam  clutter  which  can  appear  as  targets.  Ways  to  reduce  phase 
noise  and  increase  component  stability  are  the  primary  areas  where  pulse 
Doppler  research  is  being  conducted. 

A  number  of  loss  factors  are  present  that  limit  the  improvement 
factors  that  can  be  obtained,  such  as  eclipsing  and  range-gate 
straddling,  multipath  loss,  target  fluctuation  loss,  clutter  power  from 
range  ambiguous  areas,  CFAR  loss,  and  beam  shape  factor. 


Because  the  receiver  Is  gated  off  during  the  transmission  time, 
targets  received  at  this  time  are  also  blanked  or  "eclipsed".  The 
multipath  loss  Is  caused  by  multipath  cancellation  effects  that  are 
present  at  small  antenna  elevation  (or  depression)  angles  when  the 
antenna  is  pointed  close  to  the  horizon  The  loss  caused  by  clutter 
power  obtained  from  range  ambiguous  areas  must  also  be  included.  For 
example,  suppose  the  target  is  located  at  a  range  of  20  km  from  the 
radar  which  Is  twice  the  ambiguous  range  for  the  search  mode.  This 
means  that  during  the  search  mode,  the  target  power  must  compete  with 
the  mainbeam  clutter  coming  from  both  a  range  of  10  km  and  20  km. 

Another  limitation  of  a  pulse  Doppler  system  is  that  monopulse  is 
very  difficult  to  mechanize.  This  is  caused  by  the  problem  of  matching 
multiple  receiver  channels.  These  receivers  each  must  contain  clutter 
rejection  filters  and  must  be  phase-matched  to  permit  proper  tracking. 

2.2.4.  LOW  ANGLE  TRACKING 

A  number  of  techniques  are  available  for  the  angle  tracking  of 
targets.  The  major  catagories  are: 

1)  Amplitude  comparison  monopulse, 

2)  Simultaneous  lobing, 

3)  Sequential  lobing, 

4)  Conical  scan. 

Each  method  produces  a  target  angle  estimate  relative  to  an  antenna 
boresight  position.  This  estimate  is  usually  used  to  drive  antenna 
gimbals  (or  phase  shifters)  so  as  to  reduce  the  error  between  the  target 
direction  and  boresight.  The  first  two  methods  involve  multiple  squinted 
receive  channels  where  the  receive  signal  from  multiple  receive  beams  is 
simultaneously  processed  and  an  angle  estimate  formed.  The  second  two 
methods  are  sequential  in  nature  where  returns  separated  in  time  are 
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processed  from  two  or  more  squinted  beam  positions  and  the  error  angle 
obtained. 

Thus,  the  second  two  methods  are  more  sensitive  to  target 
cross-section  fluctuations  which  might  occur  between  pulses  used  to  form 
the  angle  estimate.  The  radar  has  no  way  of  knowing  whether  a  reduction 
In  received  power  between  pulse  positions  Is  due  to  a  cross-section 
fluctuation  or  Is  due  to  the  target  being  further  off  boreslght. 

A  monopulse  system  has  both  the  theoretical  and  practical  capa¬ 
bility  to  make  more  accurate  measurements  than  the  other  angle  measure¬ 
ment  schemes.  The  variance  of  the  estimated  angle  is  approximately  3  dB 
less  under  conditions  of  equal  Integrated  si gnal-to-noise  ratio  and 
optimal  beam  cross-over  points.  The  monopulse  technique  also  has  an 
Inherent  capability  for  high-precision  angle  measurement  because  its 
feed  structure  is  rigidly  mounted  with  no  moving  parts.  Most  other 
systems  move  the  feed  to  produce  the  required  squinted  beams. 

The  theoretical  standard  deviation  of  all  four  estimates  (for  large 
signal-to-noi se  ratios)  has  the  form 
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where  t  is  a  constant  depending  on  the  type  of  system  (varying  between  1 
and  3);  0g  Is  the  one-way  half  power  bearnwidth,  k$  is  a  constant 
determined  by  the  slope  of  the  antenna  gain  pattern  at  the  beam 
cross-over  points,  and  S/N  is  the  integrated  signal-to-noise  ratio  over 
the  number  of  pulses  used  to  form  the  estimate. 
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All  four  types  of  angle  trackers  are  almost  equally  susceptible  to 
glint  and  multipath.  Here,  the  term  glint  refers  to  the  wander  of  the 
apparent  center  of  reflection  of  the  target  due  to  relative  phasing 
effects  between  the  scattering  centers  making  up  the  target.  All  four 
methods  are  vulnerable  to  angle  capture  which  means  that  the  boresight 
will  tend  to  track  the  largest  scatterer  (if  it  is  significantly 
larger).  The  effect  of  target  glint  is  to  add  an  additional  variance 
factor  to  the  angle  estimate. 

Finally,  all  four  systems  are  vulnerable  to  multipath  which  is 
caused  by  the  receive  signal  being  received  from  two  paths  -  one  direct 
and  one  via  a  bounce  from  the  surface  of  the  earth.  This  also  adds  an 
additional  variance  and  bias  terms  on  the  angle  estimate  which  depends 
on  (1)  beamwidth,  (2)  ground  reflectivity,  and  (3)  relative  sidelobe 
attenuation  between  the  two  paths. 

A  tremendous  number  of  papers  have  been  written  on  the  effect  of 
multipath  and  on  various  schemes  for  reducing  multipath.  For  certain 
ground-based  radars  (such  as  FAA  radars,  for  example),  the  best  schemes 
appear  to  be  a  combination  of  very  careful  radar  placement  along  with 
the  use  of  clutter  fences  and  ground  planes. 

Theoretically,  one  should  be  able  to  use  multiple  receive  ports  to 
sort  out  the  direct  and  reflected  paths.  These  methods  are  based  on 
treating  multipath  as  a  second  radar  target  and  developing  a  two-target 
signal  processor  to  analyze  the  combined  waves  and  sort  the  direct  from 
the  reflected  signal.  The  normal  monopulse  system,  which  is  a 
one-target  processor,  cannot  develop  sufficient  information  to  separate 
the  two  waves  on  a  single  pulse  basis  (monopulse).  References  10-13 
contain  a  description  of  the  major  concepts  and  simulation  results. 
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The  most  modern  method  that  Is  being  tried  to  minimize  multipaths 
Is  the  use  of  very  high  frequencies  (35  or  94  GHz)  for  relatively  short 
range,  ground  based  radars.  The  narrow  beamwldths  that  result  will  very 
efficiently  reduce  multipath  effects  (unfortunately,  such  frequencies 
have  a  large  attenuation  factor  In  rain).  Many  organizations  are 
working  on  millimeter  wave  monopulse  feeds  and  receivers. 

Multipath  is  by  far  the  major  Item  which  limits  monopulse 
performance.  When  it  Is  not  present,  tracking  accuracies  of  1/15  to 
1/20  of  a  beamwidth  are  routinely  accomplished.  The  use  of  a  wide 
bandwidth  waveform  that  Is  linear  or  square-law  detected,  plus 
non-coherent  Integration,  will  reduce  the  effect  of  glint. 

2.3  CONCLUSIONS 

It  was  shown  that  the  detection  of  targets  in  clutter  is  well 
studied  and  therefore  understood  for  contentlonal  radar  designers.  Both 
system  level  constraints  and  fundamental  limitations  were  presented  in 
the  context  of  an  airborne  MTI  and  ground-based  pulsed  Doppler  radar 
design.  Performance  improvements  needed  to  detect  next  generation 
threats  will  require  large  processing  gains  (improvement  factors)  that 
range  from  60  -  80  dB.  These  gains  will  not  come  about  without 
significant  hardware  improvements.  Thus,  it  seems  that  an  investigation 
of  alternate  radar  architectural  paths,  such  as  NSR,  does  seem 
warranted.  The  next  section  develops  a  generalized,  wideband  radar 
model  to  formally  investigate  the  problem  of  target  detection  in 
clutter. 
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3.0 

THEORY  OF  WIDEBAND  RADAR  DETECTION 
3.1  SECTION  SUMMARY 

A  radar  detection  system  is  modeled  with  attention  paid  to  modeling 
all  frequency-dependencies  of  devices,  antennas,  propagation,  and 
scattering.  The  receiver  and  signal  are  selected  optimally  according  to 
the  Ne.yman-Pearson  criterion.  All  operations  are  assumed  linear  and 
time-invariant  (except,  of  course,  the  decision  rule  implementation  in 
the  receiver/processor). 

The  signal  design  problem,  in  this  context,  is  to  choose,  as  it 
turns  out,  the  most  favorable  signal  spectrum  modulus;  this  was  done 
under  the  assumption  of  constrained  signal  energy.  The  actual 
computation  of  the  signal  spectrum  modulus  in  a  particular  instance 
reouires  numerical  procedures;  the  method  of  bisection  was  implemented 
to  do  this. 

Based  on  the  one-shot  (single  pulse)  results  for  a  particular  set 
of  radar  system  parameters  and  models  for  the  scattering  by  object  and 
background,  it  appears  that,  at  modest  but  acceptable  (in  performance) 
energy  levels,  the  optimum  signal  spectrum  modulus  tends  to  have  a 
narrow  support,  near  a  "resonance"  of  the  scatterer.  At  rather  high 
energy  levels,  apparently  well  beyond  that  required  for  satisfactory 
detection  performance,  the  signal  spectrum  modulus  can  be  of  a  broad 
support,  possibly  in  separate  "pieces"  (sets)  corresponding  to 
"resonances"  of  the  scatterer.  However,  this  result  is  strongly 
dependent  on  the  assumed  antenna  and  target  response  characteristics. 
Additional  energy/pulse  duration  constraints  were  included  for  optimum 
design  of  periodic  waveforms,  with  the  same  general  trends  still 
remaining  true.  Finally,  propagation  media  dispersion  was  shown  to  have 
a  strong  effect  on  the  overall  optimum  signal  spectrum. 
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A  clutter/background  model  was  Incorporated  as  derived  from  a 

"two-scale"  electromagnetic  approximation.  Assuming  only  a  "five-scale" 

5 

background,  the  clutter  spectral  density  for  this  sytem  follows  a  k 
law;  this  has  a  dominant  effect  on  the  optimum  signal  spectrum  modulus 
nature. 

A  singularity  expansion  method  (SEM)  based  model  was  used  for  the 
scatterer  wavenumber  -  dependent  scattering.  The  Information  given  by 
the  SEM  method  does  not,  however,  completely  describe  the  scattering 
phenomenon. 

Ir,  Appendices  A  and  B,  it  is  recalled  that:  (1)  complex  exponential 
(sinusoids)  are  the  natural  invariants  of  linear,  time-invariant  systems 
and  hence  the  Fourier  integral  (transform)  representation  of  a  signal  is 
a  naturally  powerful  analysis  tool;  (2)  the  specific  choice  of  a 
"carrier"  frequency  is  completely  arbitrary;  however,  once  chosen,  the 
resulting  complex  envelope  representation  is  then  unique,  under  a 
certain  condition;  (3)  more  generally,  linear,  time-variant  systems 
("operators")  each  have,  generally,  their  own  "generalized  eigenvectors" 
and  appropriate  or  natural  "operational  calculus',  e.g.,  a  certain 
time-variant,  linear,  ordinary  "differential  equation"  operator  is 
naturally  associated  with  the  Walsh  functions. 

3.2  RECEIVED  SIGNAL  MODEL 

We  consider  a  radar  system  geometry  as  sketched  in  Figure  3-1.  It 
is  a  monostatic  configuration,  with  the  "radar  set"  and  "object"  in 
their  mutual  far-fields.  We  will  review  a  model  derivation  for  the 
received  signal  waveform  at  the  antenna  terminals,  given  a  transmitted 
waveform  of  arbitrary  shape.  All  transformations  of  fields  and 
waveforms  are  assumed  linear  and  time-invariant;  therefore,  sinusoids 
are  invariant,  being  changed  only  in  relative  phase  and  magnitude,  upon 
passage  through  the  system  (Appendix  A). 
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Figure  3-1.  Radar  System  Geometry 
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This  received  waveform  has  In  fact  already  been  described  In  [14], 
Ignoring  some  frequency  dependent  aspects  that  we  do  now  wish  to 
Include.  From  page  63,  Equation  3.6  of  [14],  the  received  waveform,  for 
a  transmitted  waveform  exp(-iu)t),  Is  (we  consider  at  this  point  only 
scalar  fields) 

u(t;  oj  )  =  u(o;  o) )  e"^ 


where 

u(o;  w)  = 


(4^-)  JJ dx  dy  exp[i2ky  sin  60]$(x,y;  k)g(x,y;  k)  (3-1) 


and 


[x,  y;  k)  =  exp[ik(x2  +  (y  cos  <S0)Z)/2y‘]  E^(kx/y\  ky/y';  w) 


and  y'  =  RQ  +  y  sin  <Sq,  with  k  =  2ir/X. 

Here  EQ  is  the  (two-dimensional)  Fourier  transform  of  EQ,  an  antenna 
aperture  illumination  and  g  is  the  "reflectivity  density"  characterizing 
the  scene.  For  proper  energy  normalization,//  |E  (x,y)|  dx  dy  =  1. 
The  reflectivity  density  g  is  defined  slmpfy  as  the  multiplicative 
relation  between  the  incident  and  scattered  fields.  This  appears  to  be 
a  sufficiently  robust  concept  and  is  determined  generally  by  the  methods 
of  electromagnetic  scattering  theory. 


In  the  above  development,  the  transmission  medium  is  assumed  to  be 
free  space  and  no  transmitting  and  receiving  antenna  frequency 
dependencies  are  explicitly  given,  though  they  a>"e  implicit  in  EQ's 
depenoence  on  ui .  The  result  will  later  be  extended  to  include 
dispersive  media  in  Section  3.6. 
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We  may  imagine  the  aperture  field  Eq  being  created  formally  as  in 
Figure  3-2.  We  assume  that  a  plane  PP '  is  definable,  across  which  we 
denote  the  radiating  field  as  EQ(F,  t;  to).  When  v(t)  =  exp(-1u>t),  then 
using  the  notation  of  Appendix  A 

EQ(r,  t;  to)  -  EQ(r,  o;  to)  e_1lot,  re  PP'.  (3-2) 

In  order  to  proceed,  we  make  the  crude  but  useful  approximation 

EQ(r,  o;  to)  ■  a t ( co )  EQ(r).  (3-3) 

That  is,  we  assume  that  factoring  the  frequency  and  spatial  dependencies 
is  possible.  Conceptually,  the  analysis  is  easier  but  at  very  large 
relative  bandwidths,  this  may  be  untenable.  Similarly,  upon  reception, 
the  received  scattered  field  will  produce  a  voltage  at  the  antenna 
terminals,  the  mapping  again  depending  on  a  frequency-dependent  transfer 
function,  which  may  be  denoted  ( co ) .  We  incorporate  these  frequency 
dependencies  by  setting 

-iw(t-2R  /c) 

u(t;o  )  =  a(w)  u(o;  to)  p  ,  (3-4) 

where 

a(to)  =  ar(a>)  a ^. ( (jo ) . 

With  such  a  model,  in  any  specific  instance,  the  quantities  EQ,  g, 
ar,  and  at  must  be  determined  in  some  manner. 

An  alternative,  more  rigorous  description  has  been  derived  and  is 
presented  in  Appendix  C.  Instead  of  assuming  an  aperture  is  definable, 
the  far  field  is  described  in  terms  of  vector  moments  of  electric  and 
magnetic  antenna  currents  distributed  over  a  volume,  which  led  to  the 
notions  of  a  "pattern  factor"  and  a  "vector  effective  height"  for 
describing  the  transmitting  and  receiving  transfer  functions  of  an 
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Figure  3—2.  Antenna  Aperture  Field  Generation  Model 
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antenna.  Furthermore,  rather  than  employ  a  reflectivity  density,  a 
(dyadic)  scattering  matrix  is  used  in  the  derivation.  It,  too,  is  a 
linear  pointw  se  relation  between  incident  and  scattered  rields. 
Finally,  the  formulation  in  Appendix  C  is  in  vector  form,  while  the 
results  of  this  section  represent  a  single  scalar  component.  This  more 
detailed  model  of  antenna  reception  and  transmission  transfer  functions 
is  developed  by  considering  Thevenin  equivalent,  circuits  for  the  source 
generator  and  load  impedance.  Implicit  in  the  model,  nonetheless,  are 
the  same  types  of  calculations  required  for  specific  antennas  as  used 
herein.  It  is  probably  fair  to  say  that  the  "quasi-optical "  aperture 
field  approach  will  prove  more  directly  applicable  at  higher  frequencies 
for  aperture  antennas,  while  the  "equivalent  circuit/  antenna  current" 
approach  may  be  more  amenable  in  describing  low-frequency,  wire-like 
structures.  Further  research  is  required  to  develop  insight  into  the 
differences  between  the  two  modeling  approaches. 

3.2.1  NON-SINUSOIDAL  TRANSMITTED  WAVEFORM 

If  an  arbitrary  waveform  f(t),  with  Fourier  integral/transform 
representation  +00 


-oo 


is  actually  transmitted  Instead  of  a  complex  exponential,  then,  by  the 
linearity  assumed  everywhere,  the  corresponding  received  "non- 
sinusoidal"  waveform  is  given  simply  by  the  superposition  principle, 
namely 


Jm 


where  c  Is  the  free  space  speed  of  light.  (Note  that  this  is  not  the 
conventional  Inverse  Fourier  transform  used  In  communication  theory, 
which  employs  the  "kernel"  exp(1u>t),  not  exp(-lwt).)  Similarly,  If  F(t) 
Is  a  "complex  envelope  representation",  corresponding  to  a  real 
waveform 

(  iwrt  ( 

f(t)  =  Re  \  F(t)  e  c  f  , 


where  <*>  is  a  "carrier"  frequency,  the  corresponding  real  received 


waveform  is 


;(t)  =  Re  |  u(t;  u)  e  c  j 


with  the  complex  envelope  representation 


s(t) 


■H 


doj  e 


iu(t-2R  /c) 


{a(w)f  (co)u(o;  uj)  > 


(3-7) 


-(w  +  u>c) 


3.2.2  SCATTERING  REFLECTIVITY  MODELS 


We  suppose  that  g  is  the  superposition  of  two  parts,  a  gQb 
describing  a  scattering  object  whose  presence  (or  not)  is  to  be  decided, 
and  a  background,  or  clutter,  described  by  g  Correspondingly, 


u(o;  w)  =  uQD(o;  w)  +  ucl  (o;  u>). 


(3-8) 


3.2.2. 1  CHARACTERIZATION  OF  TARGET  OBJECT 


As  is  usually  true,  and  assumed  herein,  the  spatial  extent  of  the 
object  is  sufficiently  restricted  so  that 
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where  r'  Is  the  distance  from  the  radar  to  the  point  r'  on  the  object's 
surface. 


Recalling  the  normalization  of  E.  just  mentioned, 

c 


5,(0.  o)  = 


// 


En(x.y)  dx  dy 


antenna 

aperture 
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dx 
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where  jd  is  the  effective  aperture  area. 
Therefore 


W7el2kr  90b '  *  k)  d?'-  (3-10) 
Sob 

It  may  be  convenient  at  this  DOint  to  mention  the  connection  with 
the  radar  scattering  cross-section  (RCS),  normally  defined  for 
historical  reasons,  for  a  scattering  object.  As  in  many  other  radar 
systems,  the  RCS  alone  does  not  sufficiently  characterize  a  scatterer. 
Here,  we  shall  f~nd  that  more  information  is  required  to  construct  the 
optimal  receiver,  specifically,  the  matched  filter.  However,  it  does 
turn  out  that  the  optimal  signal  and  the  system  detection  performance 
depend  only  on  the  modulus-square  of  the  (frequency-dependent)  scattered 
field.  In  this  case,  a  direct  relation  to  the  RCS  is  possible  and 
useful  because  of  the  empirical  data  available  concerning  this 
parameter. 
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The  RCS  Is  defined  to  be 


a  =  4irR0 


(3-11) 


where  Ef  Is  the  far-field  Incicent  on  the  object  and  E$  the  received 
scattered  field.  It  Is  easily  seen  from  the  calculations  summarized 
above  that 


(o,o) 


9 ob(r'i  k)  dr' 


(3-12) 


and  hence  that  the  RCS 
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k)  dr 1 |  • 


(3-13) 


As  Inferred,  we  shall  refer  to  this  form  later. 

Now  suppose  that  the  object's  surface,  SQb,  Is  defined  by  the 
relation 

Fs(x,y,z)  =  z  -  fs(x,y)  -  0,  (x.y)  e  Pcb, 


where  P  ^  is  the  projection  of  SQ^  onto  the  (x,y)-plane  (see  Figure 
3-3). 


Then  completing  the  x~ integration  in  (3-lo),  we  obtain 
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uob(o;  co)  =  1 4irR 
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i2ky  sin  5, 


r  (y;  k)  dy 


(3-14) 


where 


r(y;  k)  = 


geff(x,  y;  k)  sec  y(x,  y). 


2 

If,  over  a  range  of  k  of  Interest,  (1k/4irRo)  and  r  do  not  vary 
significantly  from  the  value  at,  say,  kQ,  then  the  relation  for 
%b^0;  ^  first  above  is  a  Fourier  transform  relation  and  we  have 


U0b(°l  t)  - 


4).^)r(T 


ct 

sir 


-k. 


(3-15) 


However,  this  form  is  not  of  Immediate  Interest  here,  as  we  definitely 
want  to  observe  the  effect  of  the  k-dependtnce  of  r  and  g. 


The  singularity  expansion  method  (SEM)  [15,16]  can  be  of  some 
assistance  In  characterizing  the  scattering  from  objects  of  interest 
here.  It  has  been  asserted,  on  the  basis  first  of  empirial  evidence  and 
later  some  theoretical  support,  that  a  class  of  scatterers  can  be 
characterized  by  their  set  of  "poles"  associated  with  their  transient 
response.  This  is  a  useful  categorization  method  because  these  poles 
appear  to  be  rather  independent  of  aspect  angle. 

Viewing  the  transient  response's  transform  as  a  partial  fraction 
expansion,  the  poles  are  explicit;  however,  the  weights  given  to  each 
such  contribution  (term)  determine  the  zeros  of  this  transform  and  these 
weights,  hence  zero?,  apparently  generally  do  depend  upon  ?spect  angle. 

Since  here  we  are  interested  in  the  wave-number,  or  frequency, 
response  of  a  scatterer,  we  must,  know  both  pole  and  i.ero  locations; 
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therefore,  knowledge  of  the  poles  alone  provides  Insufficient 
Information  for  us  here.  However,  for  simplicity,  we  will  later  assume 
a  specific  pole/zero  pattern  In  order  to  obtain  a  numerical  frequency 
dependence  for  a  target's  cross  section. 

3. 2. 2. 2  CHARACTERIZATION  OF  CLUTTER 

In  order  to  quantify  the  effect  of  clutter,  we  must  examine  the 
clutter  contribution  to  u(o;  w),  namely 

/  it  \2  rr  sin  6n 

ucl(o;  w)  ■  ^4^-j  yy  dx  dy  e  $(x,y;  k)gc1(x,y;  k).  (3-16) 

There  are  two  major  problems  Involved.  First,  the  clutter  reflectivity, 
gcl  must  be  characterized  using  electromagnetic  scattering  theory,  over 
a  very  broa>d  range  of  k,  for  interesting  background  scenes.  (This  Is 
viewed  as  fundamental ly  Important.)  Second,  there  is  a  requirement  to 
clarify  the  effect  of  antenna  pattern  variation,  via  $  ,  over  a  very 
large  range  of  k.  (This  is  considered  a  tedious  analysis  problem.) 

To  see  how  the  often-employed  "convolution  model"  is  obtained  is 
initially  helpful,  though  it  is  not  sufficiently  general  for  us  here. 

If  the  spectrum  of  the  transmitted  waveform  is  sufficiently  narrow,  say 
about  kQ,  then,  in  uc-|(o;  to  )  above,  we  may  replace  k  by  kQ  everywhere 
except  in  the  phase  2ky  sin  <5  .  We  then  get  the  Fourier  transform 
relation 


Ik. 


Ct 


2  sin  5. 


(3-17) 
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where 

+00 

G(y';  kQ)  -  J  dx'  rf(x'.y*  ;  kQ)gcl  (2  sfrhs^  k)’ 

-00 

Therefore, 

/Ik 

ucl(0i  t)  =  sir's"  G  (2  s(£  s0;  ko)-  (3-18) 

Using  (3-18)  and  (3-7),  we  now  have  a  simple  linear  system 
(convolutional)  model  as  Illustrated  In  Figure  3-4. 

In  order  that  G  be  a  sample  function  from  a  wlde-sense  stationary 
random  process,  a  very  convenient  property  for  subsequent  analysis,  we 
require  the  following  correlation  (covariance)  to  be  a  function  of  (tj, 
t2)  only  as  |tj  -  t? |: 

E[s(y';  k0)G*(y”;  k„)]  -  JJ dx'  dx"  £(x',y'  ;  k„) 

g(x", y";  k0)Rg(x'-x”,  y’-y”;  kQ)  (3-19) 

p  +■ 

where  we  have  made  a  simplification  of  variables  y  =  ■= — -----  . 

r  J  2  sin  6Q 

Here  we  assumed  the  clutter  field  9C](X'»  y* 1  kQ)  is  homogeneous  with 
covariance  function  R^. 

Suppose  that  jgc-|}is,  for  wavenumbers  of  interest,  "white",  that 
is. 


r)n  2  sin  6 

R  (•)  =-9 - - — —  5 (x 1  -x " ,  trt2) 
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(3-20) 


Transmit/Receive 
Frequency  Dependence 


Waveform 

Spectrum 


Figure  3-4.  Linear  Model  for  Clutter  Reception 


then 


E[S(-).  G*( -)] 


2  sin  s„  f 

-^7dx 


i 


|#(x',  ctx/2  sin  50;  kQ)|2  fiftj-tg)  (3-21) 


Recalling  that 


l«(x.y;  k0)|2  - 


(3-22) 


there  Is,  of  course,  now  no  phase  of  $  appearing  and  so  this  form 
essentially  plays  the  role  of  a  "window"  (or  truncation)  on  the  sample 
function  of  the  received  clutter  process.  As  usual,  this  window  Is 
broad  relative  to  the  time  extent  of  the  transmitted  pulse  and  certainly 
the  time-extent  of  the  object,  so  we  may  Ignore  Its  windowing  effect. 
(Of  course,  the  associated  gain  is  not  neglected.)  That  is  to  say, 

p 

ignoring  the  y  and  t  dependence  of  \$  \  ,  we  have,  as  the  spectral 
density  of  ucl(o;  t)  , 


M 

ucl 


(x»o  j  kg) 


Vko> 


(3-23) 


where  we  have  indicated  that  Hg  could  have  wavenumber  dependence. 

Examining  the  integral  in  this  last  expression,  and  using  (3-1)  and 
(3-22),  we  obtain 


4 


=  T  /*  ko)  i4  (3"24) 
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If  we  assume,  reasonably,  a  uniform  aperture  Illumination,  then  doing 
the  Integrals,  one  finds 

/d*|S(x,Oi  kQ)|2  =^(^*)d2  (3-25) 

where  D  and  D  are  the  antenna  aperture's  horizontal  and  vertical 
x  y 

aperture  dimensions. 

Now  we  may  write,  on  combining  terms. 


(3-26) 


Again,  It  Is  convenient  to  relate  the  spectral  density  level  of  the 
"white"  clutter  process'  reflectivity  density  to  the  "historically 
defined"  RCS  density,  aQ.  The  only  time  this  is  possible  is  when  the 
clutter  field  (or  process)  is  indeed  white  --  for  only  then,  can  a 
single  parameter  yield  a  useful  characterization.  It  may  be  shown  that 


Vko>  ■  7?  °o  <ko>. 
ko 


There  is  available  the  "two-scale  model"  for  rough  surfaces 
composed  of  a  slowly  undulating  "large-scale"  fluctuations,  treatable  by 
physical  optics,  and  a  more  rapidly  fluctuating  but  small  "fine-scale" 
variation,  treatable  by  a  "boundary  perturbation."  (All  relative  scales 
are  with  respect  to  the  RF  wavelength  X  ).  Assuming  that  there  is  in 
fact  no  large  scale  present,  one  has 


Terjm 


9ci (x»  y;  k)  ■  2  cos  60  [1  -  12k  cos  6Q  £  (x,  y)].  (3-27) 
where  £(x,  y)  Is  the  scattering  surface  variation. 

This  model  applies  under  reasonable  conditions.  In  particular,  at 
Intermediate  Incident  angles  where  there  will  be  no  "specular"  (only 
"Bragg")  scattering  the  "1"  term  may  be  neglected.  So  that 


9cl ( x,  y;  k)«  -1(2  cos  6Q)*  k  S(x,  y) 


(3-28) 


and  now 


n  (k)  =  (2  cos  5q)4  k2  n5- 


(3-29) 


Equivalently,  the  RCS  density 


I,2  (2  COS  6  r  A 

CT0(k)  =  TH  ng(k)  = - <R -  k  n 


(3-30) 


where  is  the  spectral  density  of  the  scattering  surface  perturbation. 

We  may  expect  such  a  model  to  apply  for  rather  flat  terrain  with 
small  scale  variations  due  to  rocks/gravel  or  grass,  e.g.,  to  fre¬ 
quencies  up  to  somewhere  below  X-band.  In  fact,  it  is  observed 

empirically  that  over  such  types  of  terrain,  a  (k)  does  indeed  have  a 
4  0 

k  -dependence  on  wavelength.  This  does  allow  a  means  of  estimating  the 
fundamental  parameter  ,  by  using  the  above  relation  and  empirical 
measurements  of  of  clutter  at  x-band.  See  Table  3-1  [26]. 
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TABLE  3-1.  ESTIMATES  OF  TERRAIN  SPECTRAL  DENSITIES 


Terrain  Type 
"dry,  smooth  desert" 
"golf  fairway" 


"Typical" 

0o 


Estimate  of 

_ 2£_ 

1.633  x  10" 12 
1.633  x  10"11 


When  the  reflectivity  density  gc^  depends  on  k  only  as  a  factor, 
the  above  argument  may  be  reviewed  and  it  may  be  seen  that  the  only 
assumption  necessary  to  obtain  the  final  result  is  a  neglect  of  the 
windowing  effect  of  the  antenna  pattern  in  y. 


Combining  the  above  discussions,  the  model  for  the  spectral  density 
of  the  clutter  return  yields 


ucl 


(4  7 

cos  *0  °A,  ) 

6.  sin  V"/'"*} 

A  ' 


(3-31) 


It  is  interesting  to  recall  how  the  k"  dependence  comes  about  in  the 
clutter  spectral  density.  First,  there  is  a  k^  term  because  of  the  k4 

A  o 

propagation  factor  ( k/ 4 ttR  )  ,  decreased  to  x  by  integration  over  the 

O  p 

x-direction  of  the  far-fleld  antenna  pattern.  There  is  an  additional  k 
because  of  the  small-scale,  rough  surface  scattering,  consistently  with 

4 

a  k  RCS  density,  aQ  dependency. 

As  an  example,  suppose  that  aQ  =  tt/4,  Dx  =  1  m,  *  1/3  m, 

Rq  =  10  Km,  =  1.633  x  10-12; 

then 


&  (u)  =  1.0210*10" 18*  (£)? 

ucl 
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(It  should  be  noted  that  this  value  of  was  arrived  a t  using  a  k  - 
dependent  relation  between  and  empirically-derived  aQ,  thus,  any 
error  In  aQ  could  be  magnified  at  this  point.) 

We  have  therefore  established  the  model  as  described  In  Eqs.  3-10 
and  3-31  and  illustrated  in  the  block  diagram  shown  in  Figure  3-5. 

The  major  limiting  assumptions  made  in  the  course  of  the 
development  are  that:  First,  the  "windowing",  in  the  range  dimension, 
of  the  antenna  pattern  is  long  relative  to  the  pulse  length  and  hence 
practically  ignorable,  and  second,  the  "fine-scale"  background  is 
"white"  over  the  bandwidth  of  interest.  Although  other  assumptions 
could  have  been  made,  it  is  felt  that  the  model  presented  by  Figure  3-5, 
along  with  the  supporting  equations,  is  a  reasonable  first-cut  for  a 
broadband  radar  system  model. 

3.3  OPTIMUM  PROCESSOR  STRUCTURE  AND  RESULTING  PERFORMANCE 

In  order  to  formulate  a  detection  decision  problem,  we,  first  of 
all,  acknowledge  that  s(t)  will  always  be  received  along  with  thermal 
noise  so  the  data  actually  available  for  processing  is 

z(t)  "  s(t)  +  n(t) ,  teT  (3-32) 

where  T  is  some  reception  interval,  possibly  delimited  by  the  antenna 
pattern.  The  thermal  noise  n(t),  teT  is  a  Gaussian,  zero  mean, 
possibly  white  process.  (However,  the  white  assumption  is  not  necessary 


; aTOITkTU  JfVWWVW TCVW. r/.n^’Twci#  f'la  Tx,-n»ffTHirnv*vwB\i 


* 


r, 

» 


v, 


mm 


for  us  to  proceed  here.)  Also,  It  Is  assumed  that  the  clutter-induced 
u„i(o,  t)  Is,  over  T,  a  wide-sense  stationary  random  process  of  zero 
»an,  and  that  It  Is  Gaussian.* 

We  have  the  following  decision  problem:  given  the  data  z( t), 

teT  ,  we  must  decide  whether  z(t)  Is  due  to  clutter-and-noise  or  due  to 

signal  (object)  plus  cl  utter-and-noise.  Because  of  the  assumDtlon  of 

normal  distributions,  It  Is  well  known  that  under,  say,  the 

Neyman- Pearson  criterion,  the  optimum  decision  device  Is  a  "matched 

filter"  followed  by  a  "sampler",  the  sample  being  compared  to  a 

threshold  [21,22],  The  performance  Is  uniquely  determined  by  a 

o 

signal-to-noise  ratio  (SNR)  parameter,  denoted  as  "d  "  in  the  following. 
2 

Here  d  is  the  normalized  "distance"  between  the  ^wo  hypothesized 
signals. 

Without  loss  of  generality,  we  may  allow  a  "pre-receiver11  with 
transfer  function  H(u).  Also,  we  should  note  that  the  predominant 
thermal  noise  at  these  frequencies  is  likely  due  to  that  radiated  by  the 
viewed  terrain.  Thus,  one  arrives  at  the  sketched  system  model  shown  in 
Figure  3-6. 

To  abbreviate  notation  at  this  point,  set 

B(t)  =  uQb(o;  t),  C(t)  *  u c  1  ( o ;  t). 


Also,  denote  as  s(t)  and  n(t)  the  resultant  signal,  due  to  B,  and 
thermal  (effective)  noise,  due  to  C,  n^er,  and  nrcv,  respectively. 


*TT  we  do  not  assume  Gaussian  distributions  and  then  adopt  the  softer 
maximum  signal-to-noise  ratio  (SNR)  criterion,  we  are  led  to  precisely 
the  same  results. 


60 


It  is  well  known  [21]  that  the  optimum  (“matched")  filter  is  given 


K(«)  •  -§^ 


(BfaH) 


(3~33) 


(■V,)l^1  *\,)151  ^|Mr 


where  &  (u>)  Is  the  power  spectrum  to  the  total  received  noise, 
n 


The  detection  performance  is  completely  characterized  by  the  parameter 


'  >'■/  : 


_ .-UMlisc _ ,3 ... 

^  l?r|2^n  +  ^|fa|2]|H|2  <3  34> 

ter  "rev  J 


.  **v  .9  2 

We  notice  |H|  multiplies  all  and  hence  cancels  out  of  the  d  expression 


where  it  is  not  zero.  Thus,  we  ignore  it  in  what  follows,  except  we  may 


wish  to  restrict  f(u>)  to  some  w-set,  e.g.,  for  practical  reasons. 


3.4  OPTIMUM  SIGNAL  DESIGN 


We  now  have  arrived  at  a  position  where  one  can  address  the  problem 
of  special  interest  here:  "How  should  the  signal  spectrum  f(w)  be 
chosen?" 


We  observe  immediately  that  the  only  dependence  of  dc  on  f  is 

i  ■*w/  \  1 2 

through  |fU))|  ,  that  is,  in  choosing  f  to  maximize  performance  that  is 


2  ^  2 
done  by  maximizing  d  ,  we  are  only  concerned  with  the  choice  of  |f(u>)|  . 


As  is  well  known  in  radar  work,  this  may  leave  freedom  to  choose  f  to 
meet  other  design  restraints,  e.g.,  Doppler  "resolution." 
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We  will  choose  f  under  an  energy  constraint,  as  will  be  true  in 
practice  and  as  a  constraint  is  necessary  to  avoid  impractical 
solutions.  We,  therefore,  have  the  following  problem: 


max 


,  subject  to  the  constraint: 


E 


f 


The  solution  to  this  problem  is  given  on  page  95  of  [14]:  for  all 


permitted  n> , 


(3-35) 


The  Lagrange  multiplier  parameter  is  to  be  chosen  so  as  to  meet  the 
energy  restraint.  This  apparently  must  be  done  by  an  approximation 
technique.  Here,  we  use  an  algorithm  method  based  on  the  technique  of 
"bisection."  Given  appropriate  choices  of  lower  and  upper  bounds  on  the 
frequency  range  of  interest,  the  algorithm  always  converges  because,  as 
may  be  observed,  the  energy  increases  monoton ically  (but  not  strictly 
so)  with  increasing  frequency  interval  length. 

This  solution  has  a  nature  especially  interesting  here:  As  the 

.  ,'*v  O 

energy  constraint  is  varied,  the  ^-support  of  [  f (uj ) [  can  dramatical  ly 
change,  e.g.,  from  "narrow-band"  to  "broad-band".  The  variation  is  much 
more  complex  than  a  simple  scale  change.  This  is  demonstrated  in 
Section  3.4.1. 

3.4.1  NUMERICAL  EVALUATION 


For  the  calculations  to  follow,  we  suppose  that  iTt(w)  =  a  ( w)  and 

v  r 

that  the  terrain  thermal  noise,  with  spectral  density  nQ,  dominates  the 
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receiver  thermal  noise.  Using  such  assumptions,  (3-36/  can  be  written  as 


jf  (oj)  !2  -  max  (  0 ,  j \ 


a 


^  |a!;’|3i2|a-|nt 


/ 


-yia  h. 


(3-36) 


In  addition,  it  will  be  numerically  cgirl/enient  to  consider  all  terms  as 
a  function  of  the  frequency  v,  expressed  in  Megaherz.  That  is,  we  use 


k  *  ~  *  2.0944  x  10“2  v  (v  in  MHz) 


(3-37) 


On  substituting  (3-37),  and  changing  variables  in  (3-36),  we  obtain 


f ( v) |  =  max  7  0, 


a(v)|v5  L 


c0v>/|a(v)|o(v)  - 


(3-38) 


where  cQ  is  chosen  to  meet  the  energy  constraint  and 

.3 

sin  6  R  D  n 

A  =  -  o  -2*  x  1.559  x  10  . 

1  ros4 «  jtZ  \ 


(3-39) 


For  completeness,  the  same  assumptions  yield  an  energy  expression 

Hf=106/  |f(v)|2dv  (3-40) 


and 


I2  =  ^y’)f(v)l2la(v)i2  v2.tfv)  dv 


+  |f'(v)Ha(v)  |2  v5/A1 


(3-41 ) 


where 


4  U 

A-  =  f - -T—  (2.0944  x  10’2)' 

H-R0)4n0 


=  1.759  x  10 


-2  jd2 


Rono 
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3.4. 1.1  TARGET  CROSS  SECTION  MODEL 

Of  course,  the  specific  frequency  dependence  of  the  target  cross 
section  will  have  a  strong  effect  on  the  optimum  signal  spectrum.  We 
model  the  cross  section  dependence  by 

a(v)  =  om  W(v)  (3-42) 


where  max  {W(v) }  =  1. 
v 

Moreover,  we  assume  a  specific  W(v)  dependence  as  sketched  in  Figure 
3-7.  It  Is  motivated  by  the  pole  pattern  for  a  MIG  19  obtained  by  the 
SEM  method  [15,16].  (In  later  computations,  we  will  approximate  the 
response  by  a  suitably  chosen  "Butterworth"  filter  characteristic.) 

3. 4. 1.2  ANTENNA  FREQUENCY  DEPENDENCE 

It  Is  also  obvious  that  the  antenna  response  will  strongly  drive 
the  signal  design.  For  the  numerical  examples  to  follow,  we  use  the 
transfer  function  for  a  reslstlvely  loaded  horn,  as  reported  by  Kanda 
[17]  and  sketched  in  Figure  3-8.  It,  too,  was  modelled  by  a  suitable 
Butterworth  response  In  the  analysis  to  follow. 

3. 4. 1.3  SPECIFIC  NUMERICAL  CALCULATIONS 

We  begin  by  taking  Aj  =  10®  and  A^  *  0  ^2  =  ®  x  10®  and  suppose 
that,  for  whatever  practical  reasons,  the  signals  spectral  support  is 
restricted  to 


v  •  =  135  and  vmav  =  175  (MHz), 

min  max 
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V(MHz) 


Figure  3-8.  Assumed  Antenna  Frequency  Response 


Jm _ _ _ 

Recall  that  from  Figure  3-7,  this  frequency  interval  contains  the  higher 
frequency  resonance  of  the  subject  target. 

The  resulting  normalized  optimum  signal  spectrum*  is  presented  in 

2  4 

Figure  3-9  for  energy  constraints  of  1,  10  and  10  Joules, 

respecti vely.  In  addition,  the  figure  presents  the  resulting  SNR 
2 

parameter,  d  .  It  may  be  observed  that  for  E*.  =  1,  the  resultant 

|f(o>)  I  is  quite  narrow  relative  to  the  permitted  bandwidth  and  that  the 

2 

system  performance  would  be  quite  good,  since  d  is  roughly  60  dB. 

Increasing  the  allowable  energy  does  result  in  a  signal  which  fills  the 

available  bandwidth  and  produces  an  even  larger  d  . 

Additional  computations  were  made  to  examine  the  effects  of  the 

frequency  dependence  of  the  clutter  spectrum  S  ( k ) .  (Recall  that  the 

5  c 

original  dependence  was  shown  to  be  a  k  dependence.)  A  clutter 

3  0 

spectrum  dependence  of  k  and  k  was  used  with  a  wide  variety  of  signal 

energy  constraints.  In  all  cases,  the  signal  spectrum  was  quite  narrow 

and  centered  at  150  MHz.  Moreover,  for  fixed  energy,  the  system 

2 

performance  measure  d  increased  as  the  clutter  "power  law"  decreased, 

4  2 

as  would  be  expected.  For  example,  with  =  10  ,  d  went  from  90  dB  to 
90.3  dB  as  the  clutter  spectrum  went  from  a  cubic  to  a  zero  exponential 
dependence. 


It  seemed  that  little  additional  frequency  dependence  would  exist 
without  including  a  broader  allowable  interval  to  include  the  additional 
target  resonance.  Thus,  we  started  with 


vm.  =  35  and 
mi  n 


v 

max 


200 


and  again  assumed  a  fifth  law  clutter  spectrum. 


~  *The?  opt i mum? spectrum  is  plotted  in  a  normalized  way  such  that 
|fQ(a))|  /max|f0|  is  unity. 
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Figure  3-9.  Effect  of  Energy  Constraint  on  Particular  Optimum 


Signal  Design  ( ''P*ina,*135,'y(Bax*i75) 
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Allowing  for  a  large  signal  energy,  =  10®  Joules,  we  obtained 
the  signal  spectrum  presented  In  Figure  3- 10a.  Note  the  signal  Is  now 
matched  to  the  target's  lower  resonance  (compare  to  Figure  3*-9c). 

•  .9 

We  would  guess  that  |f(co)|  should  generally  broaden  as  A^  Is 
decreased  because  Xj  must  be  larger  to  meet  the  energy  constraint.  This 
can  be  easily  accomplished  by  Increasing  the  clutter  reflectivity  by  a 
factor  cf  ten  (see  Eq„  3-39).  The  resultant  signal  spectra  is  shown  in 
Figure  3- 10b;  indeed  It  has  broadened  as  expected. 

The  constant  Aj  may  be  decreased  further,  and  Ai,  increased  by 
decreasing  the  range  R  to  1  Km  (see  Eqs.  3-39  and  3-41).  In  that 

in 

instance  Aj  becomes  10"  and  A^  is  increased  to  5  x  10  .  For  the 
constraint  of  F^  =  10®,  the  resultant  signal  is  the  same  as  that  plotted 
in  Figure  3  10b.  If,  however,  we  use  E^  =  100,  the  jf^)!^  has  perhaps 
narrowed  some,  as  shown  in  Figure  3-10c.  The  narrowing  v>  even  more 
pronounced  for  E^  =  10,  as  snown  in  Figure  3- 11a.  Note  also  there  Is  a 
slight  signal  component  centered  at  the  other  object  resonance.  If  the 
energy  is  further  decreased  to  E^  =  0.01,  then  rather  surprisingly,  the 
support  of  |?’(<d)|  jumps  to  a  location  around  the  larger  resonance  of 
the  object,  as  shown  in  Figure  3- lib.  As  E,  is  further  decreased  to 
10"  ,  a  relative  narrowing  of  the  |f(w)|*  again  occurs. 

The  above  results  were  obtained  with  Aj  =  ?0?  and  A^  =  5  x  10i0. 
If  we  again  return  to  a  longer  range  of  Rft  =  10  Km  and  a  smoother 
terrain  oackground  so  that  Aj  =  10  and  Ai>  =  5  x  10  and  use  Ef  =  10  , 
the  resulting  sharply  narrowed  spectrum  is  obtained,  as  shown  in  Figure 
3-12. 
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Figure  3—10.  Optimum  Signal  for  Larger  3andwidth  Interval 
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Figure  3-12.  Optimum  Signal  With  Reduced  Energy  Constraints 
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3.5  MULTIPLE  PULSE  (PERIODIC)  SIGNAL  DESIGN 

It  is  frequently  true  in  radar  systems  that  the  waveform  modulation 
transmitted  should  be  pulsed.  This  nay  be  because  an  area  is  to  be 
scanned  or  that  peak  power  and  resolution  requirements  are  such  that 
adequate  detection  and  parameter  measurement  performance  can  only  be  met 
by  utilizing  a  number  of  pulse  returns. 

Thus  we  assume  now 

N-l 

f(t)  -  £  f  (t  -  nT),  (3-43) 

n=0 

a  repetition,  with  period  T$t  of. a  transmitted  waveform  fQ.  We  see  that 
the  waveform  spectrum  (Fourier  transform) 

?(<*»)  =  Q(w)  f0(u>) 


where 


^  N-l  -inT  w 

QM  =  E  e  5 

n=0 


-1(N-l)Tsu/2  jsin(NTsu/2)  j 

|TTFTT^r)J 


(3-44) 


In  principle,  it  is  straightforward  to  employ  the  signal  optimiza¬ 
tion  already  discussed  in  Section  3.4.  It  is  readily  seen  that,  again 
assuming  thermal  noise  from  the  terrain  dominates  that  from  the 

-  2 

receiver,  the  optimum  pulse  modulation  f  is  specified  insofar  as  |  fQ|  , 
which  is  given  by 


f 

o 


2 


=  max 
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(3-45) 
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Additionally,  u  is  restricted  to  be  in  the  set  u  e  ,  where 

(arM  |2  t  0. 

Implied  in  the  above  discussion  is  that  the  optimum  filter  K  will 
also  contain  this  "comb  structure"  Q(co),  see  Equation  3-34  with  f  = 
Q  f  .  This  "comb  structure"  Q  fluctuates  very  rapidly,  of  the  order  of 
.jjj/N,  and  hence  so  would  the  optimum  |fQ(w)l  .  For  example,  if  T  =  1 
ms,  then  f$  ■  1/Ts  *  1  KHz  and  this  is  quite  small  relative  to  the 
allowed  bandwidth  interval  ft*. 

Conversely  |f  (w)|  '  will  fluctuate  very  rapidly  general  ly  with 
respect  to  its  permitted  bandwidth  ft*.  This  results  in  a  tendency  for 
fQ(t)  to  have  a  long  time  extent. 


For  example,  if  the  required  energy  E  is  relatively  quite  small, 
and  £PC,  |B|  ,  and  | a t !  are  uniform  over  ft*,  then  !fQ(w)|2  tends  to 

have  its  support  on  narrow  intervals,  say  of  width  A,  centered  on  k^, 
where  k  is  such  that  k..  e  ft*  and  also  a«j  /N.  That  is,  roughly, 

“  a 


U) 


2 


£ 


(3-46) 


Suppose  ft*  =  { Lo :  -ft/ 2  +  <^<0^  +  ft/2  .}  Then,  if  we  assume  f  is  real 

(this  results  in  minimal  time  extent  for  f),  we  find 


f  !t)  *  IHLilit/A) 

0  irt 


sin  [(K  +  1 )  u,st/2] 
sTn  T*s’t/2) 


(3-47) 
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If  &  «  gj  /N,  then  2 tt/a  »  NT  and  there  Is  an  Inconsistency  with  the 
desire  to  have  a  periodic,  pulse-like  modulation. 

In  order  to  find  a  desirable  fQ(t),  it  Is  necessary  that,  e.g.,  Its 
time-extent  be  constrained.  In  order  that  It  represent  a  possible  pulse 
to  be  periodically  repeated.  This  problem  Is  discussed  In  the  next 
section. 

(It  may  be  seen  that,  as  -*■  0,  f  tends  toward  |Q|~*  (in  its 
dependence  of  Q).  That  Is,  £  t  and  fQpt  bear  an  inverse  relationship 
with  respect  to  |Q|.  The  "overall"  dependence  of  the  system  on  |Q| 
tends  to  disappear.  This  is  actually  reasonable.  It  is  only  the  total 
energy  of  f  that  is  constrained  and  SNR  becomes  independent  of  f  . ) 


3.5.1  FORMATION  WITH  ENERGY  AND  TIME-EXTENT  CONSTRAINTS 

When  a  periodic  pulse  modulation  Is  imposed  for  reasons  such  as 
discussed  above,  we  also  often  want  to  constrain  the  time  extent  of  the 
periodically  repeated  pulse,  typically  to  a  small  fraction  of  the 
period.  Thus,  analytically,  we  would  straightforwardly  impose  another 
constraint.  We  will  discuss  two  ways  this  can  be  done  in  a  reasonably 
tractable  manner,  though  only  one  is  considered  at  length  here. 

The  most  straightforward  method,  perhaps,  is  to  use  a  "radius  of 
gyration"  measure,  t,  of  the  time  extent  of  the  pulse  modulation, 
defined  by 


A  -  j  dt  (fj 

—  00 

assuming  that 


(3-48) 
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=  0 


l.e.,  the  centroid  of  |f (t)  j2  is  at  zero.  The  denominator  is  the 
energy,  Ef,  of  f  which  is  also  to  be  constrained.  Also  by  Parseval's 
relation, 


which  Is  therefore  a  reasonable  form  to  constrain. 
Forming  the  Lagranglan 


*<?>  •& 


r  IBJfl2  du, 

J  p„  *pclfl2 


f'U\2  du  (3-50) 


where  we  have  used  IbJ2  -  I B ! 2  latJ^»  pn  =  ^  »  and  Pc  =  ,9£|at  |2- 

Furthermore,  x^  and  x^  are  Lagrange  multipliers*!!^  f  is  to  be  chosen  to 
maximize  while  satisfying  the  two  constraints  on  energy  and  time 

extent. 


In  a  classical  way,  we  are  led  to  the  Euler-Lagrange  necessary 
condition  for  f  to  be  a  stationary  point  of  j2?(f).  If  fp  is  the  real 
part  of  f,  then  fp  must  satisfy 


X2  fR  +  X1  fR  ' 


!3, |P„  f 

_ 1 1  n  r 

pn  +  Pc*?!' 


=  0 


(3-51) 
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The  same  relation  holds  for  fj.  This  is  a  nonlinear,  ordinary 
differential  equation  for  which  suitable  boundary  conditions  are  to  be 
Imposed,  e.g.,  f  having  support  confined  to  Q  may  be  desirable.  The 
solution  of  this  problem,  which  can  only  be  done  numerically,  guided  by 
fairly  sophisticated  analysis.  Is  beyond  the  scope  of  this  study. 

To  Impose  a  constraint  on  the  modulation  time-extent  that  Is  more 
tractable  while  still  reasonably  meeting  the  constraint  requirement,  we 
note  that,  If  f  has  support  jt|  <T,  then  4>  ,  where  <T  =  |f|2,  has  support 
confined  to  \  t  \  <  2T.  So  we  impose  the  constraint  directly  on  <j>  .  To 
avoid  a  nonlinear  differential  equation,  we  use  the  "equivalent 
rectangle"  definition  of  time  extent: 


t2 


/l*(t)|2  dt 
Mo)  |2 


(3-52) 


The  denominator  is  again  recognized  as  the  energy,  E^,  of  f,  which  will 
also  be  constrained.  We  are  therfore  led  to  the  two  constraints: 


_1_ 

2tt 


2tt 


2 


duj  =  tE 


2 

f 


and  also,  of  course, 


<t(ui)  >0,  U)  6  Q. 


(3- 53a) 


( 3- 53b ) 


(3- 53c) 
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We  are  now  formally  led  to  the  problem:  maximize  by  the  choice  of 

where 


I B  I  ^  <f  dca 


Pn  +  P„  4> 
n  c 


(3-54) 


We  must  first  be  concerned  by  the  possible  conflict  of  ( 3- 53a )  and 
( 3- 53b ) .  $  which  satisfy  3-53a  and  3- 5 3b  will  be  called  "admi ssable" . 
Intuitively,  from  the  usual  theory  of  the  Fourier  transform,  we  know 
that  any  function  and  its  transform  cannot  be  simultaneously  "too 
constrained"  in  time  extent  and  bandwidth.  Here,  the  consequences  are 
readily  seen.  Dividing  (3-53b)  into  the  square  of  (3-53a),  we  have 


(3-55) 


by  the  Schwarz  inequality. 


p 


J{  1)^  du)  •  _/" $ ( u>)2  dw 

ft _ ft _  _  yi(ft) 


(3-56) 


where  y(  •)  is  the  measure  (interval  length)  of  the  set  ft  and  equality 
holds  if  and  only  if  <f(oj)  =  constant  for  w  e  ft.  If  the  support  of  ?,(w) 
is  confined  to  ft*  c  ft,  the  same  bound  occurs,  with  ft*  replacing  ft. 
Thus,  if  the  "bandwidth"  u(ft*),  the  (Lebesque)  measure  of  the  set  ft*, 
ft*c  ft,  is  too  small,  i.e.,  if 
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u(ft*)  <  2u/t,  (3-57) 

then  there  are  no  admissible  functions  T.  Further,  often  we  are 
Interested  In  meeting  another  constraint  not  explicitly  mentioned: 

corresponding  to  a  peak  power  restraint.  This  constraint  Is 
conventionally  met  by  use  of  "large-time  bandwidth  product"  signals, 
l.e.,  here  we  may  well  want 

x  »1,  ft*c  si.  (3-58) 

Note  also  that  If  5* Is  a  constant  over  a  set  ft*  c  ft  ,  then  necessarily 

p  =  =  i/x  (3-59) 

ci\ 

therefore,  such  T do  not  seem  to  be  all  that  interesting,  as  again,  they 
will  be  of  "unity  time-bandwidth  product." 

3.5.1  SOLUTION  FOR  OPTIMUM  SIGNAL 

Considering  the  maximization  of  (3-54),  we  see  that  the  Lagrangian 
is  the  form 


2 


=  h  f  F(?’  ^  ^ 


(3-60) 


where  suppressing  its  aj-dependence,  F  functionally  depends  on  $  as 


c,x 


F(x)  “ITS  -  C3X  -  c< 


(3-61) 


where  Cj  -  |Bj|  /Pn,  c2  =  pc/pn»  c3  =  xi/pn  and  c4  =  X2^Pn’  we  assume  pn 


80 


Jm 


is  Independent  of  u.  It  Is  clear  that  9S  is  maximized  over?'  by 
choosing  each  ?  value  to  maximize  F(?(w),  w  )  for  each  w  e  n*  c  n  . 
That  is,  we  wish  to  find  x  >  0  maximizing  F(x)  for  a  set  (cj,  c2,  c3, 


Me  compute 

C.  ~  (1  +  CnX)  (Co  +  2C-.X) 

F  1  (x)  \ - 2- 

(1  +  C2x)^ 

the  equation  for  stationary  points, 

F  * ( x)  =  0, 


(3-62) 


has  the  same  solution  (x  f  -1/c)  as  the  equation 

(c3  +  2c4x)(l  +  c2x)2  -  Cj  *  0.  (3-63) 

This  Is  a  cubic  equation  whose  roots  may  be  found  by  a  well  known 
formula . 


We  can  demonstrate  that  the  solution  to  Equation  (3-63)  may  have 
one  positive  solution  and  two  negative  solutions  (possibly  the  same). 
We  also  see  that  it  has  one  (at  most)  positive  solution,  and  then  only 
if  c3  <  Cj.  Therefore,  if  c3  >  Cj,  we  take  xQpt  =  0.  Otherwise,  we 
solve  the  cubic  equation  for  its  single  positive  root.  We  notice  also 
that,  keeping  in  mind  Cj  represent  j  B ^ |  2/ Pn »  which  depends  on 

' r  >  n'  -  for  which  a  nonzero  x 


smaller  we  make  c3(c3  >  0),  the  more 


opt 


the 
i  s 


found.  Further,  making  c^  very  large  tends  to  make  xQpt  very  small. 
These  two  observations  elucidate  how  a  broad  bandwith,  hence  small 


time-extent,  solution  can  occur. 
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3.5.2  SPECIFIC  NUMERICAL  RESULTS 

It  Is  now  straightforward  to  numerically  find  =  |f  |  ^  for  a  given 
set  (Cj,  c^).  We  will  employ  numerical  values  already  given  in 
Section  3.  4. 


The  energy  constraint  is 


( v ) 


the  time-extent-related  constraint  is 

«*{v) 

and  the  SNR  is 


(3-64) 


(3-65) 


2  106  f  I  c.(v)?(v)  ) 

d  =_2 V  J  |l"+  c-(v)f(v)  '  c3»(v)  •  V(v)  dv- 
n*(v) 


A  series  of  optimal  signal  spectra  were  computed  for  numerous 
energy/signal  duration  extents,  E^/t.  Selected  examples  are  presented 
in  Figures  3-13  and  3-14.  It  was  determined  that,  generally  speaking, 
the  time  extents  were  of  a  "practical"  nature,  ranging  from  about  5  x 
10  ^  sec  to  314  x  10"®  sec.  The  narrower  signal  extents  are  associated 
with  rather  large  energies,  providing  SNRs  much  larger  than  would  be 
required  for  adequate  detection  probability.  Moreover,  as  the  signal 
energy  constraint  decreases,  the  optimum  signal  spectra  "gravitates" 
from  the  lower  end  of  n*,  (near  the  maximum  (resonant)  response  of  the 
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b)  Energy  Constraint  -  0.4  Joules,  Pulse  Extent-  1 .57x1 0"7  sec,  d2-43.9dB 

If  Ml 2 


Figure  3-14.  Effects  of  Lower  Energy  Constraints  on  Periodic 
Signal  Design  (^min  -25,  ^^-200) 
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reflecting  object),  to  an  Increasingly  narrow  support  near  the  higher 
frequency  (resonant)  response  of  the  object. 

The  results  of  this  computation  illustrate  that  there  is  a  strong 
spectral  dependence  of  the  optimal  signal  in  the  assumed  energy/pulse 
duration  constraint.  Undoubtedly  other  target  cross  section  and  antenna 
response  functions  would  have  yielded  different  signal  spectra,  although 
the  trend  would  undoubtedly  be  to  follow  the  resonance  of  the  target 
response,  if  it  were  within  the  overall  system  passband.  (It  is 
surprising,  however,  that  the  spectra  gravitates  towards  the  higher 
resonance,  despite  of  the  v5  clutter  law  dependence.)  Further  results 
should  be  computed  for  specific  targets  that  are  of  interest  for 
military  applications. 

3.6  OPTIMUM  SIGNALS  FOR  DISPERSIVE  PROPAGATION  MEDIA 

3.6.1  INTRODUCTION 

In  some  situations,  the  propagation  path  may  contain,  perhaps  in 
part  of  its  extent,  a  material  consistency  that  is  described  by  a 
complex  dielectric  constant  that  depends  significantly  upon  wavenumber, 
or  frequency.  Such  constituents  may  be  due  to  suspended  molecules  (such 
as  oxygen  «nd  water  vapor  normally,  present  in  the  atmosphere),  suspended 
particulates  (water  droplets  and  aerosols,  e.g.,  in  the  atmosphere),  and 
also  peculiar  propagation  media  such  as  plasmas  (e.g.,  the  ionosphere) 
and  soils  (e.g.,  when  objects  of  Interest  may  be  buried). 

Because  of  the  causality  property  of  Maxwell's  equation  for  the 
electromagnetic  field,  the  real  and  imaginary  parts  of  a  complex 
dielectric  constant  are  necessarily  related.  That  is,  dispersion  and 
absorption  (attenuation)  are  encountered  together.  Of  course,  the 
dispersion  effect  is  to  degrade  the  pulse  shape  and  hence  detection  and 
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resolution  capabilities  of  a  radar.  The  effects  of  attenuation  are  to, 
in  effect,  negate  radiated  power  and  hence  degrade  the  detection 
performance.  When  present  with  sufficient  influence,  these  effects  can 
dramatically  effect  the  optimal  signal  design. 

To  exemplify  the  impact  on  the  signal  design,  we  will  consider 
dispersion  due  to  a  molecular  constituency  present  in  the  propagation 
path.  Further,  we  shall  consider  the  effect  of  a  single  "resonance"  or 
"absorption  line."  (The  extension  of  these  results  to  multiple 
resonances  is,  in  principle,  straightforward. )  For  simplicity,  we  shall 
assume  either  that  the  entire  propagation  path  is  homogeneous  or  that  a 
"slab"  of  dispersive  media  is  present  in  such  a  way  that  diffraction  (at 
the  boundaries)  need  not  be  considered.  Also,  it  would  often  be  the 
case  that  the  density  (in  gm/cm  )  of  the  molecular  constituents  would  be 
low,  so  that  the  resultant  wavenumbers  (for  plane  waves)  do  not  differ 
greatly  from  those  for  an  in- vacuo  path. 

A  review  of  the  previous  development  shows  that,  generally  speak¬ 
ing,  all  results  carry  forward  with  the  wavenumber  k  replaced  by  k(w), 
as  determined  by  the  dispersion  relation.  The  end  effect  is  that  the 
"block  diagram"  for  the  overall  system  model  will  have  an  additional 
"box",  a  linear,  coordinate  invariant  transformation  characterizing  the 
dispersion  and  attenuation  effects.  Further,  then,  as  we  have  shown 
that  only  the  modulus-squared  of  the  transfer  function  of  these  "boxes" 
enters  into  the  signal-to-noise  ratio  formula,  given  that  the  processor 
is  the  optimal  matched-filter,  as  discussed  previously.  Thus,  it  will 
be  only  the  attenuation  that  will  enter  into  the  optimum  signal  design 
problem.  We  re-emphasize  that,  among  the  a  priori  data  that  is  assumed 
known,  will  now  be  k(j),  the  dispersion  and  attenuation  relations. 
This,  in  a  particular  situation  may,  or  may  not  be,  a  reasonable 
assumption. 
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It  might  be  that  only  an  a  priori  knowledge  of  certain  statistical 
distributions  about  k(w)  are  reasonably  assumed.  In  such  an  instance, 
the  "signal"  and  "terrain-emitted  noise"  will  be  correlated  as  both 
traverse  the  same  dispersive  medium.  Then,  the  likely  difficult  optimal 
approach  must  follow  the  path  generally  discussed  above,  or  more  simply 
and  sub-optimally,  we  may  simply  use  an  averaged  k(w)  in  place  of  the 
unknown  random  k(w). 


We  begin  by  briefly  reviewing  a  standard  molecular  absorption  model 
that  will  serve  us  "canonically"  here  [24,25]. 

3.6.2  A  SIMPLE  DISPERSION/ABSORPTION  MODEL 


One  considers  the  propagation  of  plane  waves  in  an  isotropic  medium 
containing  a  dielectric  constant  and  it  is  found  that  the  phase  velocity 
of  a  plane-wave  can  depend  upon  an  index  of  refraction  dependence  upon 
frequency.  The  fundamental  model  for  this  phenomenon  depends  upon  a 
dynamical  model  for  the  molecular  constituents  in  the  medium.  The 
simplest  such  dynamical  model  is  a  simply  resonant  oscillator  possessing 
a  free  period  which  is  excited  by  the  t ra ve r s i n g/ i nc i d e n t 
electromagnetic  field.  The  resulting  complex  wavenumber  that  is 
associated  with  a  permissible  solution  of  the  wave  equation  may  have  a 
form  such  as 


-  k(oi)  = 

0) 


(3-66) 


p 

where  ^  is  the  resonant,  frequency,  bj  is  proportional  to  the  number  of 
such  resonant  molecules/oscillators  per  unit  volume,  and  b^  accounts 
for  the  dissipative  effect  of  molecular  collisions.  (The  in-vacuo 
wavenumbers  permitted  are  k  =  ^/c  as  noted  and  used  above.)  The 
determination  of  these  constants  in  any  specific  instance,  is  an  object 
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of  theory  and  experiment.  Resolving  such  a  k(w)  into  its  real  and 
Imaginary  parts  results  in,  respectively,  the  "dispersive  part" 


U)  “  GO 


-  Re{k(cj) }  =  n  =  1  + 

'jj 


UJ_ _ 

2,2  2x2  ,  2.2 

-  a)  )  +  a)  D 2 


_  C 

=  —  a 

to 


(3-67) 


and  the  "attenuation  part" 


ykun 


bl  b2 

2c  ,  2  2x2  .  2.  2 

(u)  -  u)  J  +  uj  D0 


=  6 


(3-68) 


3.6.3  EXTENSION  OF  RADAR  SYSTEM  MODEL 

If  we  review  the  derivation  of  the  basic  system  model  employed 
here,  under  the  assumption  that  k(w)  does  not  depart  grealy  from  k  = 
w/c,  we  see  that  we  need  only  replace  k  by  k(w)  in  all  forms.  Thus,  the 
form  (Equation  3-8)  for  the  received  signal  waveform  is  now 

i  f  ^  -lut+21  R  k(io)  _ 

s(t)  =  /  doo  f (cu)  e  0  Q(oj)u(0 ;u)  (3-69) 

where,  as  stated,  in  the  representation  of  u(o;  w)  in  Equation  3-1,  k  is 
everywhere  replaced  by  k(w). 

Generally,  if  the  dispersive  effects  are  "weak",  we  can  write 

k(uO  =  1  +  kD(u.)  (3-70) 

In  the  specific  instance  of  molecular  absoprtion  mentioned  above  this  is 
the  case,  with  the  obvious  definition  of  k0(u>).  Thus,  we  may  write 
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a(io)u(0;  oj) 


2pL 


s(t  +  V  "  7J  7d“ 

(2 

provided  tq  =  2RQ/c.  The  "new11 
Figure  3-15. 


i(i)t 


i  2RQkD  (</j) 


system  block  diagram  is  as  shown  in 


It  is  clear  that,  under  the  assumptions  stated  above,  that  any 
region  traversed  by  the  ray  path  does  not  cause  diffraction  effects,  if 
a  subpart  of  the  ray  path  is  the  only  part  in  the  dispersive  region,  and 
is  of  "thickness"  r  .  Then  we  can  replace  ( 2RQ)  by  (2r  )  in  the  above. 

We  see  that  all  of  the  earlier  development  now  carries  over  if  we 
simply  replace  at(^)  everywhere  by  at(a))  hp(w),  where  h  D  U )  = 
exp[i2r0kD(u)]* 

We  thereby  can  readily  include  the  dispersion  and  attenuation  due 
to  a  frequency-dependent  dielectric  constant  characterizing  (part  of) 
the  propagation  path. 

The  optimal  receiver  contains  a  filter  matched  to  the  received 
waveform,  s(t).  Assuming  that  hg(^),  that  is,  kg(ui),  is  known,  then  the 
overall  system  transfer  function  and  the  SNR  maximized  over  receiver 
choice  depends  functionally  only  upon 

I hD(oj)  | 2  =  exp[-4rok0I(u)]  (3-71) 
where  kgj(w)  =  Im  (kpj(oj)}  *  8(u>) 


*To  correctly  model  the  thermal  noise,  that  component  due  to 
terrai r,-emi tted  thermal  radiation  should  also  be  modified  by  its 
propagation  through  the  dispersive  medium.  Thus  we  should  replace 
by  &  |  expfiR  k  (<*>)]  I2  where  R  is  the  length  of  the  one-way  nter 

ter  0  u  0 

propagation  path. 
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Figure  3-15.  Radar  System  Model  For  Weak  Dispersive  Media 
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Of  course,  it  also  remains  true  that  the  SNR,  maximized  over 
receiver  choice,  depends  only  upon  |f(w)|  ,  the  magnitude-square  of  the 
signal  spectrum.  Thus,  we  may  utilize  the  same  algorithms  for  this 
optimization  problem  solution  by  just  replacing  lazuli  by 

|StMI2  lhD(»)l2. 

In  particular,  we  shall  now  carry  over  the  choice  of  an  optimum 
signal  under  energy  and  time  extent  restraints. 

3.6.4  SPECIFIC  NUMERICAL  CALCULATIONS 


To  exemplify  these  matters,  we  shall  employ  a  single,  simple 
resonance  with  parameters,  attenuation  and  bandwidth,  chosen  in  a 
reasonable  way,  namely  to  be  comparable  to  those  of  the  oxygen 
absorption  "line"  for  a  "normal"  atmosphere.  However,  we  shall  choose  a 
resonant  frequency  considerably  lower,  in  the  general  band  of 
frequencies  we  have  considered  in  the  principal  example  studied  in 
Section  3-5.  We  model  as 


8(w) 


8 


0 


+ 


(3-72) 


with  s(w  )  =  max  [b(uj)]  =  10  ,  corresponding  to  a  10  dB/Km  attenuation. 

0  uj  ^ 

The  width  Af  =  Au)/2ir  of  the  oxygen  line  is  about  b  x  10  Hz.  We  shall 
take  Af  here  a  few  MHz  for  a  "cleaner"  effect. 


In  view  of  the  specific  object  (scatterer)  "resonance"  curve  chosen 
in  Section  3-4,  we  here  will  find  it  especially  interesting  to  choose  f 
=  w  /2  ff  equal  to  one  of  those  resonance  lines,  namely,  150  MHz.  In 
terms  of  the  normalized  frequency  v  =  f/ 10  •  Equation  3-72  becomes 

j3 ( v )  =  -7 - — - To -  1  v,  AV  in  MHz  (3-73) 

[(v  -  v0)/avJ^  +  1 
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For  the  10  Km  path  length  RQ  used  in  above  examples,  we  find 


exp[-4R03(v)]  =  exp< 


(3-74) 


This  is  a  rather  severe  attenuation  and  we  may  find  it  convenient  to 
replace  (-40)  by,  e.g.,  -4,  corresponding  to  a  less  dense  dispersive 
constituency  or,  e.g.,  a  1  Km  thick  medium. 


Typical  results  of  the  signal  optimization  are  shown  in  Figure 
3-16.  As  expected,  the  signal  modulation  is  chosen  as  in  the 
non-dispersi ve  case,  except  that  no  signal  spectrum  is  "wasted"  where 
the  attenuation  is  too  great.  Also,  there  seems  to  be  a  tendency  toward 
a  curious  increased  weighting  near  the  edge  of  the  heavily  attenuated 
region. 


3.7  FURTHER  DESIGN  OBSERVATIONS 


The  preceding  developments,  in  particular  the  design  examples,  show 
that  the  optimum  signal's  spectrum  modulus  can  assume  a  rich  variety  of 
forms,  depending  upon  the  signal  energy  and  pulse  duration  constraint, 
the  maximum  permitted  bandwidth,  the  nature  of  the  path  dispersion,  and 
the  remaining  free  system  parameters.  Additionally,  the  spectral 
modulus  is  strongly  influenced,  as  was  seen,  by  the  form  of  the 
frequency  response  of  the  object  to  be  detected,  denoted  |w(v)|2. 

This  being  the  case,  the  question  naturally  arising  next  concerns 
how  one  might  design  such  a  detection  system  for  a  realistic  practical 
application.  While  the  designer  desires  an  optimal,  or  near  optimal, 
design,  the  latter  should  also  possess  a  certain  robustness  with  respect 
to  some  degree  of  departure  of  parameters  from  their  a  priori  assumed 


92 


50 


100 


*-►  >p(hhz) 


a)  Energy  Constraint 


136  Joules,  Pulse  Extent— 


1.8x10 


sec 


b)  Energy  Constraint 


106  Joules,  Pulse  Extent - 


2.6x10 


sec 


If  HI 2 


c)  Energy  Constraint  —  30.2  Joules,  Pulse  Extent-  3.0x10  7  sec 
Note:  Maximum  attenuation  of  media  is  at  150  MHz. 


Figure  3-16.  Effect  of  Energy  Contraint  on  Optimum  Signal 
For  Dispersive  Media  ('Pmin*25' 
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values.  The  most  straightforward  procedure  is  to  assume  an  a  priori 
distribution  is  known,  or  can  be  estimated,  for  such  parameters.  Then, 
"standard"  statistical  procedures  can  be  employed.  We  will  sketch  this 
out  in  the  case  of  variable  object  (target)  frequency  responses. 
Clearly,  there,  if  we  design  a  system  for  certain  resonances  of  one 
object,  the  detection  performance  may  be  very  poor  relative  to  other 
objects  with  differing  resonances.* 


For  example,  It  is  conceivable  that  a  particular  radar  system  might 
be  dedicated  to  detecting  an  object  which  is  a  member  of  a  small  number 
of  possible  classes.  It  may  also  be  reasonable  to  assume  estimates  of 
the  a  prior  probability  of  a  given  class  can  be  made  available.  Roughly 
speaking,  we  might  expect  that  an  appropriate  system  design  would  then 
be  concerned  with  an  "average  object  response",  in  place  of  a  single 
object  class  response.  This  will  be  the  case  under  the  "soft"  maximum 
signal-to-noi se  ratio  (suitably  defined)  criterion  as  we  shall  see. 
However,  this  is  likely  not  an  optimal  procedure  when  the  optimality 
criterion  is,  e.g.,  maximizing  the  decision  error  probability. 


We  now  revisit  the  basic  hypothesis  testing  problem.  We  fix  a,  the 
class,  and  actually  repeat  the  considerations  implicit  in  the  prior 
discussion,  finally  averaging  over  a.  That  is,  the  likelihood  ratio, 
for  fixed  a  and  random  phase  e,  is  well  known  [22,23]  to  be 


,  Fe19  r  ,  * 

#.  ( z  j  8  ,  a)  =  exp<(Re  {  —  -  J  z(t)S  (t)  dt 

Ln°  T 


J 

2n" 


Jiyt) 


dt 


(3-75) 


where  z  denotes  the  data 


z(t)  s  e1'0  S  (t)  +  rift),  t  e  T, 

a 


(3-76) 


*Happi ly,  however,  resonances  are  strongly  related  to  an  object's 
size,  so  that  certain  classes  (sizes)  of  military  objects  (e.g.,  manned 
bombers)  are  likely  to  have  similar  responses. 
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where  n(t)  is  a  sample  function  of  white  Gaussian  noise,  and  a  signifies 
class  membership.  Carrying  out  the  average  over  6,  one  obtains  the  well 
known  form  [22,23] 

t(z|o)  =  'ofe  <<tj)  e'E“A'°  (3-77) 

where  IQ  is  the  modified  Bessel  function  of  the  first  kind,  order  zero, 
and  the  real  signal  s  (t),  corresponding  to  the  complex  signal  S„(t), 
has  energy 


E 

a 


/iyy2 

T 


dt  . 


(3-73) 


Suppose  that  the  number  of  classes  is  countable,  with  a  priori 
probabilities  P^,  then  the  desired  likelihood  ratio,  which  specifies  the 
optimal  processor,  is 


l(z) 


(3-7S) 


Thus,  this  optimum  processor,  for  the  detection  purpose,  is  in  fact  the 
approximately  weighted  sum  of  all  the  subprocessors  required  in  order  to 
perform  object  classification.  Generally,  because  of  the  nonlinear 
forms,  a  simpler  form  is  not  assumed.  In  principle,  the  error 
performance  can  now  be  determined.  The  signal  is  then  chosen  to 
minimize  the  error. 
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For  example,  suppose  that  the  probability  of  error,  given  a,  has  a 
bound  and  approximation  of  the  form 

P(Eja)  <  1/2  exp[-d2/4]  (3-80) 

as  is  precisely  the  case  when  "signal  plus  noise"  and  "noise  alone"  are 
equally  likely  [22,23],  Averaging  over  the  a  prior  probability  of  the 
classes,  the  desired  error  bound  is  of  the  form 

P(L')  <  Pa  exP  C-<£/4].  (3-81) 

a 

Here  d^  is  given  by  (3-34),  with  [  B |  replaced  by  |B  |  .  The  optimal 
signal  spectrum  modulus  is  now  reasonably  chosen  to  minimize  the 
bound/approximation,  e.g.,  with  an  energy  restraint.  This  appears  to  be 
a  more  difficult  optimization  problem  than  that  already  considered  above 
and  we  do  not  pursue  it  further  here. 


A  simple  approach  to  which  it  may  be  necessary  to  resort  anyway 
when  the  distribution  laws  of  the  "noises"  are  non-Gaussi an  or  even 
unknown,  is  to  choose  the  signal  in  order  to  maximize  a  signal -to-  noise 
ratio  (SNR).  We  carry  over  the  definition  already  made  above,  except 
additionally  we  average  over  the  a  prior  distribution  of  a.  As  we  have 
seen,  for  fixed  a,  we  have  the  SNR 


j2  _  J_ 

a  2tt 
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B  |2 
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It  is  therefore  clear  that  the  SNR 

dn  =  Er 

a  a 
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2  2 
has  the  same  form  as  d  employed  earlier,  except  that  |B  |  Is  replaced 

by  the  average 


Therefore,  the  signal  spectrum  modulus,  optimum  according  to  this  SNR 
criterion.  Is  first  the  solution  we  found  before,  except  |B|  is  now 
replaced  by  Ea{lBa|2|>  the  expected,  or  averaged  object  frequency 
response  (modulus-squared).  Of  course,  this  is  an  intuitively  agreeable 
result. 

3.8  CONCLUSIONS 

The  foregoing  development  is  typical  of  the  analysis  that  can  and 
should  be  made  for  any  radar  system.  The  basic  procedure  is  to,  first, 
establish  as  comprehensive  an  analytic  model  as  possible  for  the  radar 
system  in  its  entirety  and  then,  second,  establish  a  design  criterion  by 
which  the  processor  and  waveform  modulations,  and  perhaps  other  system 
parameters,  are  chosen.  Naturally,  the  success  of  this  approach  will 
depend  upon  the  adequacy  of  the  model,  the  appropriateness  of  the 
criterion,  and  the  ability  to  perform  the  requisite  analysis  including 
optimization. 

Here  a  "monostatic"  radar  detection  system  was  modeled  rather 
completely  with  respect  to  a  matter  of  principal  concern:  the  system's 
overall  wavenumber,  or  frequency,  dependence.  The  employed  receiver  and 
decision  structure  is  optimal,  according  to  the  Neyman- Pea rscn 
criterion,  under  normal  distributions  and,  additionally,  in  the  absence 
of  normal  distributions,  is  optimal  under  the  softer  maximum  signal-to- 
noise  ratio  criterion.  The  optimality  criterion  dictated  how  the 
signal,  or  modulation,  should  be  designed.  Under  a  constraint  on 
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allowable  signal  energy,  this  optimization  problem  has  a  known  solution 
which  can  be  found  numerically  by  a  "bisection  type"  algorithm.  It 
turned  out  that  solely  the  modulation's  spectral  magnitude  is  thereby 
specified,  leaving  other  attributes  free  to  perhaps  meet  other  desirable 
properties  of  the  modulation,  e.g.,  Doppler  resolution.  An  approximate 
"synthesis"  procedure  was  given  to  realize  such  waveforms  in  the  large 
time-bandwidth  product  case  (Appendix  B). 

The  foremost  Importance  of  the  work  given  here  in  the  present 
context  is,  then,  the  optimal  signal/modulation  spectrum  modulus.  Its 
"support",  i.e.,  frequencies  where  its  non-zero  values  reside,  is 
thereby  given  Indirectly,  without  any  a  priori  explicit  restrictions. 
The  radar  system  parameters  entered  the  signal  optimization  problem  as 
two  parameters  plus  the  energy  constraint  parameter.  As  these  were 
varied,  it  was  observed  that  the  signal  spectral  modulus  exhibited  a 
wide  range  of  interesting  behavior,  occupying  location  at  or  near 
"resonances"  of  the  object  to  be  detected  and  which  lay  within  the 
system  passband,  and  varying  from  narrow  to  broad  with  even  disjoint 
supports.  Not  all  of  these  variations  could  be  predictable  a  priori, 
apparently,  demonstrating  the  value  of  this  development. 

Additional  results  were  obtained  for  periodic  (pulse)  modulation 
formats  as  well  as  including  the  effects  of  weak  dispersion  in  the 
propagation  path.  It  is  concluded  that  solutions  to  the  rich  signal 
design  problem  for  detection  of  targets  In  clutter  does  Indeed  yield,  in 
some  cases,  signals  that  are  of  wide  relative  bandwidths.  Thus,  future 
applications  analyses  for  specific  targets  of  interest  should  be 
conducted  to  quantify  the  expected  utility  of  such  an  approach. 

The  next  section  critically  reviews  the  available  literature  on  NSR 
techniques  and  on  related  topics  which  strongly  depend  on  relative 
wideband  waveforms. 


4.0 

DISCUSSION  OF  SELECTED  WIDEBAND  RADAR  (WBR)  RESEARCH 


Section  2  outlined  the  principal  radar  applications  for  which  WBR 
systems  offer  potential  for  improved  performance  and  summarized  the 
capabilities  and  limitations  of  current  narrowband  radars  in  these 
areas.  In  Section  3,  a  general  theory  for  detection  of  targets  in 
clutter  was  presented  for  systems  of  arbitrary  bandwidth.  This  theory 
Is  relevant  to  several  of  the  radar  applications  identified  in  Section 
2,  and  serves  as  a  means  for  evaluating  present  research  efforts  into 
WBR  methods. 

This  section  reviews  selected  areas  of  current  WBR  research  and 
compares  the  potential  of  the  techniques  with  the  capabilities  of 
present  narrowband  radar  systems.  Although  the  review  is  by  no  means 
exhaustive,  the  areas  considered  herein  are  representative  of  research 
which  most  significantly  deviate  from  the  standard,  narrowband  systems 
approaches.  Specifically,  they  were  chosen  for  having  satisfied  one  or 
more  of  the  following  criteria:  (1)  a  novel  approach  which  cannot  be 
realized  with  conventional,  narrowband  systems,  (2)  a  controversial 
method  for  which  a  unified,  end-to-end  model  is  required  for  evaluation, 
(3)  a  technique  which  was  specifically  identified  for  consideration  as 
part  of  the  objectives  of  the  NSR  Techniques  research  program. 

To  this  end,  a  pair  of  topical  areas  in  WBR  research  are  discussed 
in  the  sections  that  follow.  The  first  encompasses  work  in  techniques 
based  on  the  singularity  expansion  method  (SEM),  which  relates  the 
transient  response  (radar  return)  of  a  target  to  a  set  of  poles  in  the 
complex  frequency  domain  [27-31].  These  poles  depend  solely  upon  the 
shape  and  material  parameters  of  the  target,  and  are  independent  of  the 
aspect  or  form  of  the  incident  radiation.  This  representation  of  a 
target's  radar  signature  suggests  several  potential  approaches  toward 
improved  performance  in  the  applications  outlined  in  Section  2.1;  these 
are  described  in  Section  4.1. 


A  second  research  area  has  been  expounded  almost  exclusively  by  H. 
Hannuth  and  Involves  the  use  of  what  have  been  called  "carrier-f ree"  or 
large  relative  bandwidth  signals  for  improving  performance  in  virtually 
all  radar  applications  [32-52],  This  work  has  been  the  subject  of 
extensive  controversy  [53,54],  and  a  need  to  resolve  these  issues  and 
identify  the  merit  of  Harmuth's  claims  have  been  explicitly  included  as 
part  of  the  objectives  of  the  program,  A  discussion  of  Harmuth's  work 
is  presented  in  Section  4.2. 

Section  4.3  compares  the  approaches  discussed  in  Sections  4.1  and 
4.2  with  the  capabilities  of  present  narrowband  radar  systems  on  the 
basis  of  technical  merit,  potential  performance  improvement,  and 
practical  implementation  limitations. 

4.1  SEM- BASED  METHODS 

The  singularity  expansion  method  was  first  applied  to  transient 
electromagnetic  problems  by  C.  Baum  and  his  colleagues  [27,28]  in  an 
effort  to  characterize  the  currents  induced  on  targets  by  a  nuclear 
electromagnetic  pulse  (NEMP).  More  recently,  the  SEM  has  been 
identified  as  a  useful  method  in  radar  applications,  particularly  for 
target  discrimination  and/or  classification  [55-67],  The  essence  of  the 
SEM  and  its  utility  in  target  ID  can  be  summarized  briefly  as  follows. 

Let  f(t)  represent  the  temporal  impulse  response  of  a  target,  i.e., 
f(t)  is  an  arbitrary  component  of  the  vector  field  scattered  by  a  target 
under  illumination  by  an  impulsive  plane  wave  incident  field  of 
arbitrary  polarization  and  direction  of  incidence.  For  t  >  t,  where  t 
is  a  sufficiently  large  delay  (late-time  response),  f(t)  can  be  expanded 
as  a  sum  of  complex  exponentials  [55,56],  namely 
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(4-1) 
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where  the  an  are  Independent  of  time.  The  natural  frequencies  sn 
correspond  to  poles  in  the  complex  frequency  response  F(s)  of  the 
target,  which  is  obtained  via  a  two-sided  La  Place  transform  on  f(t). 
The  sn  have  been  shown  to  be  independent  of  the  form,  direction  of 
Incidence,  or  polarization  of  the  incident  field.  Conversely,  they 
depend  quite  strongly  on  the  shape  and  material  composition  of  the 
scattering  body.  As  such,  they  are  characteristic  of  a  particular 
target  and  are  invariant  to  the  geometry  in  which  f(t)  is  measured.  All 
effects  of  the  geometry  and  incident  field  are  contained  in  the 
coefficients  an  [27-31]. 

This  invariant  property  of  the  natural  frequencies  is  the  principle 

upon  which  all  SEM-based  radar  applications  are  founded.  Through  either 

direct  or  indirect  measurement  of  a  subset  of  the  s  from  the  received 
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signal  of  an  impulse-like  radar,  discrimination  and  perhaps  even 
classification  of  the  target  which  produced  the  return  can  potentially 
be  achieved.  In  the  subsections  that  follow,  a  review  of  the  research 
programs  in  SEM-based  radar  at  specific  institutions  is  presented. 

4.1.1  RESONANCE  REGION  RADAR  (R3)  -  Research  at  General 

Research  Corporation  and  The  Naval  Postgraduate  School 

Perhaps  the  most  advanced  application  of  the  SEM  to  radar  is  the 
work  of  M.  L.  VanBlaricum  and  his  associates  at  General  Research 
Corporation  (GRC)  and  M.  A.  Morgan  and  his  co-workers  at  the  Naval 
Postgraduate  School  (NPS)  [29,55-57,65].  Beginning  with  VanBlaricum' s 
thesis  work  at  the  University  of  Illinois  [29],  a  systematic 
theoretical,  numerical,  and  experimental  research  program  has  evolved 
for  performing  target  identification  by  extraction  of  the  natural 
resonances  from  the  target's  time  domain  response  to  broadband 
excitation.  A  brief  outline  of  the  methodology  follows. 


The  transient  response  of  a  target  is  obtained  by  measuring  the 
scattered  field  produced  by  an  incident,  pulse-like  waveform  from  an 
i mpu 1 se.  generator  excitating  a  broadband  horn  antenna.  A  similar 
antenna  is  used  to  receive  the  scattered  field.  Alternatively,  a 
transfer  function,  which  is  the  Fourier  transform  of  the  transient 
response,  can  be  measured  by  replacing  the  Impulsive  incident  field  with 
a  series  of  CW  fields  produced  by  a  swept  frequency  transmitter.  In 
either  case,  calibration  data  from  a  known  scatterer  (sphere)  and  from 
the  environment  (chamber)  are  used  to  remove  the  effects  of  the  clutter, 
transmitter,  and  receiver  responses  from  the  measured  target  response. 
The  measurement  system  is  illustrated  in  Figure  4-1  and  the  calibration 
procedure  in  Figure  4-2.  The  calibration  is  performed  in  the  frequency 
domain. 

The  calibrated  data  represent  the  target  impulse  response  (or 
transfer  function)  to  an  excitation  with  a  uniform  spectrum  over  all 
frequencies  transmitted.  Typical  frequency  and  time  domain  responses 
for  a  generic,  simplified  target  model  are  shown  in  Figure  4-3.  This 
sampled  response  is  input  to  an  algorithm  which  estimates  the  natural 
resonance  frequencies  representing  the  poles  in  the  complex  spectrum.  A 
modified  version  of  Prony's  method  is  used  in  performing  the  required 
computations  [29,57]. 

The  target  di scrimi nation  function  as  demonstrated  in  the  most 
recent  GRC/NPS  results  [55-57]  is  achieved  using  a  two  step  process. 
First  theoretical  and/or  experimental  impulse  response  data  from  a 
variety  of  aspects  are  processed  using  Prony's  method  to  build  up  a 
histogram  of  pole  occurrences  in  the  complex  frequency  plane  for  each 
target  of  interest.  Based  on  the  pole  clustering  that  occurs  in  the 
histograms,  a  set  of  circular  regions  are  identified  about  each  pole 
cluster.  In  the  second  step,  the  poles  of  the  response  from  an  unknown 


Figure  4-1.  Experimental  facilities  for  time-domain  (TD)  and 

frequency-domain  (FD)  target  backscatter  measurements 
(after  [ 57 ]) - 


Figure  4-2,  Calibration  procedure  for  scattering  data  from  the  system  in 
Figure  4-1  (after  [57]). 


Figure  4-3.  Frequency  and  time  domain  responses  compared  to  theory  for  generic  model,  nose 
incidence  (after  [57]). 
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target  are  obtained  and  characterized  as  falling  either  inside  or 
outside  the  regions  for  each  of  the  targets  “cataloged"  in  the  first 
step.  A  positive  identification  is  assumed  when  the  number  of  poles 
falling  In  the  cluster  regions  of  a  particular  target  exceeds  the  number 
for  any  other  target.  Preliminary  results  from  the  laboratory  indicate 
reasonable  success  at  discrimination  between  two  classes  of  targets 
[65].  Failures  occurred  primarily  with  data  collected  at  nulls  in  the 
radiation  pattern  of  the  natural  resonances  (aspects  for  which  an  =  0, 
Eq.  (4-1)),  or  for  data  within  insufficient  signal-to-noi  se  ratio  (<  20 
dB). 


While  the  GRC/NPS  research  has  successfully  demonstrated  the 
applicability  of  the  SEM  to  aspect- i ndependent  target  identification  in 
the  laboratory,  there  are  several  important  considerations  which  must  be 
addressed  before  a  practical  vesion  of  a  resonance  region  radar  can 
become  operational.  These  are  discussed  in  Section  4.3. 

4.1.2  RADAR  WAVEFORM  SYNTHESIS  -  Research  at  Michigan 
State  University 

The  use  of  the  SEM  in  radar  applications  is  not  restricted  to  the 
explicit  calculation  of  natural  resonance  poles  from  radar  backscattered 
signals.  Indeed,  VanBlaricum  has  suggested  that  the  eventual 
operational  radar  system  for  target  identification  based  on  the  research 
described  in  Section  4.1  would  not  rely  directly  on  pole  extraction  from 
each  transmitted  pulse  [57].  Instead,  a  data  base  of  known  pole 
distributions  would  be  used  to  perform  a  "correlation-  prediction"  with 
the  return  from  an  unknown  target  to  arrive  at  a  "most  likely"  target 
identification.  In  this  approach,  the  incident  waveform  would  not  be 
tailored  to  a  specific  target  or  class  of  targets,  except  perhaps  to  the 
extent,  of  selecting  the  appropriate  spectral  content  of  the  transmitting 
signal  so  as  to  maximize  the  sigral-to-noi se  ratio  in  the  backscattered 
return  (see  Section  3). 


106 


An  alternative  approach  can  be  taken  when  It  is  necessary  to 
discriminate  a  single  target  or  set  of  targets  from  a  larger  class.  In 
this  case,  the  transmitted  and/or  received  signal  can  be  tailored  to 
result  In  a  return  which  Is  characteristic  of  the  target  of  interest  and 
at  the  same  time  significantly  different  from  the  returns  due  to  other 
scatterers.  This  concept  of  "Radar  Waveform  Synthesis"  has  been  the 
subject  of  significant  research  activity  at  Michigan  State  University 
under  K.  M.  Chen  and  his  colleagues  [58-60].  The  foundation  of  their 
approach  is  again  based  on  the  SEM.  In  this  case,  a  transmitted 
waveform  is  designed  such  that  a  single  natural  resonance  frequency  is 
excited  on  the  target  of  interest.  The  late-time  scattered  waveform 
then  appears  as  a  single  damped  sinusoid.  For  targets  significantly 
different  from  that  for  which  the  transmitted  waveform  was  designed,  the 
return  signal  is  distorted  due  to  the  presence  of  more  than  one  natural 
resonance.  In  this  manner,  the  desired  target  can  conceptually  be 
discriminated  from  the  others.  Of  course,  the  technique  implicitly 
requires  a  knowledge  of  the  natural  resonant  frequencies,  i.e.,  the  SEM 
poles,  in  order  to  synthesize  the  proper  transmitted  waveform,  but  this 
is  probably  not  a  restrictive  assumption. 

Numerical  results  for  the  radar  waveform  synthesis  method  have  been 
reported  for  the  case  of  a  high-Q  target  (a  thin  wire  at  both  normal 
[58]  and  oblique  [60]  incidence),  as  well  as  for  a  low-Q  target  (a 
sphere)  [59].  The  results  include  examples  of  the  synthesized  waveforms 
required  to  excite  either  the  first  or  third  natural  resonance  with 
either  a  zero  or  maximum  initial  value  for  the  late  time  return,  as  well 
as  illustrations  of  the  scattered  waveforms  for  both  the  desired  target 
and  targets  whose  characteristic  dimensions  differed  by  5  to  20  percent 
from  those  for  which  the  incident  waveform  was  designed  [58-60],  While 
the  latter  do  contain  some  deviation  from  a  single  damped  sinusoid, 
there  is  an  issue  as  to  whether  the  observed  differences  are  sufficient 
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to  provide  robust  discrimination  in  the  presence  of  noise.  This  concern 
is  further  complicated  by  the  lack  of  similar  data  for  more  complex 
targets  such  as  the  scale  model  aircraft  used  by  VanBlaricum.  For  these 
classes  of  targets,  the  differences  in  target  geometry  are  more 
pronounced  than  the  simple  scaling  considered  by  Chen,  et  al.,  in  his 
calculations,  and  hence  the  variations  in  the  return  waveform  may  be 
more  significant.  This  advantage  may  be  offset,  however,  by  the  lack  of 
a  precise  knowledge  of  the  location  of  the  SEM  poles  for  noncanonical 
scatterers.  The  poles  must  instead  be  obtained  through  experimental 
and/or  numerical  techniques  such  as  Prony's  method  as  described  in  the 
previous  section.  The  generation  and  radiation  of  the  required  radar 
waveform  is  also  a  potential  limiting  factor  in  implementing  this 
technique. 

This  latter  concern  has  been  addressed  recently  by  Rothwell ,  et 
al.,  [ (56]  and  Webb  and  Chen  [67].  They  suggest  a  method  wherein  a 
somewhat  arbitrary  waveform  is  transmitted  and  the  scattered  return  is 
processed  by  convolution  with  stored  versions  of  the  synthesized 
waveforms  for  single  mode  target  excitation.  This  process  effectively 
achieves  the  same  result  as  transmission  of  the  synthesized  waveforms 
themselves;  the  convolved  signal's  late-time  response  will  consist  of  a 
single  damped  sinusoid  provided  the  return  originated  from  the  desired 
target.  The  original  transmitted  waveform  need  only  contain  the 
appropriate  frequency  content  so  as  to  insure  excitation  of  all  relevant 
natural  resonances  with  sufficient  signal  power  in  the  presence  of 
noise.  It  is  important  to  note  that  this  frequency  content  could  be 
determined  by  a  process  similar  to  that  discussed  in  Section  3. 
Additional  considerations  in  the  practical  implementation  of  the  radar 
waveform  synthesis  method  are  presented  in  Section  4.3. 


V — I 

Jm 


4.1,3  The  K-PULSE  CONCEPT  -  Research  at  the  Ohio  State 
University  (OSU) 

One  of  the  earliest  and  most  significant  contributors  to  radar 
target  characterization  using  time-domain  scattering  data  was  E.  M. 
Kennaugh  and  others  at  Ohio  State  University  [61-64].  He  and  his 
colleagues  were  Instrumental  in  laying  the  groundwork  for  utilizing  both 
the  early-time  forced  response  and  late-time  natural  resonance  response 
of  the  total  time  domain  scattered  field  for  target  identification  and 
di scrimi nation . 

Research  regarding  the  former  centered  around  the  concept  of  the 
ramp  response  waveform  technique  for  target  shape  reconstruction 
[66,67].  Under  the  physical  optics  approximation,  it  was  shown  that  the 
field  scattered  by  a  target  as  a  function  of  time  when  illuminated  by  a 
ramp  waveform  is  directly  proportional  to  the  cross-sectional  a^ea  of 
the  target  as  a  function  of  range  over  the  illuminated  portion  of  the 
target  surface.  This  result  can  be  shown  to  be  a  time  domain 
representation  of  the  Bo jarski -Lewis  physical  optics  inverse  scattering 
theory  [68,69].  Hence,  a  reconstruction  of  the  target  profile  can  be 
obtained  by  collecting  ramp  response  waveforms  over  several  aspect 
angles  (provided  the  requirements  of  the  physical  optics  approximation 
are  satisfied).  3ecause  the  physical  optics  response  is  zero  when  the 
incident  field  is  zero,  it  follows  that  the  ramp  response  consists  only 
of  that  portion  of  the  signal  corresponding  to  the  time  during  which  the 
incident  ramp  waveform  passes  over  the  target  (early-time  response). 
This  is  not  particularly  a  disadvantage,  however,  because  a  large 
portion  of  the  total  energy  in  the  scattered  field  is  contained  in  the 
early-time  [65].  Experimental  demonstration  of  econstructions  of  the 
basic  shapes  of  military  aircraft  and  missiles  has  been  reported  [63]. 

As  Kennaugh,  Moffatt,  and  others  at  OSU  continued  their  work  in 
transient  scattering,  it  became  apparent  that  the  approximations 
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inherent  in  the  physical  optics  theory  limited  the  overall  applicability 
of  the  ramp  response  technique  when  used  strictly  as  described  above. 
In  order  to  extend  the  potential  of  their  methods  for  target 
discrimination  into  the  low-frequency  and  resonance  scattering  regimes 
so  as  to  include  the  late-time  response,  the  concept  of  natural 
resonances  (i.e.,  the  SEM)  and  the  geometrical  theory  of  diffraction 
(GTD)  were  brought  to  bear  on  the  problem.  By  integrating  the  ramp 
response  concept  with  the  waveform  predicted  from  a  knowledge  of  the 
natural  resonances  for  a  variety  of  targets,  a  predictor-correlator 
technique  was  formulated  from  which  target  discrimination  capability  was 
successfully  demonstrated  experimentally  [62,63].  The  method,  like 
VanBl ari cum1 s  proposed  classification  scheme,  requires  a  catalog  of  the 
complex  SEM  poles  in  order  to  predict  the  response  from  a  "candidate" 
target.  The  classification  is  achieved  by  comparing  the  measured  ramp 
response  with  the  predicted  waveforms  and  determining  the  normalized 
mean-square  error.  The  effects  of  noise  are  a  major  concern  [62]  and 
are  a  topic  of  further  research,  particularly  in  the  prediction  portion 
of  the  algorithm.  Similarly,  by  combining  the  ramp  response  with  GTD 
contributions  such  as  creeping  waves  and/or  edge  diffracted  fields,  the 
total  target  transient  response  can  more  accurately  be  approximated. 
Examples  for  spheres,  cones,  and  cones  spheres  have  been  reported  which 
validate  this  approach  [61]. 


More  recently,  Kennaugh  has  proposed  an  approach  to  target  discrim¬ 
ination,  based  almost  exclusively  on  the  SEM,  known  as  the  kill-pulse  or 
K-pulse  concept  [64].  This  technique  can  be  viewed  as  a  special  case  of 
the  radar  waveform  synthesis  method,  although  it  predates  Chen's  work  by 
several  years.  Like  the  radar  waveform  synthesis  method,  the  K-pulse 
concept  requires  illumination  (effectively)  of  the  target  with  a 
specific  incident  field  whose  complex  frequency  spectrum  is  chosen  to 
null  out  al  1  the  complex  poles  of  the  target's  response.  Thus,  the 
radar  return  consists  only  of  the  early- time  forced  response,  which  is 
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by  nature  of  finite  duration.  Because  several  such  Incident  waveforms 
exist,  the  K-pulse  Is  defined  as  the  Incident  waveform  of  minimum 
duration  having  the  above  properties. 

The  principal  advantage  of  K-pulse  excitation  Is  that  It  produces  a 
time-limited  response  from  the  desired  target  regardless  of  both  the 
aspect  and  range  ( 1 . e . ,  near-field  or  far-fleld)  of  the  target.  It  may 
therefore  prove  useful  In  situations  where  Insufficient  data  Is 
available  for  Imaging  or  other  more  data  Intensive  Identification 
schemes.  Like  many  of  the  other  techniques  described  In  the  previous 
section,  a  catalogue  of  the  complex  natural  resonances  of  targets  of 
Interest  is  required  for  the  K-pulse  concept.  The  number  of  data 
entries  required  for  each  target  is  small,  however,  because  of  the 
aspect-independence  of  the  complex  natural  resonances.  This  compares 
favorably  with  imaging  methods,  where  keys  must  be  developed  for  a  wide 
range  of  aspect  direction. 

4.2  "CARRIER- FREE"  RADAR  CONCEPTS 

A  second  area  of  research  in  wideband  radar  theory  and  applications 
has  been  conducted  almost  exclusively  by  H.  Harmuth  and  his  students  at 
Catholic  University  [32-52].  The  work  centers  around  the  use  of 
carrier-free  or  large  relative  bandwidth  signals  as  opposed  to  the 
conventional  radar  waveforms  which  impress  the  "useful"  signal 
Information  on  a  sinusoidal  carrier,  the  bandwidth  of  the  former  being 
only  a  few  percent  of  the  frequency  of  the  latter. 

Although  Harmuth' s  research  has  continued  for  over  a  decade,  it  has 
been  the  subject  of  intense  controversy  [53,54]  and  has  to  some  extent 
motivated  the  research  reported  herein.  A  survey  of  the  published  work 
in  carrier-free  radar  suggests  a  division  Into  two  principle  areas:  (1) 
utilization  of  carrier-free  radar  concepts  in  present-day  narrowband 


radar  applications  [33-46,52],  and  (2)  analysis  and  design  of  broadband 
components  for  the  generation,  transmission,  radiation,  reception,  and 
processing  of  carrier-free  radar  signals  [32,47-52].  A  review  and 
assessment  of  each  of  these  areas  are  presented  in  the  subsections  that 
follow. 

4.2.1  CARRIER-FREF  RADAR:  APPLICATIONS  STUDIES 

Harmuth's  earliest  papers  addressed  the  application  of  carrier-free 
radar  principles  to  a  series  of  generic  radar  and/or  radio  communication 
problems  [33-46], 

The  subject  of  spread-spectrum  communications  was  addressed  in  an 
article  [33]  which  suggested  the  use  of  periodic  nonsinusoidal  waveforms 
as  carriers  for  modulation  by  the  baseband  spread-spectrum  signals.  It 
was  shown  that  the  choice  of  broadband  carriers  leads  to  spectrum 
spreading  over  several  times  the  bandwidth  achievable  when  the  signal  is 
modulated  on  a  sinusoidal  carrier.  This  result  should  not  be  at  all 
surprising.  In  addition,  the  bandwidth  of  the  signal  is  on  the  order  of 
its  center  frequency,  i.e.,  the  transmitted  signal  is  of  large  relative 
bandwidth,  measured  in  terms  of  [33] 

fH  -  fL  (4-2) 
"  '  fH  *  fl 

where  f^  ^  are  the  upper  and  lower  frequency  limits  of  the  signal 
spectrum. 

Sinusoidal  carrier  systems  for  spread-spectrum  communications  are 
limited  by  hardware  considerations  (transmitters,  antennas,  and 
receivers)  to  values  in  the  range  of  n  <  0.1.  Thus,  large  spectrum 
spreadings  requires  a  correspondingly  high  carrier  frequency  for 
transmission.  Attenuation  due  to  atmospheric  losses  can  limit  the 


operating  distance  over  which  acceptable  communications  can  take  place. 
Harmuth  addresses  these  issues  in  a  cursory  manner  by  describing 
conceptual  hardware  implementations  of  receivers  which  selectively 
respond  to  nonsinusoidal  waves  of  a  particular  period  [32,33].  He 
argues  that  these  equipment  would  permit  the  use  of  large  relative 
bandwidth  signals  (having  n  approaching  unity)  which  would  not  suffer 
from  the  atmospheric  attenuation  present  at  higher  frequencies  [33].  It 
is  not  clear,  however,  how  the  distortions  introduced  by  dispersion  and 
frequency-dependent  attenuation  of  the  nonsinusoidal  carrier  signal 
would  affect  the  performance  of  the  proposed  receiver  implementation. 

Application  of  nonsinusoidal  or  carrier-free  signals  to  low-angle 
radar  tracking  is  the  subject  of  another  of  Harmuth's  articles  [34].  As 
discussed  in  Section  2.2,  the  principle  limitation  in  these  radar 
systems  are  the  effects  of  multipath  reflections  from  the  earth's 
surface  as  the  tracking  angle  gets  closer  to  the  horizon.  The  material 
in  Ref.  34  concentrates  primarily  on  addressing  this  limitation.  In 
that  article,  Harmuth  demonstrates  that  the  reflection  coefficient  of 
both  bare  ground  and  sea  water  at  near  grazing  incidence  is  very  nearly 
-1  from  100  MHz  to  10  GHz.  Similar  arguments  are  applied  to  the 
metallic  surface  of  targets.  With  these  results,  Harmuth  argues  that 
the  direct  reflection  will  be  reversed  in  polarity  with  respect  to  the 
transmitted  waveform,  while  the  multipath  signal  will  not.  He  then 
claims  that  the  use  of  a  nonsinusoidal  waveform  (a  video  pulse  of  about 
1  nsec  duration)  would  allow  this  polarity-reversal  to  be  used  to 
discriminate  between  the  two  returns  more  easily  than  a  conventional 
narrowband  waveform.  However,  the  basis  for  his  recommendation  of 
nonsinusoidal  waveforms  is  qualitative  and  tenuous  at  best.  Using 
simplistic  forms  for  the  return  from  a  complex  target  such  as  an 
aircraft  (see  Fig.  8,  [34]),  Harmuth  asserts  that  the  sum  of  the  direct 
and  indirect  return  results  in  a  more  "detectaDle"  signal  in  the 
presence  of  multipath.  It  is  evident  from  the  experimental  data  of  Van 
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Blaricum,  Fig.  4-3,  that  such  a  model  for  the  return  is  not  valid. 
While  Harmuth's  other  arguments  regarding  lower  atmospheric  losses  and 
reduced  noise  may  favor  nonsinusoidal  waveforms,  it  is  not  possible  to 
accept  its  use  for  low-angle  tracking  until  a  performance  analysis  of 
the  tracking  error,  similar  to  that  presented  in  Section  2.2.4  for 
narrowband  signals,  is  carried  out. 

Analytical  deficiencies  such  as  that  pointed  out  above  are  unfortu¬ 
nately  typical  of  many  of  Harmuth's  applications  oriented  articles. 
Often  a  great  deal  of  effort  is  spent  addressing  a  single,  sometimes 
secondary,  element  of  the  total  analysis,  while  the  remainder  is  left  to 
heuristic  qualitative  arguments.  A  second  example  of  this  occurs  in  a 
paper  advocating  the  use  of  large  relative  bandwidth  signals  for 
over-the-horizon  (OTH)  radar  [35].  Harmuth  properly  recognizes  the  fact 
that  OTH  systems  are  limited  to  signals  with  carrier  frequencies 
occupying  the  10  to  30  MHz  range,  thus  ensuring  reflection  off  one  of 
several  Ionospheric  layers.  It  is  also  true  that  in  order  to  achieve 
even  modest  range  resolution,  the  required  signal  bandwidths  are 
sufficiently  large  (>5  MHz)  such  that  the  relative  bandwidths,  Eq. 
(4-2),  are  not  characteristic  of  narrowband  waveforms  (n  <  0.05).  It  is 
well  known  that  signals  of  such  bandwidths  will  suffer  significant 
dispersion  upon  reflection  by  the  ionosphere,  along  with 
frequency-dependent  absorption  [53].  These  effects  will  severely 
distort  the  transmitted  waveform  and  make  reception  of  the  return  signal 
using  methods  based  on  a  priori  knowledge  of  the  incident  field  nearly 
impossible.  Rather  than  address  this  important  issue,  Harmuth  devotes 
nearly  the  entire  article  to  a  discussion  of  modulation  and  demodulation 
technique  for  signals  of  the  prescribed  relative  bandwidth  which  are 
founded  on  such  a  priori  information.  The  question  of  the  utility  of 
large  relative  bandwith  signals  for  OTH  radar  applications  remains 
therefore  largely  unanswered. 
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Harmuth  suggests  the  use  of  nonsinusoidal  signals  for  synthetic 
aperture  radar  (SAR)  applications  in  a  series  of  ten  papers  [36-45], 
each  addressing  a  particular  aspect  of  the  SAR  concept.  Unfortunately, 
his  treatment  of  these  topics  suffers  from  two  serious  shortcomings:  a 
misunderstanding  of  the  nature  of  conventional  SAR  theory  and 
processing,  and  an  emphasis  on  processing  techniques  for  nonsinusoidal 
SAR  which  are  insufficiently  robust  to  noise  and/or  signal  distortion 
for  use  in  practical  system  implementations.  It  is  regretable  that 
these  limitations  obscure  what  may  be  a  valid  application  for  large 
relative  bandwidth  signals  (although  not  necessarily  the  signals 
suggested  by  Harmuth).  Indeed,  the  improved  performance  of  target 
detection  in  clutter  suggested  by  the  results  in  Section  3  for  low 
frequency  carrier  (and  hence  large  relative  bandwidth)  signals  would 
certainly  pertain  to  SAR  applications. 

Harmuth's  most  serious  misconception  regarding  the  principles  of 
SAR  is  his  statement  that  a  SAR  system  relies  upon  the  existence  of  a 
Doppler  shift  due  to  relative  motion  between  the  sensor  and  the  scene  in 
order  to  synthesize  a  large  aperture  from  which  improved  resolution  is 
achieved.  He  (correctly)  suggests  that  a  nonsinusoidal  cAR  radiating  a 
narrow  pulse  can  achieve  similar  improvements  without  the  need  for 
exploiting  any  Doppler  effects  [36,37].  From  this  he  concludes  that  the 
standard  SAR  configuration  of  a  single  transmi tter/recei ver  which 
synthesizes  an  aperture  by  radiating  pulses  periodically  as  it  moves 
along  the  flight  path  can  be  replaced  by  a  stationary  array  of 
transmi tter/recei  vers  located  at  the  points  where  the  pulses  are 
radiated.  It  is  further  stated  that  with  such  an  array,  Doppler  shifts 
could  be  used  to  detect  target  motion,  since  they  are  not  being  used  to 
form  the  synthetic  aperture.  All  these  conclusions  are  indeed  true. 
But,  as  stated  above,  they  are  also  true  for  a  SAR  system  radiating  a 
conventional  "chirped"  signal.  While  the  use  of  Doppler  concepts  are 
convenient  in  conveying  the  principles  of  synthetic  aperture  radar  [70], 
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It  Is  well  known  that  the  Doppler  effects  which  occur  during  the 
reception  of  a  single  pulse  must  be  neglected  In  order  to  properly 
formulate  the  SAR  theory  [71,  Eq.  (4)].  Thus,  a  system  which  could 
achieve  a  measurement  of  the  scattered  pulses  from  several  points  along 
a  synthetic  aperture  while  remaining  stationary  would  actually  perform 
better  than  one  In  motion.  This  could  obviously  be  achieved  using 
Harmuth's  stationary  array.  Such  an  implementation  has  in  fact  been 
suggested  by  Farhat  [72]  using  a  conventional  chirped  waveform. 

It  is  also  worth  noting  that  Harmuth's  claim  that  such  an  array 
would  require  an  element  spacing  of  half  a  wavelength,  while  the  spacing 
of  a  nonsinusoidal  array  would  be  much  greater,  (implying  less  elements 
and  hence  less  complexity  [36]),  is  also  false.  Both  arrays  would  have 
spacings  dictated  by  the  rate  of  change  of  target  distance  from  one 
element  to  another,  so  as  to  unambiguously  sample  the  entire  scene  as 
defined  by  the  beamwidth  of  the  individual  array  elements.  In  both 
cases,  this  distance  is  many  times  the  wavelength  of  the  carrier  and  is 
given  by  the  limits  on  the  conventional  SAR  pulse  repetition  frequency 
[73].  It  can  thus  be  concluded  that  aperture  synthesis  can  be  obtained 
from  any  system  capable  of  measuring  distance  to  a  target  as  a  function 
of  aspect  direction.  Certainly  a  radar  using  either  a  chirped  signal  or 
a  nonsinusoidal  pulse  satisfies  this  requirement.  Therefore,  the 
arguments  relating  to  the  use  (or  non-use)  of  Doppler  shifts  made  by 
Harmuth  cannot  justify  a  preference  for  nonsinusoidal  waves  in  SAR 
applications. 

A  second  claim  which  Harmuth  expounds  in  several  papers  [36-45]  in 
the  "Nonsinusoidal  SAR"  series  is  that  the  use  of  nonsinusoidal  waves 
leads  to  angular  resolution  which  depends  both  on  the  signal  bandwidth 
and  the  radiated  power;  a  property,  he  claims,  which  is  not  true  of 
narrowband  SAR  systems.  There  are,  however,  several  fallacies  in  his 
arguments. 
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To  show  this,  consider  a  classical  SAR  operating  at  a  wavelength  xQ 
with  a  bandwidth  per  pulse  of  Af.  It  is  Important  to  distinguish  the 
difference  between  range  and  cross-range  resolution,  given  by  [70] 


k  c 
Pr  "  2Af 


(4-3) 
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respectively,  where  k  >  1,  R  is  the  range,  and  L  is  the  synthetic 

r  >  a 

aperture  length.  Note  that  Pr  exhibits  the  well  known  inverse 
dependence  on  bandwidth,  while  pa  is  achieved  by  coherently  tracking  the 

a 

phase  of  the  carrier,  hence  the  dependence  or.  xQ. 


For  Harmuth's  system,  a  nonsi nusoi dal  pulse  of  duration  AT  is 
assumed  to  be  radiated.  Using  a  very  idealistic  model  for  the  receiver 
transfer  function,  the  received  signal  is  assumed  to  be  a  triangular 
pulse  of  duration  2aT.  The  classical  range  resolution  from  such  a  pulse 
is  simply 


k  cAT  k  c 
_r _  r 

2  =  4Af 


(4-5) 


where  Af  is  the  effective  one-sided  bandwidth  of  the  receiver  lowpass 
transfer  function.  Note  that  Eqs.  (4-3)  and  (4-5)  are  in  complete 
agreement  provided  that  the  two-sided  lowpass  bandwidth  is  used  in  Eq. 
(4-5).  Since  no  carrier  is  available  for  tracking,  the  achievable 
cross-range  resolution  is  a  result  of  triangulation  of  the  return  for 
various  points  along  the  synthetic  aperture  length  as  a  function  of 
range.  This  can  be  shown  to  lead  to  a  cross-range  resolution  of  the 
form 
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which  again  is  completely  analogous  to  Eq.  (4-4).  This  is  also  con¬ 
sistent  with  Brown's  interpretation  of  SAR  in  the  Fourier  domain  [74]. 
Eqs.  (4-4)  and  (4-6)  express  the  well  known  result  that  the  azimuth 
resolution  of  a  SAR  improves  with  increasing  carrier  frequency;  for  a 
lowpass  (no  carrier)  system,  It  is  limited  by  the  maximum  frequency  of 
the  signal.*  Thus,  while  the  use  of  nonsinusoidal  waves  provides 
comparable  range  resolution  to  a  classical  narrowband  SAR  provided  the 
same  absolute  bandwidth  is  transmitted,  it  will  in  general  produce 
degraded  cross-range  resolution,  since  typically  Af <  f  . 

In  order  to  overcome  this  limitation  associated  with  conventional 
processing  of  unconventional  (i.e,,  nonsinusoidal)  waveforms,  Harmuth 
advocates  the  use  of  "slope  processing"  of  the  returned  signals.  Under 
the  assumption  that  the  scattered  signal  is  the  sum  of  many  triangular 
pulses  with  different  amplitudes  and  delays,  he  shows  that  measurement 
of  the  slopes  of  these  pulses  leads  to  Improvement  over  the  resolution 
given  Eq.  (4-6).  Specifically,  he  asserts  that 


2k  cR 


where  yrP/Pn  is  the  voltage  signal-to-noise  ratio  (SNR)  at  the  receiver 
[36,44].  While  such  a  scheme  may  be  valid  theoretically  under  rather 
idealistic  assumptions,  it  is  not  robust  enough  to  be  considered  for 
practical  implementation.  This  is  because  the  unknown  distortion  of  the 
pulse  Introduced  by  antenna  and  atmospheric  dispersion  as  well  as  the 
complex  scattering  by  the  target  (c.f.  Fig.  4-3)  makes  any  processing 


*This  result  is  also  observed  in  computed  tomography  [75]. 
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which  depends  on  a  specific  pulse  shape  Invalid.  Furthermore,  It  Is 
likely  that  for  Imaging  applications  yP/P  <  1,  making  the  use  of  Eq. 
( 4“ 7 )  a  moot  point.  Even  under  the  circumstances  for  which  slope 
processing  is  valid  ( 1 . e . ,  high  SNR  and  Ideal  pulse-shape  preservation), 
the  Improvements  suggested  by  Eq.  (4-7)  are  not  restricted  to 
nonsinusoldal  waveforms.  Similar  tradeoffs  exist  between  bandwidth  and 
SNR  in  the  classical  SAR  for  a  variety  of  superresolution  techniques. 
All  lack  the  robustness  required  for  operation  under  real  world 
conditions,  particularly  for  distributed  targets,  which  is  the  typical 
mode  of  operation  for  a  SAR  when  producing  a  radar  image. 

In  fairness,  it  Is  worth  mentioning  a  point  made  by  Harmuth  early 
in  his  treatment  of  SAR  [36],  wherein  he  states  that  the  theoretical 
resolution  limit  of  a  SAR  can  be  achieved  for  all  points  at  arbitrary 
ranges  provided  one  uses  the  differences  of  squares  of  time  (o 
distance)  Instead  of  the  differences  of  time  alone.  This  is  an 
important  observation,  and  is  equivalent  to  reconstruction  by  back- 
projection  along  spherical  surfaces  [76]  instead  of  straight  lines  [77]. 
Such  a  method  is,  however,  extremely  computationally  intensive  and  is 
hence  less  preferred  than  Fourier  domain  methods  for  range-curvature 
correction  [78].  It  is  also  true  that  such  processing  can  be  applied  to 
conventional  as  well  as  nonsinusoidal  signals,  and  hence  does  not 
provide  justification  for  use  of  the  latter. 

The  final  application  of  nonsinusoidal  waves  expounded  by  Harmuth 
is  their  use  in  "anti-stealth  radar,"  i.e.,  for  detection  of  targets 
treated  with  absorbing  materials  for  reduction  of  their  radar  cross- 
section.  In  the  subject  paper  [46],  Harmuth  correctly  observes  that 
most  radar  absorbing  material  (RAM)  coatings  are  restricted  to 
relatively  high  frequencies  due  to  thickness  limitations  (for  ohmic  RAM) 
or  to  small  relative  bandwidths  (for  tuned  or  resonant  RAM). 
Furthermore,  the  application  of  RAM  coatings  is  generally  localized  to 


portions  of  a  target  associated  with  prominent  "scattering  centers"  such 
as  edges,  corners,  etc.  This  concept  Is  by  nature  valid  only  at 
wavelengths  much  shorter  than  the  target  dimensions  (i.e.,  It  Is  a  high 
frequency  model).  As  a  result,  the  localized  RAM  offers  little 
potential  for  RCS  reduction  when  the  wavelength  Is  Increased  to  match 
the  size  of  the  target;  In  this  case,  the  target  resonates  as  a  whole. 
In  fact,  these  resonances  are  the  basis  of  the  SEM  concepts  discussed  In 
Section  4.1. 

It  is,  therefore,  entirely  appropriate  for  Harmuth  to  conclude  that 
radar  signals  with  large  relative  bandwidth  are  strong  candidates  for 
use  in  detecting  RAM-treated  targets.  While  the  conclusions  are  valid, 
however,  the  reasons  presented  for  justifying  them  are  at  best  overly 
simplistic.  Specifically,  Harmuth  discusses  the  situation  of  an 
infinite  metallic  half  space  coated  with  a  lossy  dielectric/magnetic 
material  and  illuminated  at  normal  Incidence  by  a  baseband  video  pulse 
of  finite  duration.  His  arguments  proceed  as  follows.  By  properly 
choosing  the  thickness  and  material  properties  of  the  coating,  the 
reflected  waveforms  from  the  front  and  back  surfaces  can  be  made  tc 
cancel  exactly  for  a  sinusoidal  signal  at  a  given  frequency.  For  a 
nonsinusoidal  signal,  on  the  other  hand,  the  front  and  back  surfaces 
would  produce  independent  reflections  which  would  not  overlap  provided 
the  pulse  was  of  sufficiently  short  duration.  While  all  this  is  true, 
it  is  also  irrelevant  in  practical  applications.  Even  if  the  pulse 
shape  could  be  preserved  after  propagation  through  the  atmosphere,  it  is 
not  possible  to  ignore  the  fact  that  a  complex  target  would  not  preserve 
the  shape,  and  more  importantly,  the  duration,  of  the  incident  pulse 
(once  again,  see  the  waveforms  of  Fig.  4-3).  This  is  because  the 
incident  wave  is  being  scattered  continuously  from  points  along  the 
range  dimension  of  the  target,  and  hence  would  consist  of  the 
superposition  of  a  continuously  delayed  replicas  of  the  incoming  signal. 
This  interpretation,  which  accounts  from  the  early-time  return  (see 


Section  4.1),  Is  further  complicated  by  the  late-time  resonances  which 
follow.  Thus,  even  if  one  could  describe  the  scattered  waveform  as  the 
sum  of  returns  from  the  Inner  and  outer  layers  of  the  coating  (this 
Itself  Is  questionable),  a  system  which  relies  on  their  "disjointness" 
in  time  would  fall.  What  Is  amazing,  as  mentioned  above,  Is  that 
despite  the  shortcomings  In  Harmuth's  reasoning,  his  conclusions 
regarding  the  utility  of  large  relative  bandwidth  signals  for  detecting 
low  observable  targets  Is  still  essentially  valid. 

4.2.2  CARRIER- FREE  RADAR:  COMPONENTS  STUDIES 

In  a  more  recent  series  of  articles  (save  one),  Harmuth  addresses 
the  realization  of  various  critical  components  required  for  operating  of 
nonsinusoldal  radar  systems,  namely  transmitting  and  receiving  antenna 
elements  and  arrays  [47-49],  selective  receivers  [32],  waveguides  [50], 
and  resonant  cavities  [51].  Each  of  these  topics  is  discussed  below. 

The  need  for  efficient,  directive  and  distortionless,  transmitting 
and  receiving  antennas  is  perhaps  the  most  important  requirement  in 
implementing  practical  nonsinusoidal  radar  systems.  While  a  great  deal 
of  research  has  been  performed  over  the  years  on  frequency  independent 
antennas  (see  Chapter  6,  [79],  for  example),  the  definition  and  analysis 
of  radiating  systems  specifically  for  large  relative  bandwidth  video 
signals  is  essentially  new  [80,81],  Harmuth's  work  in  this  area  has 
emphasized  the  use  of  a  radiating  element  which  he  calls  a  "large 
current  radiator  (LCR),"  shown  in  Figure  4-4  [47],  The  radiator  is 
designed  to  operate  much  like  a  short  electric  dipole,  and  as  such, 
radiates  at  a  single  frequency  with  the  familiar  figure-eight  elevation 
pattern  and  isotropic  azimuth  pattern  characteristic  of  dipole  current 
elements.  Directive  radiation  patterns  are  then  achieved  by  combining 
several  LCRs  into  an  array  [49].  The  differences  between  the  LCR  and 
conventional  short  electric  dipole  are  an  attempt  to  overcome  the 
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Figure  4-4.  Schematic  Diagram  of  the  Large  Current 
Radiator  (after  [47]). 


limitations  of  toe  latter  for  radiation  of  low  frequency  signals. 
Because  the  electric  dipole  is  by  nature  a  capacitive  structure,  it 
cannot  support  large  currents  at  low  frequencies,  even  with  top  loading, 
since  the  current  literally  has  "nowhere  to  go"  [47-49].  In  order  to 
Increase  the  current  flowing  along  the  dipole  at  low  frequencies,  the 
LCR  feeds  the  dipole  at  its  endpoints,  and  thus  "closes  the  loop," 
giving  the  current  a  place  to  flow.  Should  the  return  path  of  the  loop 
be  symmetric  with  the  dipole  element,  however,  the  total  structure  will 
radiate  in  the  less  efficient  quadrapole  mode.  In  order  to  emphasize 
the  dipole  mode  of  operation,  the  loop  Is  intentionally  made  asymmetric 
by  first  isolating  the  return  path  from  the  dipole  via  a  conducting 
shield,  and  then  coating  the  shield  with  absorbing  material  to  eliminate 
the  "image"  of  the  dipole  formed  by  the  conductor  (see  Figure  4-4). 
Further  asymmetry  can  be  introduced  by  using  different  conductor  geome¬ 
tries  for  the  dipole  and  the  return  path.  Arrays  of  LCRs  of  this  type 
have  been  successfully  used  to  radiate  and  receive  nonsinusoidal  signals 
for  into-the-ground  radar  probing  applications  [52]. 

Harmuth's  analysis  of  the  LCR  attempts  to  address  its  performance 
as  a  transmitter  and  receiver  [47,48],  both  alone  and  in  arrays  [49],  in 
terms  of  standard  measures  for  sinusoidal  time  dependence  such  as  radia¬ 
tion  resistance,  input  impeoance,  gain,  directivity,  efficiency,  and  so 
on.  Each  of  these  quantities,  under  the  conventional  definition,  is 
frequency  dependent,  and  hence  for  time-varying  signals  of  nonvanishing 
bandwidth,  each  becomes  a  time- varying  quantity.  As  a  result,  it  is 
usually  appropriate  to  define  them  in  an  average  sense  (with  respect  to 
time),  so  as  to  remove  any  dependence  on  the  instantaneous  variation  of 
the  waveform  [80]. 

It  is  well  known  from  Maxwell's  equations  that  any  antenna  immersed 
in  a  linear  stationary  medium  can  be  modeled  as  a  linear,  time-invariant 
system.  Thus,  the  output  of  an  antenna,  for  a  given  input,  is 


determined  strictly  from  a  knowledge  of  the  transfer  function  (in  the 
frequency  domain)  or  the  impulse  response  (in  the  time  domain)  of  the 
system.  Reciprocity  guarantees  that  this  fact  holds  regardless  of 
whether  the  antenna  is  used  for  transmission  or  reception.  (For 
details,  see  Appendix  C.)  Once  the  transfer  function  or  Impulse 
response  is  known,  all  of  the  above-mentioned  performance  parameters 
(ga'n,  etc.)  can  be  determined.  For  wire-like  antennas,  such  as  the 
LCR,  the  transfer  function  is  most  straightforwardly  defined  in  terms  of 
the  current  distribution  along  the  wire.  While  it  is  not  possible  in 
general  to  determine  analyt1call>  the  exact  current  distribution  in 
closed  form,  there  exist  well  established  "exact"  numerical  methods 
(such  as  the  method  of  moments)  for  evaluating  the  current  in  either  the 
frequency  [82]  or  the  time  [83]  domains.  It  then  follows  that  the 
antenna  performance  parameters  can  likewise  be  numerically  calculated. 


Rather  than  apply  this  standard  method  of  analysis,  Harmuth  uses  a 
rather  unconventional  approach  in  his  examination  of  the  LCR.  In  an 
attempt  to  provide  an  analytical  characterization  of  the  performance  of 
the  LCR,  Harmuth  resorts  to  a  series  of  approximations  which  are 
justified  at  best  heuri stically,  thus  leaving  many  of  his  results  open 
to  question.  He  begins  by  assuming  the  current  is  spatially  uniform 
along  the  length  of  the  LCR  [47].  This  is  probably  a  reasonable 
assumption  for  short  wire  lengths  and  symmetric  excitation.  Harmuth's 
second  approximation  is  critical  to  his  results,  and  at  the  same  time, 
subject  to  debate.  First,  a  specified  antenna  current  1(t)  is  assumed 
to  flow  along  the  LCR.  The  instantenous  power  delivered  to  the  far  zone 
is  shown  to  be 
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The  essense  of  the  approximation  is  to  equate  Pj  with  the  instantaneous 
power  at  the  antenna  terminals,  i.e., 
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P^t)  *  u(t)i(t),  (4-9) 

where  u(t)  Is  the  terminal  voltage.  Given  i(t),  Pj  cen  be  obtained  from 
Eq.  (4-8)  and  u(t)  from  Eq.  (4-9).  By  modeling  u(t)  as  being  generated 
by  a  Norton  equivalent  current  source  with  driving  current  iQ  and 
internal  resistance  R^ ,  it  follows  that  the  required  time-varying  source 
current  iQ(t)  can  be  found. 

Unfortunately,  this  procedure  comes  under  scrutiny  as  to  the  vali¬ 
dity  of  Eq.  (4-9),  which  neglects  the  contributions  to  the  terminal 
voltage  from  ohmic  losses  in  the  wire  and  reactive  fields  in  the  near 
zone  of  the  antenna.  Harmuth  presents  heuristic  arguments  for 
neglecting  the  latter  when  the  current  i(t)  has  a  sufficiently  rapid 
rise  time,  but  this  begs  the  question  since  the  required  rise  time  may 
be  limited  by  the  neglected  portions  of  the  driving  voltage,  i.e.,  by 
the  terminal  inductance  of  the  ICR  wire  loop.  The  ohmic  losses  are 
completely  ignored  by  Harmuth,  and  yet  it  is  known  they  represent  a 
significant  portion  of  the  input  resistance  of  a  conventional  dipole 
[79].  Whether  the  same  is  true  for  the  LCR  has  yet  to  be  determined. 

The  use  of  the  LCR  as  the  fundamental  building  block  in  an  antenna 
array  for  directive  radiation  of  nonsinusoidal  waveforms  is  the  topic  of 
Ref.  49.  Harmuth  defines  the  directivity  in  terms  of  peak  amplitude, 
peak  power  (amplitude  squared),  and  energy  (integrated  power)  patterns, 
the  latter  being  the  most  robust  concept  and  one  which  has  found 
acceptance  elsewhere  [80].  Furthermore,  from  a  target  detection  point 
of  view,  the  energy  pattern  is  the  most  relevant  as  shown  in  the 
analysis  presented  in  Section  3  of  this  report.  It  is  important  to 
point  out  a  specific  definition  of  directivity  is  valid  only  to  the 
extent  that  it  is  required  by  a  more  fundamental  performance  speci¬ 
fication  such  as  detection,  resolution,  or  estimation.  It  is  in  this 
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context  that  a  fourth  directivity  measure,  the  slope  pattern,  introduced 
by  Harmuth  as  a  means  of  improving  the  directivity,  must  be  questioned. 
The  ability  to  measure  slope,  as  mentioned  in  Section  4.2.1,  is 
extremely  sensitive  to  noise  ar.d  waveform  distortion,  and  hence  is  not 
likely  to  provide  a  robust  detection,  resolution,  and/or  estimation 
capability. 

Despite  these  shortcomings,  the  experimentally  proven  utility  of 
the  LCR  for  short  distance,  into-the-ground  probing  justifies  its 
consideration  as  a  candidate  transmitting  and  receiving  antenna  for 
nonsinusoidal  radar  applications,  provided  further  analysis  supports  the 
preliminary  results  demonstrated  to  date  [47-49,52]. 

Another  set  of  radar  system  components  investigated  by  Harmuth  for 
operation  with  large  relative  bandwidths  is  the  lossless  waveguide  and 
waveguide  cavity  resonator  [50,51].  These  are  important  elements  for 
the  transmission  and  selective  reception,  respectively,  of  sinusoidal 
waveforms  in  conventional  radar  systems.  In  what  is  perhaps  the  most 
significant  result  of  his  research  in  nonsinusoidal  waves,  Harmuth  has 
shown  that  certain  signals  with  large  relative  bandwidth  can  propagate 
and  resonate  di s tort ion- free  in  rectangular  waveguides  [50]  and  cavity 
resonators  [51].  This  finding  goes  against  the  intuitive  view  of 
waveguides  and  resonators  as  dispersive  components,  and  opens  up  the 
possibility  of  extremely  simple  structures  for  the  transmission  and 
selective  discrimination  of  periodic,  nonsinusoidal  waves,  an  issue 
which  has  been  raised  in  arguments  against  the  adoption  of  such 
waveforms  [53].  Harmuth 's  result  is  based  on  the  fact  that  there  exist 
a  denumerably  infinite  number  of  conventional  waveguide  modes  and 
harmonically  related  frequencies  for  which  the  waveguide  impedance  and 
phase  velocity  are  independent  of  frequency.  Any  superposition  of  these 
mode-frequency  combinations  constitute  a  valid  distortion-free  propa¬ 
gating  signal.  because  the  frequencies  are  harmonically  related,  it 
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follows  that  the  resulting  signal  is  periodic  with  period 

T  -  (4-10) 

While  the  significance  of  this  discovery  cannot  be  downplayed,  there 
remains  additional  work  to  be  done  before  the  utility  of  these 
components  in  nonsinusoidal  radar  applications  is  fully  established. 
Most  importantly  is  a  need  to  establish  the  information  "bandwidth"  that 
can  be  propagated  with  each  of  the  nonsinusoidal  modes  in  the  waveguide 
at  some  level  of  acceptable  distortion,  in  the  same  way  the  conventional 
modulation  is  used  to  convey  information  on  each  of  the  sinusoidal 
modes.  Inat  is,  if  the  nonsinusoidal  modes  are  perturbed  slightly 
through  some  modulation  scheme  (see  [39],  for  example),  does  the 
resulting  waveform  still  propagate  without  unacceptable  dispersion  down 
the  guide?  Excitation  and  detection  of  nonsinusoidal  signals  in  the 
waveguides  and  cavity  resonators  must  also  be  addressed. 

Having  reviewed  two  principal  areas  of  research  in  wideband  radar 
concepts  (the  SEM,  Section  4.1,  and  carrier-free  or  nonsinusoidal  radar. 
Section  4.2),  the  next  section  evaluates  these  concepts  as  to  their 
technical  merit  and  potential  for  improvement  over  conventional, 
narrowband  radar  systems.  Practical  limitations  in  implementing  the 
preferred  approaches  are  discussed. 

4.3  COMPARISON  TO  NARROWBAND  RADAR  SYSTEMS 

It  is  evident  from  the  discussions  in  the  previous  sections  that 
the  use  of  wideband  radar  signals  offers  potential  for  improvement  in 
many  aspects  of  radar  system  performance.  It  is  also  true  that  there 
exist  certain  practical,  and  perhaps  even  fundamental,  limitations  on 
the  implementation  of  radar  systems  using  such  signals  with  present  day 
technology. 


The  discussions  In  Section  2  of  this  report  suggest  that  Improve¬ 
ments  In  the  performance  of  conventional  narrowband  radar  systems  could 
be  obtained  over  a  wide  range  of  applications  through  advances  In  the 
following  generic  areas: 

1.  Increased  system  slgnal-to-nol se  ratio  (SNR). 

2.  Finer  range  and  angular  resolution  at  all  ranges. 

3.  Clutter  suppression. 

4.  Aspect  and  Illumination  Independent  target  classification. 

The  potential  of  wideband,  1.e.»  large  relative  bandwidth,  signals 
for  achieving  advances  in  nearly  all  of  the  above  areas  has  been 
identified  in  various  places  throughout  this  report.  The  signal 
scattered  by  a  target,  particularly  one  treated  with  RAM,  can  be  greatly 
enhanced  by  operation  in  the  resonance  region  of  the  target  scattering 
domain,  particularly  with  signals  spanning  several  octaves.  This  fact 
has  been  recognized  In  both  the  SEM  (see  Section  3)  and  nonsinusoldal 
radar  [46]  research.  Harmuth  has  also  pointed  out  that  the  external 
noise  in  the  0.5  -  10  GHz  region  is  considerably  lower  than  at  higher 
frequencies  [52],  making  operation  in  this  region  desirable,  all  other 
considerations  being  equal.  These  two  facts  suggest  improved  SNR  as  a 
significant  consequence  of  radar  systems  operating  at  frequencies  below 
10  GHz.  Section  3  showed  that  enhanced  clutter  suppression  is  also 
possible  with  operation  at  lower  frequencies. 

Inasmuch  as  range  resolution  is  fundamentally  determined  by  the 
signal  bandwidth,  it  becomes  clear  that  fine  range  resolution  and  low 
frequency  operation  imply  the  use  of  large  relative  bandwidth  signals. 
As  discussed  in  Section  4.2.1,  angular  resolution,  on  the  other  hand,  is 
governed  primarily  by  the  absolute  frequency  (not  bandwidth)  of  the 
signal,  whether  real  or  synthetic  aperture  areas  are  used  (neglecting 
non- robust  superresolution  methods).  This  suggests  a  tradeoff  between 
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relative  bandwidth  and  the  need  for  both  fine  angular  resolution  and 
high  SNR  and/or  low  clutter.  It  should  be  pointed  out  that  low  angle 
tracking,  while  often  considered  to  be  limited  by  angui&r  resolution, 
can  also  be  improved  through  better  range  resolution  because  of  the 
delay  In  the  multipath  signal. 

Finally,  regarding  point  (4)  above,  perhaps  the  single-most 
emphasized  advantage  of  the  SEM  approach  is  the  ability  to  completely 
characterize  the  scattering  from  a  target  with  a  finite  number  of 
parameters  (the  SEM  poles)  regardless  of  aspect  cr  type  of  illumination. 
Of  course,  the  excitation  of  many  poles  with  significant  energy  requires 
a  broadband  waveform  center  about  the  resonant  frequency  range  of  the 
target,  again  suggesting  the  use  of  large  relative  bandwidth  signals. 

Given,  then,  the  obvious  potential  of  wideband  radar  signals  for 
improved  performance,  it  is  necessary  to  quantitatively  assess  the 
extent  to  which  this  potential  can  be  realized  in  a  particular  radar 
application.  This  can  only  be  achieved  with  an  end-to-end  systems 
analysis  so  that  the  tradeoffs  between  advantages  and  disadvantages  can 
be  fully  ascertained.  Such  an  analysis  for  the  problem  of  target 
detection  in  clutter  plus  noise  was  performed  in  Section  3.  The  results 
of  that  analysis  clearly  show  the  variability  in  the  optimum  waveform 
for  differing  system  constraints.  They  also  allow  the  effects  of 
improvements  in  wideband  components  such  as  antennas,  receivers,  etc., 
to  be  incorporated  in  the  performance  prediction.  These  components,  as 
discussed  later  in  this  section,  are  a  major  practical  limitation  in 
developing  wideband  radar  systems.  Perhaps  the  most  important 
observation  to  be  made  from  Section  3  is  that  one  cannot  simply  argue 
for  wideband  versus  narrowband  signals;  instead,  each  application  must 
be  addressed  individually. 


Regarding  the  use  of  the  SEM  for  target  classification,  Van 
Blarlcum  and  others  readily  admit  the  need  for  further  studies  before 
Its  promising  potential  can  be  realized  [65].  In  addition  to  Improve¬ 
ments  In  the  classification  algorithm  regarding  pole  extraction  and 
recognition  and  the  inclusion  of  early-time  response  data.  Van  Blaricum 
clearly  states  that  an  end-to-end  system  analysis  Is  required  to  address 
practical  issues  such  as  the  operating  scenario,  component  design,  etc. 

Harmuth,  on  the  other  hand,  is  less  willing  to  acknowledge  the  need 
for  further  extensions  of  his  cursory  studies  of  nonslnusoidal  radar 
applications  [32-46].  His  simplifying  and  often  piecemeal  approach, 
coupled  with  his  tendency  to  consider  only  particular  types  of 
waveforms,  makes  it  difficult  to  assess  the  performance  gain  which  might 
be  obtained  in  each  of  the  applications  he  considers  should  a  more 
rigorous  analysis  be  used.  It  is  important  to  emphasize,  however,  that 
despite  the  shortcomings  of  his  analyses,  the  applications  suggested  by 
Harmuth  and  discussed  in  Section  4.2.1  may  still  be  enhanced  through  the 
use  of  wideband  waveforms. 

A  major  issue  in  implementing  radar  systems  with  large  relative 
bandwidth  is  the  existence  of  components  capable  of  operating  over  wide 
ranges  of  frequencies.  Harmuth' s  treatment  of  the  principle  limiting 
factors,  namely  transmission,  radiation,  sensing,  and  selective 
reception  of  nonsinusoidal  waves  [47-51],  is  much  more  satisfying  than 
his  work  in  WBR  applications.  This  is  an  important  distinction,  since  a 
component  concept  developed  by  Harmuth  could  prove  useful  in 
applications  expounded  elsewhere,  such  as  using  the  SEM  for  target 
classification.  While  his  work  in  NSR  components  also  requires 
additional  refinement,  Harmuth's  orginal  contributions  in  the  areas  of 
antennas,  waveguides,  and  cavity  resonators  represent  a  significant 
first  step. 
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It  is  quite  common  in  discussions  regarding  the  potential  of 
wideband  radar  systems  to  call  out  "fundamental  limitations"  inherent  in 
the  use  of  signals  having  large  relative  bandwidths  [53,54].  One  must 
carefully  distinguish,  however,  those  limitations  which  are 
"fundamental"  from  those  which  are  "practical",  l.e.,  those  limited  by 
current  state-of-the-art  technology.  Perhaps  the  most  significant 
fundamental  limitation  In  the  use  of  wideband  radar  signals  is  the 
distortion  Introduced  by  (unknown)  dispersion  and  frequency-dependent 
attenuation  due  to  the  atmosphere,  ionosphere,  and  to  a  lesser  extent, 
terrain.  Any  WBR  system  concept  must  be  sufficiently  robust  to  tolerate 
this  distortion,  either  by  minimizing  it  a  priori  or  correcting  for  it  a 
posteriori.  A  second,  potentially  "fundamental"  limitation  concerns  the 
efficient  radiation  and  reception  of  baseband  (or  nearly  baseband) 
signals.  While  it  is  true  that  radiation  efficiency  depends 
fundamentally  on  the  acceleration  rate  (and  hence  frequency)  of  charges, 
it  is  not  necessarily  true  that  a  reasonably  efficient  radiator  at  lew 
frequencies  could  not  replace  a  more  efficient  high  radiator  should  t.he 
difference  be  made  up  in  reduced  noise  or  increased  sca'tering  from  the 
target.  What  is  currently  needed  is  a  quantification  of  just  how 
"reasonably  efficient"  a  radiator  for  large  relative  bandwidths  can  be. 
To  this  end,  it  is  Important  to  investigate  not  just  conventional 
antenna  designs,  intended  for  signals  with  small  instantaneous 
bandwidths,  but  to  consider  unconventional  antennas  whose  design  is 
motivated  specifically  by  the  desire  to  radiate  wideband  waveforms  (such 
as  Harmuth's  LCR). 

The  above  discussion  is  not  meant  to  suggest  that  "practical" 
limitations  are  of  lesser  significance.  Nonetheless,  it  is  true  that  in 
today's  technological  environment,  a  concept  which  may  not  be 
practically  implemented  in  the  near  term  should  not  necessarily  be 
disregarded  for  long  term  application,  particularly  when  the  potential 
payoff  is  high.  What  is  important  is  to  employ  sound  research  methods 
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to  both  quantify  this  payoff  and  Identify  those  critical  technologies 
where  Improvements  are  required.  Only  then  can  the  cost  versus  risk 
yersus  payoff  be  confidently  established.  Unfortunately,  such  an 
approach  has  not  always  been  pursued;  Instead,  emotional  and  defensive 
reactions  have  often  colored  what  were  to  be  otherwise  scientific 
conclusions. 

In  summary,  the  theoretical  potential  of  wideband  radar  systems  for 
improved  performance  over  conventional  narrowband  radar  systems  has  been 
demonstrated  repeatedly  in  the  literature  (see  references).  This 
potential  is  tempered  by  both  fundamental  and  practical  limitations  in 
implementing  such  systems.  Two  major  and  distinct  areas  of  research  in 
wideband  radar  were  reviewed  in  Sections  4.1  and  4.2.  The  advantages 
and  disadvantages  of  each  were  discussed.  It  was  concluded  that  in 
general,  the  system  concepts  based  on  the  SEM  were  more  robust  to 
unknown  perturbations  in  the  radiated  waveform,  target  shape  and  aspect, 
and  resolution,  than  the  techniques  using  baseband  video  pulses  or 
"nonsinusoidal "  waves  suggested  by  Harmuth.  On  the  other  hand,  the 
former  have  often  not  addressed  the  limitations  imposed  by  present  day 
technology;  In  this  regard,  the  latter  has  introduced  novel  concepts  for 
overcoming  these  limitations.  Future  research  should  aim  at  integrating 
both  of  these  areas  via  an  end-to-end  system  formulation  and  analysis  to 
identify  the  tradeoffs  in  payoff  and  risk  to  better  focus  development  of 
the  critical  technology  areas.  An  example  of  such  an  analysis  was 
presented  in  Section  3  for  the  problem  of  target  detection  in  clutter 
and  noise.  The  tradeoff  between  a  system  constraint  (total  radiated 
energy)  and  the  optimum  signal  spectrum  was  clearly  illustrated. 


5.0 

CONCLUSIONS  AND  RECOMMENDATIONS 


A  preliminary  Investigation  Into  the  potential  advantages  and 
perceived  limitations  of  wideband  radar  (WBR)  system  concepts  has  been 
performed.  A  wideband  radar  Is  defined  qualitatively  as  utilizing 
signals  whose  relative  bandwidth  Is  significantly  larger  than  those  used 
In  the  corresponding  conventional  narrowband  radar  system.  The 
capabilities  of  wideband  radar  systems  for  providing  Improved 
performance  over  current  narrowband  radars  was  examined.  At  the  same 
time,  the  limitations  of  each  approach,  both  fundamental  and  practical, 
were  defined.  The  most  promising  concepts  for  wideband  radar 
applications  and  hardware  components  were  Identified.  The  results  of 
the  research  effort  are  summarized  below.  Recommendations  for  further 
research  conclude  the  section. 

5.1  SUMMARY 

A  generic  set  of  radar  applications  were  Identified  In  Section  2. 
These  Include  search,  target  detection,  tracking  (particularly  at  low 
angles),  high  resolution  imaging  (SAR),  and  target 
classification/identification.  The  capabilities  and  limitations  of 
current  narrowband  radar  systems  for  several  of  these  applications  were 
discussed  in  detail.  In  general,  Improvements  In  performance  could  be 
realized  In  all  areas,  as  might  be  expected,  through  advancement  in  one 
or  more  of  the  following  categories: 

1.  Increased  system  slgnal-to-nolse  ratio  (SNR) 

2.  Finer  range  and  angular  resolution  at  all  ranges 

3.  Improved  clutter  suppression 

4.  Robust,  aspect  and  illumination  independent  target  classifi¬ 
cation 
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While  this  conclusion  may  appear  obvious.  It  Is  Important  In  that  It 
provides  a  quantitative  list  of  performance  parameters  by  which  a 
candidate  WBR  system  can  be  judged. 

Throughout  the  years,  the  utility  of  wideband  or  nonsi nusoldal 
signals  In  radar  applications  has  been  the  subject  of  (often  Intense) 
debate,  particularly  as  regards  the  work  of  H.  Harmuth  and  his 
colleagues  [32-52].  The  research  program  described  herein  was  in  part 
motivated  by  a  need  to  objectively  resolve  the  issues  behind  this 
controversy.  It  became  apparent  rather  early  In  the  program  that  the 
longevity  of  the  controversy  was  promoted  primarily  for  two  reasons:  (1) 
a  tendency  of  some  advocates  to  expound  the  virtue  of  the  proposed  WBR 
concepts  beyond  what  Is  supported  by  their  (often  oversimplified) 
analyses,  and  (2)  the  emphasis  by  both  proponents  and  opponents  on  a 
single  element  of  a  particular  WBR  system  or  application  (e.g.  antenna, 
atmospherics,  etc.)  In  order  to  justify  either  acceptance  or  rejection 
of  the  entire  concept.  It  was  therefore  concluded  that  In  order  to 
provide  an  objective  evaluation  of  a  specific  nonsi nusol dal  radar 
concept,  an  end-to-end  formulation  and  analysis  of  the  system 
performance  is  required. 

An  example  of  such  an  analysis  was  presented  In  Section  3  for  one 
of  the  above-mentioned  generic  radar  applications,  namely  target 
detection  In  clutter.  The  analysis  was  supported  by  a  formulation  of 
the  transfer  function  for  a  general  radar  system  utilizing  signals  of 
arbitrary  spectral  content.  The  details  of  the  formulation  are 
contained  In  Appendix  C.  It  was  shown  that  for  optimum  detection,  the 
shape  of  the  power  spectrum  of  the  transmitted  signal  was  highly 
dependent  upon  the  spectral  transfer  function  of  the  transmitting  and 
receiving  antennas,  target,  clutter,  and  noise.  Furthermore,  even  when 
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these  are  fixed,  the  shape  of  the  signal  spectrum  was  highly  dependent 
upon  the  total  energy  of  the  pulse,  which  was  Included  In  the  analysis 
as  a  system  constraint.  These  results  demonstrate  the  Importance  of  an 
end-to-end  system  analysis,  wherein  the  limitations  of  system  components 
(e.g.  the  antenna)  can  be  traded  off  against  the  advantages  of  enhanced 
target  scattered  and/or  reduced  atmospheric  attenuation.  It  should  also 
be  pointed  out  that  the  results  of  Section  3  show  that  only  the  power 
spectrum  of  the  transmitted  signal  enters  Into  the  definition  of  the 
optimum  receiver;  the  form  of  the  temporal  signal  which  conveys  the 
prescribed  spectrum  is  arbitrary,  thus  avoiding  the  need  to  necessarily 
radiate  ultra-short  video  pulses  to  achieve  large  relative  signal 
bandwldths  (at  least  for  the  problem  of  target  detection).  This 
requirement  has  been  one  of  the  major  focal  points  of  the  nonslnusoldal 
radar  controversy  [53,  54].  While  the  results  of  Section  3  are  not 
completely  general,*  they  do  set  the  stage  for  future  analyses  of  other 
potential  WBR  applications. 

With  this  analytical  groundwork  in  place,  the  research  effort  turned 
to  a  review  of  the  recent  work  In  unconventional,  wideband  radar 
concepts.  The  review,  presented  in  Section  4,  focused  on  research  which 
met  one  or  more  of  the  following  criteria: 

1.  A  novel  approach  which  cannot  be  realized  using  conventional, 
narrowband  systems. 

2.  A  controversial  method  for  which  an  objective,  end-to-end 
analysis  Is  required  for  evaluation. 


*The  target  transfer  function  was  assumed  fixed  and  known,  implying 
a  known  target  and  aspect  direction.  Suggestions  on  how  to  mitigate 
these  requirements  are  discussed  In  Section  3. 


3.  A  technique  which  was  specifically  Identified  for  consideration 
as  part  of  the  research  program  reported  herein. 

In  order  to  limit  the  review  to  the  scope  of  the  program,  a  pair  of 
research  areas  were  selected  which  meet  the  above  requirements: 
techniques  based  on  the  singularity  expansion  method  (SEM),  and  the 
so-called  "nonslnusoldal  radar''  studies  based  on  the  used  of  extremely 
narrow  baseband,  or  video,  pulse  waveforms.  The  former  has  been  used  by 
a  large  number  of  researchers  In  a  variety  of  applications;  those 
specifically  dealing  in  radar  include  the  work  at  General  Research 
Corporation  (GRC)  [55-57,  65]  In  "resonance-region  radar"  (R  ),  studies 
at  Michigan  State  University  [58-60,  66,  67]  on  "radar  waveform 
synthesis,"  and  the  investigations  at  Ohio  State  University  [61-64]  Into 
ramp  response  waveforms  and  the  "kill-pulse"  concept.  The  latter  has 
been  carried  out  almost  exclusively  at  Catholic  University  [32-52]  by  H. 
Harmuth  and  his  colleagues. 

The  research  activities  in  question  could  be  divided  into  studies 
of  wideband  radar  applications  or  wideband  radar  components.  Virtually 
all  the  SEM  work  concentrates  on  the  former,  and  depends  in  some  manner 
upon  the  differences  between  the  characteristic  natural  resonances  of 
dissimilar  targets  to  achieve  target  classification.  Of  significant 
importance  Is  the  fact  that  the  natural  resonances  are  Intrinsic  to  each 
target  and  independent  of  both  the  aspect  direction  and  spatial 
variations  of  the  Illuminated  wave.  This  suggests  that  SEM-based 
classification  schemes  may  provide  robust  performance  In  a  variety  of 
geometries  using  a  small  parameter  set  for  target  characterization 
(l.e.,  a  few  SEM  poles).  While  all  three  of  the  above-mentioned  groups 
have  employed  sound  electromagnetic  theoretical  formulations  to  describe 
their  concepts,  each  employs  a  somewhat  different  approach  to  achieving 
target  discrimination.  Specifically,  each  has  different  criteria  for 


136 


- \ 

;erim 

Lmmd - 


determining  the  required  transmitted  waveform  and  each  uses  different 
algorithms  for  processing  the  receiver  signal  scattered  by  the  target. 
As  Is  to  be  expected  In  preliminary  analyses  of  this  type,  simplifying 
approximations  have  been  made  to  varying  degrees  regarding  the  effects 
of  the  antenna  and  Intervening  medium  on  the  waveforms.  Nonetheless, 
each  demonstrates  the  potential  fcr  significant  Improvements  in  the 
classification  of  unresolved  targets.  Assessment  of  the  relative 
performance  (correct  versus  Incorrect  classification)  of  each  technique, 
however,  requires  an  appropriate  end-to-end  systems  analysis  similar  to 
that  presented  In  Section  3;  such  an  analysis  has  yet  to  be  performed 
for  any  of  the  candidate  methods.  In  fact,  GRC  has  also  Identified 
systems  analyses  as  the  next  step  in  their  study  [57]. 

Harmuth's  work  In  nonsl nusoidal  radar  has  addressed  both  radar- 
applications  [33-46]  and  radar  system  components  [32,  47-51],  As 
discussed  In  Section  4.2.1,  his  approach  to  the  former  Is  often 
oversimplified  and  at  times  misdirected.  Specifically,  there  is  a 
tendency  to  overemphasize  a  single,  sometimes  secondary,  aspect  of  the 
problem,  thus  leaving  the  major  Issues  unresolved.  The  analysis  Is 
almost  always  driven  to  consideration  of  a  specific  radar  waveform,  a 
baseband  video  pulse,  so  that  the  results  are  highly  dependent  upon 
maintaining  the  specific  form  of  the  waveform  after  propagation  through 
the  medium  and  scattering  by  the  target.  This  assumption  leads  to 
detection  and/or  processing  schemes  which  are  not  sufficiently  robust 
for  practical  operating  conditions.  This  is  unfortunate,  since 
Harmuth's  conclusions  that  wideband  radar  offers  significant  potential 
for  Improvement  In  the  performance  of  conventional  narrowband  systems  is 
essentially  correct.  The  deficiencies  In  his  analysis,  however,  have 
led  to  substantial  criticism  of  the  WBR  concept. 
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In  contrast  to  his  treatment  of  radar  applications,  Harmuth's  work 
In  nonslnusoldal  radar  components,  particularly  antennas  [47-49]  and 
waveguide  devices  [50,  51]  Is  thorough  and  rigorous.  Although  more 
research  must  be  done  to  assess  their  full  utility  in  WBR  applications, 
Harmuth's  novel  concept  of  the  large-current  radiation  (LCR)  as  an 
element  for  nonslnusoldal  radar  has  already  been  used  successfully  for 
Into-the-ground  probing  [52],  and  his  demonstration  that  nonslnusoldal 
waves  can  propagate  distortion-free  in  waveguides  and  resonant  cavities 
represent  important,  new  developments. 

A  summary  assessment  of  the  SEM-based  and  nonslnusoldal  radar 
research  reviewed  herein  is  contained  in  Section  4.3.  The  reader  Is 
referred  to  that  section  for  further  details. 

5.2  RECOMMENDATIONS 

The  long-term  potential  of  using  large  relative  bandwidth  signals 
and  wideband  radar  systems  in  military  applications  cannot  be  Ignored. 
Significant  improvement  could  be  realized  in  detection  of  low  observable 
targets,  low-loss  atmospheric  propagation,  robust  target  classification, 
and  high  resolution  over-the-horizon  applications,  among  others.  It  is 
also  true  that  significant  research  is  necessary  to  bring  this  potential 
to  fruition,  particularly  in  two  areas: 

1.  More  detailed  systems  and  performance  analyses  of  theoretical 
WBR  concepts  (e.g.,  resonance-region  radar  [57],  radar  waveform 
synthesis  [58],  kill-pulse  and  ramp  waveforms  [63,  64],  etc.)  to 
determine  the  optimum  implementation  of  an  SEM-based  radar  system  for 
each  candidate  application,  particularly  regarding  the  waveform  design 
and  signal  processing  architecture. 

2.  Further  development,  of  devices  and  components  for  the 
above-mentioned  wideband  radar  designs,  with  specific  emphasis  on 
transmitting  and  receiving  antennas  and  antenna  arrays,  low  nclse 
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wideband  receivers,  and  optimum  signal  processors. 

The  research  efforts  in  these  categories  must  obviously  be 
coordinated,  so  that  critical  technology  requirements  identified  in  (1) 
can  be  properly  addressed  In  (2).  Considering  the  extent  of  the 
research  performed  to  date  (see  the  References)  and  the  fact  that  the 
Department  of  Defense  has  funded  longer-term,  higher-risk  sensor 
concepts  such  as  those  envisioned  for  the  Strategic  Defense  Initiative, 
It  Is  recommended  that  continued  support  be  provided  to  the  most 
promising  WBR  concepts,  such  as  those  based  on  the  SEM,  and  components, 
such  as  the  LCR  and  the  time  domain  sensors  being  developed  at  the 
National  Bureau  of  Standards  [81],  with  the  goal  of  an  eventual 
demonstration  cf  a  prototype  WBR  system  In  a  realistic  radar  application 
scenario.  In  order  to  achieve  this  goal,  It  is  recommended  first  and 
foremost  that  the  detailed  performance  analyses  In  (1)  above  be 
conducted  on  a  comparative  basis  to  define  the  approaches  most  worthy  of 
consideration. 
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APPENDIX  A 

REMARKS  ON  LINEAR  MAPPINGS  AND  FREQUENCY  DOMAIN  ANALYSIS 

Part  of  the  controversy  over  Nonsinusoidal  radar  (NSR)  centers  on 
the  engineering  community's  use  of  sinusoids  and  Fourier  analysis  to 
describe  linear  systems.  The  purpose  of  this  Appendix  is  to  review  the 
properties  of  such  analysis  techniques  and  thus  to  demonstrate  that  the 
techniques  are  not  limiting  the  class  of  radar  systems  that  can  be 
considered.  In  addition,  it  is  shown  that  the  "complex  envelope" 
notation,  commonly  used  in  radar  systems  analysis,  is  not  limited  to 
describing  waveforms  with  narrow  relative  bandwidth. 

A. 1.  LINEAR,  TIME-INVARIANT  MAPPINGS 

Let  L  be  a  linear  mapping  of  a  space  of  functions  into  itself  or 
another  space  of  functions.  (We  proceed  informally.) 

DEFINITION:  A  linear  mapping  L  is  said  to  be  "time-invariant"  if, 

given  that 


y(t)  =  L{x(t)} 


then 


L{x(t  -  t)}  *  y(t  -  t). 
We  introduce  the  notation: 


e(t;  co  )  =  L(e1a,t}. 


Preceding  Page  Blank 
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THEOREM:  If  L  Is  a  linear,  time-invariant  (LTI)  map,  then 


e ( t ;  to)  =  e(o;  ^){e1ajt} 


Proof  -  As  L  is  LTI,  the  following  string  of  equalities  holds: 


S  iwt,  iOJT  .  r  Wt-x)-,  iut  ..  v 

e(t;  to)  =  Lie  }  =  e  L{e  '  '}  =  e  e\t  -  x;  to). 


Setting  t  =  0  and  then  making  the  change  of  variable  (-x)+t. 


e(t;  to)  =  e(0;  to)  e1a)t.  Q.E.D. 


That  is,  the  complex  exponentials  are  the  natural  invariants  of  LTI 
maps,  or  LTI  systems. 


DEFINITION:  £(00)  =  e(0;  to)  is  the  "transfer  function"  associated  with 

the  LTI  map  L. 


Since,  speaking  nonrigorously,  any  function  x(t)  has  a  Fourier 
integral  representation, 

^4-co 

x ( t )  =  ~  I  e,UJ^(o))  d03 


1  f  icot~, 
2tTJ  e  xl 


where  x(oj)  is  the  spectrum  (Fourier  transform)  of  x(t),  one  has  a 
powerful  "operational  calculus"  for  LTI  maps.  For 


y(t)  =  L { x ( t ) }  * 


if 


+00 


L{elajt}x(o>)  doj 


H 

/ 


+DO 


2^  I  elwt  l(aj)x(a>)  dw 
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(The  assertion  or  the  Theorem,  recalling  the  Introduced  notation 
and  deflntlon.  Is 


i(o>)  e(t;  co)  =  L{e(t;  «)} 

that  Is,  e(t;  oj)  =  expOwt)  is  a  "generalized  eigenfunction"  of  L  and 
^ (oj)  Is  the  corresponding  "generalized  eigenvalue. 11  The  w-support  of 
the  spectrum  of  L,  Is  generally  the  continuum. 

A  "linear  antenna"  is  a  more  general  Idea.  A  receiving  antenna  Is 
modeled  by  a  mapping  from  a  function  of  four  variables,  space  and  time, 
to  a  function  of  one  variation,  time.  The  converse  is  true  for  a 
transmitting  antenna.  In  the  latter  case,  we  may  Imagine  a  function  of 
time  v(t)  impressed  on  a  pair  of  terminals,  the  result  being,  say  a 
scalar,  electromagnetic  field  appearing  at  some  suitable  arbitrarily 
chosen  reference  plane  PP',  as  illustrated  In  Figure  A-l. 

DEPINTION:  (1)  A  linear  antenna  structure  has  a  linear  mapping  L  , 

that  is,  for  scalars  a ,  3 , 

L{aVj(t)  +  Sv2(t)}  =  aL  (Vj(t)}  +  BL{  v^C  t )} 

(2)  A  linear  antenna  structure  Is  time-invariant  if  its 
associated  mapping  L  is  such  that,  if 

Lfv( t) )  =  E^r(F,  t). 


then 


L  * v ( t  -  x ) }  =  Etr(r,  t  -  :). 
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Figure  A-1.  Defining  Structure  for  Antenna  (Radiating)  Element 
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We  now  introduce  the  notation 

e(r,  t;  w)  *  L  {eiaJt>. 

We  have 

e(r,  t;  «)  =  L{e1ajt}  =  e1WT  Lfe1^1"1^}  =  e(r,  t-x;  to) 

so  that.  In  the  same  way  as  above, 

e(r,  t;  oj)  =  e(r,  o;  oj  )  e1ajt. 

DEFINITION:  £tr(r,  “)  =  e(r»  o ;  co )  is  the  “transfer  function"  of  the 
associated  LTI  transmitting  antenna.  (Notice  that,  in  this  case,  the 
transfer  function  depends  upon  the  selected  spatial  location  r. ) 

In  a  similar  manner,  a  transfer  function  may  be  defined  for  a 
receiving  antenna  structure. 

A. 2.  LINEAR,  TIME-VARIANT  MAPS 

More  generally,  linear  mappings,  of  both  bounded  and  unbounded 
nature,  can  have  "generalized  eigenvectors  and  eigenvalues",  that  Is, 
spectral  representations.  Generally  speaking,  each  time-variant,  linear 
mapping  has  its  own  natural  spectral  representation.  While  this  can  be 
a  complicated  subject  in  functional  analysis,  even  in  the  classical 
discussion  of  Sturm-Liouville  differential  equations,  this  is  a  well 
known  phenomenon.  For  example,  the  Bessel  functions,  the  Walsh 
functions,  etc.,  are  the  natural  orthonormal  sets  associated  with 
certain  differential  equations  with  time-varying  coefficients. 
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Given  a  spectral  representation,  generally  speaking,  there  will 
exist  an  “operational  calculus"  that  provides  a  spectral  representation 
of  a  polynomial  In  the  original  mapping.  However,  unlike  the  case  of 
linear,  time-invariant  mappings,  polynomial  combinations  of  a  linear, 
time-variant  maps  are  generally  not  of  practical  interest,  and  certainly 
not  as  widely  applicable. 

A. 3.  COMPLEX  ENVELOPE  NOTATION 

As  the  preceding  discussion  has  made  clear,  when  an  overall  system 
is  linear  and  time-invariant,  the  Fourier  integral  representation  of  an 
arbitrary  waveform  as  a  weighted  sum  of  complex  exponentials  or 
sinusoids,  is  natural  and  greatly  facilitating. 

A  separate  matter  that  commonly  arises  In  conjunction  with  this 
representation  is  variously  labeled  "the  complex  envelope 
representation"  or  "In-phase  and  quadrature  representation",  or 
"envelope  and  phase  representation",  or  such.  These  are  very  convenient 
for  both  analysis  and  intuition  and,  are  typically  used  for  waveforms 
whose  spectral  support  is  small  relative  to  their  mean  frequency,  so 
called  "narrowband"  waveforms.  However,  more  generally,  the  following 
may  be  established. 

The  complex  representation  from  which  the  other  representations  may 
be  derived,  is  a  mapping  or  relation  between  a  class  of  real  bandpass 
functions  and  a  class  of  complex  functions.  In  order  to  be  useful,  the 
mapping  must  be  one-to-one.  The  necessary  condition  for  this  is  easy  to 
state  after  a  few  definitions.  First,  the  class  of  real  "bandpass" 
functions  have  their  spectral  support  in  the  frequency  set 

n(<0  a  {co:  -  co  -  <  co  <  -cd  +  cu  -  n  <  «  <  w  +  ft) 

u  u  0  0  o 
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the  collection  of  such  function  Is  denoted  R (oj _ ,  S7 ) .  Second,  the  class 

0 

of  complex  "low  pass"  functions  have  their  spectral  support  In 


=  {so;  | co |  <  ft} 


the  collection  of  such  functions  Is  denoted  C(ft). 
to  R{ojq ,  ft)  is  a  mapping  K  such  that  V  F  e  C(ft), 


KF  =  Re 


F  e 


Ioj  tl 

0  ( 


The  mapping  from  C(n) 


It  is  easy  to  show  that  the  range  of  K  Is  R(£jJ0»  and  that: 


THEOREM:  K  Is  one-to-one  onto  R(fo0,  if  and  only  If  <oc  >  Q,  . 


That  Is,  as  long  as  a  "center  frequency"  (oQ  is  chosen  so  that 
the  complex  representation  Is  unqiue.  Otherwise,  is 

arbitrary.  And,  In  particular,  note  that  the  ratio  oj q/?2  >  l  may  be 
arbitrarily  close  to  1,  that  is,  o)Q  may  be  arbitrarily  close  to  &  . 
Hence,  In  this  sense,  the  real  bandpass  signal  does  not  have  to  be 
"narrowband",  i.e.,  J2/a)Q  <<  1.  Thus  the  often  made  conclusion  that  "if 
complex  notation  is  used,  the  waveforms  must  have  small  relative 
bandwidth"  is  formally  not  correct. 
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APPENDIX  B 

THE  GENERATION  OF  SPECIFIED  SPECTRUM  MODULUS 
WITH  A  SPECIFIED  SIGNAL  TIME  ENVELOPE 

The  solution  to  the  system  design  problem  posed  in  Chapter  3 
resulted  in  the  specification  of  the  optimum  signal  modulation's 
spectral  modulus,  given  an  energy  constraint.  No  other  aspect  of  the 
signal  modulation  was,  at  least  directly,  specified.  As  a  practical 
matter,  the  signal  modulation's  time  envelope  is  usually  a  "rectangular 
pulse",  of  a  duration  T  chosen,  under  a  practical  peak  power  constraint, 
of  duration  sufficient  to  meet  the  energy  constraint. 

One  then  arrives  at  a  version  of  a  classical  problem:  How,  if 
possible  at  all,  to  specify  the  signal  modulation's  phase  modulation  to 

*v 

realize  the  specified  pair  (time  envelope,  spectral  modulus).  One  such 
procedure,  which  applies  where  the  "time-bandwidth  product,  TBP",  that 
is,  for  example,  the  product  of  T  times  the  extent  of  the  support  of  the 
spectral  modulus,  is  large  [19,20]  will  be  described  herein.  Depending 
on  the  energy  restraint  and  the  peak  power  limit,  this  can  be  true  here 
and,  if  so,  presents  a  design  problem. 

The  intuitive  idea  is  this.  Suppose  the  specified  time  envelope  is 
rectangular,  of  duration  T.  Suppose  also  that  the  signal's  spectrum 
modulus  is,  as  sketched  in  Figure  B-la,  the  solution  of  the 
above-discussed  optimization  problem.  When  the  TBP  is  large,  it  is 
possible  to  relate  the  spectral  modulus  to  the  signal's  instantaneous 
frequency  modulation,  X(t),  (the  time-derivative  of  the  signal's  phase 
modulation).  One  imagines  the  "spectral  modulus"  as  a  "density" 
distributed  in  porportion  to  the  relative  time  the  instantaneous 
frequency  spends  near  a  given  frequency.  With  this  idea,  one  may 
construct  the  X(t)  that  will  generate  the  above  spectral  modulus,  to  a 
good  approximation  when  Tn»  1.  (See  the  sketch  in  Figure  B-lb). 
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To  establish  this  approximation  scheme,  express  the  signal 
modulation  In  Its  Fourier  Integral  (spectral)  representation, 


f(t)  -  |f(t)| 


ei9(<o) 


doj 


and  then  evaluate  the  integral  by  the  method  of  stationary  phase.  In 
the  simplest,  and  usually  appropriate  case,  there  Is  one  stationary 
point  and^  the  evaluation  leads  to  a  differential  equation  Involving 
|f(t)|,  |f[X(t)]|,  and  4>  [X(t)].  Involving  the  stationary  condition  on 
<l>"[X(t)],  the  solution  Is  seen  to  be 


X(t)  =  G“1[2ttG(  t)D , 


where 


G(t) 


■S 


f (x) | 6  dr 


:-i 


and  G  Is  the  inverse,  presumed  to  be  definable,  of 

.x 


G(x 


>■/ 


f(u)r 


tn  the  present  application  the  support  of  the  optimal  signal 
modulation  was  seen,  in  the  specific  examples,  to  have  two  disjunct 
supports  (located  near  the  specifically  chosen  object's  "resonances"), 
for  a  certain  Interval  of  energy  constraints.  For  example,  a  spectrum 
similar  to  that  sketched  was  observed.  In  this  case  G(x)  will  appear  as 
sketched  in  Figure  B-2b.  Clearly  G"*  does  not  strictly  exist.  However, 
the  situation  is  readily  handled  by  dividing  the  range  of  G  into  two 
sets,  Y1  =  ^(y:  0  £  y  £  Vj)  and  Y,,  =  (y:  <  y  £  y^  ;  over  each  of 
these  sets,  G-1  exists. 
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a)  Typical  Signal  Modulus 


* 

G 


b)  Cumulative  Distribution  of  Signal  Modulus 
Figure  B-2.  Typical  Signal  Modulus  and  Cumulative  Distribution 


Thus,  a  sketch  of  ?  likely  sufficiently  general  algorithm  to  find 
such  a  frequency  modulation  and  verify  It  is  as  follows: 

1.  F^nd  |f  (cu)i  with  the  numerical  procedure  described  in 

Section  3.4.  |  fQ(w)r  is  stored  as  an  array. 

2.  (r)  Find  the  (finite  number  of)  compact  sets  that  make  up 
the  support  of  |f^w)|  .  Denote  them  A.,  A„,  ...  A^.  These  may  be 
stored  as  a  2  x  N  arvay 

(b)  Calculate  the  inae.'lnlte  integrals 

|  f(w)l‘-  du),  x  e  An,  n  *  1,  . ...  N, 


for  a  set  of  {x„  1  -  1,  ....  I)  and  store  as  an  nx  I-dimensional 

m,  1 

array. 

(c)  Calculate  the  set  of  N  Images  y 

A 

N,  and  the  set  of  images  {y_  .  =  G(x  ,)}. 

in  y  i  m  f  I 

3.  (a)  For  each  t.  e  T,  the  it,}  a  selected,  suitable  set  of 

J  1 

sampling  points  (stored  as  an  array),  calculate 

yj  =  G(ti)  -  j  (tj  +  1/2). 

(b)  Determine  which  image  Ym  contains  y^ .  Then  evaluate 


=  G(<c  ),  n  =  1 . 

r  |hi 


G(,(x)  - 


lo 


2p 


X(tj>  ‘  Sm  <*j>  +  xm|  ; 

>lo 

an  Interpolation  scheme  will  he  required,  e.g.,  find  the  (y  4}  of  the 

Ml  J  1 

array  2(c)  closest  to  y^. 

4.  (a)  Having  the  U(tj),  t^  e  T»  1  =  1,  ...»  J},  integrate  to 

find  the  phase  modulation  4>(t)  at  the  { t^> . 

(b)  As  a  verification,  form  f?(t)  =  jf(t)|  and 

compute  Its  (discrete)  Fourier  transform,  e.g.,  by  an  FFT,  and  compare 
|f>)!2  to  |f0(«.)[2. 


Appendix  C 

GENERAL  TRANSFER  FUNCTION  AND  IMPULSE  RESPONSE 
OF  A  WIDEBAND  (NONSINUSOIDAL)  MONOSTATIC  RADAR 

An  end-to-end  formulation  for  the  input-output  relationship  of  a 
radar  system  operating  with  arbitrary  wideband  signals  is  developed  in 
this  Appendix.  This  type  of  formulation  is  required  in  order  to  fully 
understand  the  tradeoffs  between  the  various  elements  of  such  a  system. 
These  tradeoffs  will  then  allow  an  objective  assessment  of  the 
advantages  and  limitations  of  using  wideband  signals  in  particular  radar 
applications. 

Section  C.l  presents  an  antenna  current-based  representation  of  the 
total  system  transfer  function  founded  on  rigorous  electromagnetic  field 
theory.  All  potential  frequency  dependent  quantities  are  identified. 
Section  C.2  discusses  how  this  representation  is  related  to  the  antenna 
aperture-based  formulation,  which  was  used  in  Section  3  to  demonstrate 
how  the  end-to-end  system  mode!  can  be  applied  to  a  particular  radar 
scenario,  namely  target  detection.  Specifically,  the  antenna  tranfer 
functions  on  transmit  and  receive  are  identified  using  reciprocity,  and 
the  equivalence  of  the  reflectivity  and  scattering  matrix  descriptions 
of  a  target  is  shown. 

C.l.  SYSTEM  FORMULATION 

A  block  diagram  of  a  general  (monostatic)  radar  system  is  shown  in 
Figure  C-l.  The  pertinent  inputs  and  outputs  of  each  major  component 
are  identified  symbolically  and  will  be  defined  in  detail  below.  It  is 
desired  for  analytic  tractabi 1 ity  that  the  overall  system  be  linear;  for 
this  reason  the  quantities  of  interest  will  be  voltages,  currents,  and 
fields,  as  opposed  to  powers,  cross-sections,  etc.,  which  are 
proportional  to  the  square  of  the  former.  Without  a  significant  loss  of 
generality,  the  system  input  is  assumed  to  be  a  voltage  waveform 
produced  by  an  appropriate  signal  generator,  while  the  output  is  also  a 
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Figure  C-l.  Block  Diagram  of  General  Radar  System 
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voltage  corresponding  to  that  which  effectively  appears  at  the  terminals 

★ 

of  the  load  resistance  of  the  receiving  antenna.  Contrary  to  claims 
made  elsewhere  [39  ],  the  analysis  is  most  straightforwardly  conducted 
in  the  frequency  (u>)  domain,  where  to  is  an  angular  frequency  variable 
associated  with  the  temporal  Fourier  transform.  This  transform  is 
defined  as 


F(w)  "  /  f(t)  e,0Jtdt,  (C-l  ) 

-  00 

and  its  inverse 

00 

f(t)  -•  2 7T  J  F(w)  e"iwtdaj.  (C-2) 

•  oo 

Some  comments  are  appropriate.  First,  f(t)  and  F(co)  can  in  general  be 
distributions,  so  that  the  existence  of  Eqs.  (C-l)  and  (C-2)  is  not  an 
issue.  Second,  the  temporal  functions  must,  satisfy  causality,  i.e. 

f(t)  »  0,  t  <  0, 

thus  allowing  the  lower  limit  in  Eq.  (C-l)  to  be  replaced  by  t  =  0. 
Although  a  LaPlace  transform  is  often  used  in  this  situation,  it  is  not 
necessary  to  do  so,  and  Fourier  transform  will  be  used  here  instead. 

The  radar  transmitter,  defined  by  ari  input  voltage  waveform  V.  (to), 
produces  an  output  current  I  (to)  which  is  applied  to  the  terminals  of 

a 

the  transmitting  antenna.  It  is  convenient  to  write 


*Some  amount  of  amplification  of  other  linear  processing  is 
allowable  between  the  load  resistance  and  the  system  output  terminals. 
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(C-3) 


v1n(u)  *  [ztio.)  ♦  za{(0)]iaU), 


where  (co )  «  ( co)  -  1  X^(oj)  =  trans impedance  of  transmitter, 


Za(w)  *  R_ ( uj)  -  I X_ (uj)  *  input  impedance  of  the  antenna. 

o  d  a 


For  convenience,  R  and  X.  are  the  radiation  resistance  and  reactance  of 

a  a 


the  antenna,  respectively.  All  other  contributions  to  the  antenna 
impedance  (in  particular,  internal  losses)  have  been  assumed  lumped  with 


*t* 


By  definition,  then. 


v, „(<■>) 

zt<“)  *  r 


(  C  -  4  ) 


Let  the  current  I  be  distributed  among  N  infinitesimal  radiators 

01 


which  constitute  the  antenna  array.*  This  distribution  is  to  occur  In 
such  a  manner  that  the  frequency  dependence  of  each  current  radiator  is 


the  same  as  I_(u>).  Specifically,  the  vector  current  distribution  of  the 

a 


antenna  array  has  the  form 


N 


V 


Ja(?,w)  =  Ia(co)  ^  ans(f  '  fn)p(rn}t 
n~-1 


(  C  -  5  ) 


where  p  is  a  unit  vector  describing  the  orientation  of  each  radiating 


element.  For  consistency,  the  units  on  -5(r  -  r  )  are  inverse  distance 


cubed  (a  density),  while  those  of  an  are  distance.  In  this  manner, 


1h 


anp(rn)  is  the  dipole  moment  of  the  n  element.  In  general,  the  an  are 


complex,  to  allow  for  a  phase  shift  between  elements.  Furthermore,  this 


*More  general  current  distributions  will  be  considered  later  in 
this  memo. 
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phase  can  be  at  most  linearly  dependent  on  frequency,  thereby  allowing  a 
time  delay  to  be  specified  between  elements  In  the  time  domain. 
According  to  the  discussion  above,  it  Is  required  that 

N 

2>J  '  !•  tC-6) 

n  =  1 

The  distribution  Ja  can  be  made  continuous  by  letting  N  +  “  in  an 

a 

appropriate  limit;  however,  for  NSR  systems,  the  expected  antenna 
designs  employ  discrete  radiating  elements  [47-49]. 

The  far-zone  electric  field  E1  ( r , oi )  generated  by  J  radiating  into 
a  linear,  homogeneous,  isotropic,  dispersive  medium  whose  permittivity, 
permeability,  and  conductivity  are  given  by  e  (u>),y(o>),  and  a(u>)  is 


E^r.w)  = 


k2M  e^Hr 
4ir1w«(ui)  p 


rxrx 


f T.tFiw)  e-,kt“)?-r' 


.3 

d  r 


(  C- 7  ) 


_  e  v  '  f  v  -  -  x  -ik (u)r.r 

-  “?7r - r -  V  )  ^  anrxrxp(rn)e  v  ;rrn, 

n = 1 


where 


Mw) 

eU) 


rr  =  far-zone  position  vector 

w^e{oj)u(u>) 


eo(u) 


+  i 


o  (to) 
0) 


(C-8) 
(  C  -  9  ) 


and  V  is  a  volume  enclosing  the  antenna.  Equation  (C-7)  is  valid  when 
r,  r‘  satisfy  the  far-zone  criterion 


cs  r  - 


(C-10) 
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r  -  r 


Equation  (C-10)  Is  generally  valid  when  the  origin  of  coordinates  is 
taken  to  lie  near  the  center  of  the  antenna  array.  Under  these 
circumstances,  r  can  be  considered  to  specify  an  angular  direction 
relative  to  the  antenna  in  the  far-field,  as  shown  in  Figure  C-2. 


Let  the  field  if*  be  Incident  upon  a  scattering  target  whose  center 
is  given  by  the  position  r.  The  target  is  bounded  by  a  surface  A,  and 
for  the  moment  Is  characterized  by  arbitrary  consti tui ti ve  parameters 
p^,  at  which  may  be  functions  of  both  position  and  frequency.  It  is 
also  assumed  that  the  target  is  sufficiently  small  that  the  antenna 
resides  in  the  far-zone  of  the  target,  thereby  satisfying  a  relationship 
analogous  to  Eq.  (C-10). 


For  all  practical  situations,  the  field  Es  backscattered  by  the 
target  can  be  described  by  a  scattering  matrix  such  that 


P(r\<-)  = 


ik(<*>)r  = 


4?r  r 


S' E1  (r,w)  e 


-ik(w)r«  r 1 


(C-ll  ) 


where  S  is  a  2  x  2  matrix  in  dyadic  form.  Note  that  S  is  defined 
relative  to  the  target  center,  that  is,  with  respect  to  the  position 
vector  r.  In  this  manner,  the  propagation  phase  of  the  scattered  field 
at  the  antenna  phase  center  ( r *  =  0),  from  Eqs.  (C-9)  and  (C-  11),  is 
2k(w)r,  as  expected. 


The  elements  of  S  depend  strongly  on  the  incidence  direction  r,  the 
target  orientation  (aspect),  and  the  frequency  u>.  The  latter  dependence 
is  particularly  significant  here,  and  when  necessary  will  be  shown 
expl icitly. 
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The  voltage  Vj  { oj)  appearing  across  the  terminals  of  the  receiver 
when  connected  to  a  load  Impedance  Zj(ui)  Is  best  explained  with  the  aid 
of  Figure  C-3.  In  Figure  C - 3 ( a ) ,  the  antenna  Is  operated  as  a 


transmitter,  producing  an  Input  antenna  current  Ia ( o>)  In  response  to  an 


Input  voltage  V^n(oj).  Equation  (C-3)  characterizes  this  relationship. 
In  Figure  C-3b,  the  antenna  is  operating  as  a  receiver,  with  the  load 
Impedance  ZL  replacing  the  Input  voltage  generator.  In  addition,  a  new 


voltage  source  V  (w),  generated  by  the  scattered  field  Es,  serves  to 


drive  the  circuit.  From  the  theory  of  receiving  antennas  [79],  the 
voltage  across  the  load  Is 


vL(«)  =  zlMil(u>) 


Z|  («)  vs(«) 


zLfc)  +  zt(«J“+  z^r 


( C  -  1  2  ) 


where 


VC(J  *  h(r,  w)  •  IV  -  0,  u>). 


(  C  -  1  3  ) 


The  parameter  h  is  the  vector  effective  height  of  the  antenna,  and  is 
given  explicitly  by 


h(f'“)  =  rb|'Vr-“)  e'U(“),:,r'<J3r' 


A 

=  D  anp(rj 


-i  k(w)r- r 


(C-i 4) 


N  =  1 


n-  *  rv 


The  units  of  h  are  length,'  whicn  is  why  it  is  referred  to  as  an 
effective  height. 
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The  final  step  in  the  end-to-end  analysis  is  the  inclusion  of  a 
linear  filter  operating  on  the  voltage  VL  to  produce  a  voltage  Vp. 
Denoting  the  transfer  function  of  the  filter  by  H(w),  it  follows  that 


VF(u>)  -  H(o3)VlU 


(  C  -  1  5  ) 


In  this  respect,  ZL  can  be  considered  the  input  impedance  of  the  filter. 
By  combining  Eqs.  (C-7)  and  (C-14),  one  can  write 


Ei(r,  w)  =  mLwleik(a))r  j  x  f  x  h(r,  u).  (C-16) 

4  tt  1  r  3 


The  scattering  matrix  S  is  transverse  to  the  direction  of  propagation  r, 
i  .e. 


r  *  S  =  S  •  r  =  0, 


(  C  -  1  7  ) 


and  since 


A  A 


r  x  r  x  h  »  r(r  *  h)  -  h. 


(  C  -  1  8  ) 


it  follows  from  Eqs.  (C-16)  through  (C-18)  that 

==  i  iu)u(u))  ik((o)rT  ,  .  -r, a  ) 

S  E  -  e  I, l u))S  n(r,w  ;• 

4  tt  r  a 


(  C  -  1  9  ) 


The  total  input-output  relationship  for  the  system  can  be  found 
from  Eqs.  (C-3),  (C-ll)  through  (C-13),  and  (C-19),  namely 


Z,  («) 


VF(w)  =  H(w)  Z~"(w')  +  zju)  +  Z  («')  (Zt(w)  +  za(w)) 

Led 


-1 


{ C - 20  ) 


IwmM  pi2k(w)r 

(4*r)2 


F(r,u»).  S-K(r,w)Vin(w) 
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The  time-domain  relationship  is  simply 


vr(w)  =  /  VpHe'^dw 


Tin(t,)9(t-t,>dt'» 


(C-21  ) 


where  v^n(t),  g(t)  are  the  Inverse  Fourier  transforms  of  V^n(w)  and 

VF(«) 

sH  5  v~T^y  •  (c-22> 

resoecti vely .  G(w)  and  g(t)  are  the  transfer  function  and  impulse 
response  of  the  radar  system.  Causality  and  the  assumption  that  v^n  ~  0 
for  t<  0  have  been  used  to  determine  the  limits  of  integration  in  Eq. 
(C-21). 

The  model  for  the  system  characterized  by  Eqs.  (C-20)  and  (C-21)  is 
sufficiently  general  to  conform  to  most  noiseless,  backscatter, 
co-polarized  radar  systems  over  their  linear  operating  region.  The 
inclusion  of  a  random,  time-varying  noise  component  can  be  Inserted 
either  before  or  after  the  filter  H(o;)  by  addition  to  or  Vp, 
respectively  (sec  Figure  C-l).  The  generalization  to  a  monostatic 
and/or  bistatic  geometry  with  arbitrarily  polarized  receiver  is 
straightforward  through  identification  of  separate  transmitter  and 
receiver  vector  effective  heights. 
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C.2  TRANSFER  FUNCTION  RELATIONSHIPS  FOR  ANTENNAS  AND  SCATTERERS 

In  this  section,  the  relationships  for  the  transfer  functions  of  a 
transmitting/receiving  antenna  are  provided.  The  role  of  reciprocity  In 
transmission  and  reception  Is  defined.  The  nature  of  these  transfer 
relationships  Is  that  they  avoid  the  need  for  explicit  knowledge  of  the 
current  or  aperture  field  distribution  of  an  antenna*  In  order  to  relate 
the  excitation/received  voltages  to  the  far-zone  radiated/scattered 
fields.  Instead  these  distributions  are  implicitly  contained  in  the 
definition  of  the  "pattern  factor"/vector  effective  height  of  the 
antenna,  as  will  be  shown  below.  These  quantities  are,  in  general,  more 
"measurable"  than  the  above-mentioned  distributions.  Thus,  they  are 
more  useful  in  systems  analysis  for  a  specified  antenna;  they  are  less 
useful,  however,  for  analyzing  or  designing  an  antenna  for  a  particular 
application  or  to  a  given  specification  because  of  their  "implicit" 
nature. 

Attention  is  also  be  given  to  target  scattering  transfer 
relationships,  in  an  attempt  to  cast  the  reflectivity  (Section  3)  and 
scattering  matrix  (Section  C.l)  models  into  a  common  framework.  The 
utility  of  each  model  will  be  discussed. 


C.2.1.  THE  ANTENNA  TRANSFER  FUNCTIONS 

The  far-zone  electric  field  radiated  by  an  antenna  occupying  a 
volume  V  is  [84] 

•  _  i  k  r  _ 

^(r,  oj  )  =  -  ik- -  Zrt(r  xrxN-rxM)  ( C-23 ) 

47Tr  0 

_ F  =  rr 

*These  explicit  distributions  were  used  in  the  formulations 
presented  in  Section  C.l  and  Section  3,  respectively. 
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where  N,  M  are  the  electric  and  magnetic  current  moments  defined  by 


(C-24) 


C  C- 25 ) 


and 


Z 


o 


intrinsic  impedance  of  the  medium. 


The  inclusion  of  the  magnetic  current  density  Jm  is  a  generalization  of 
the  results  of  Section  C.l  to  allow  the  definition  of  equivalent  current 
distributions  for  apertures,  etc.  As  before,  the  Intrinsic  impedance  Z0 
and  the  wave  number  of  the  medium  k  are  assumed  to  be  functions  of  the 
angular  frequency  ou. 


It  is  convenient  to  define  a  pair  of  quantities  F  and  h,  the  field 
pattern  factor  and  the  antenna  vector  effective  height,  such  that 


•  i  kr 

e'(F,  <0)  =-iZ0iaM 


f(r,  w) 


( C-26  a ) 


i  kr 

=  IFF**1"'  “>• 


( C- 2  5  b ) 
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The  pattern  factor  Is  a  dimensionless  quantity  which  describes  the 
angular  variations  of  the  far-zone  field  at  the  frequency  w,  while  the 
vector  effective  height  Is  a  parameter  with  dimensions  of  length  which 
represents  the  effective  receiving  length  of  the  antenna  as  a  function 
of  angle  and  frequency.  Both  quantities  also  contalR  the  full 
polarization  properties  of  the  antenna  and  its  radiated  field.  As 
before,  I  is  the  input  current  at  the  antenna  terminals. 

O 

From  Eq.  (C-26)  it  is  clear  that 

F(r,  oj )  =  kh(r,  co),  (C-27) 

which  is  simply  a  statement  of  the  Rayleigh-Carson  reciprocity  theorem 
for  a  linear  antenna*  in  a  linear,  isotropic  medium. 


If  the  antenna  excitation,  referred  to  the  antenna  terminals,  is 
characterized  by  a  Thevenin  equivalent  circuit  with  an  input  voltage 

Vin(u)  and  a  source  impedance  Z ^ p ( co) , *  it  follows  that 


E1  (r ,  w)  =  Vin(oo) 


7a^“T 


i  kr 


7-W 


4irr 


-F ( r ,  oj), 


(C-28) 


where  Z,  is  the  antenna  input  impedance.  A  transmitting  transfer 

o 

function  is  now  defined  as  [84] 


u  /'l  \  _  r£  (r,  u»)  - 

Hy(r,  w)  v.  (u)  e 

in 


i  kr 


Za  +  Zin, 


F(r,  oj)  , 


(C-29) 


*The  term  linear  antenna  is  used  in  the  sense  of  a  linear  (as  opposed  t 
nonlinear)  system.  It  does  not  mean  the  antenna  is  necessarily  wirelike. 


which  can  be  viewed  as  the  range-normalized  far-zene  field  response  to  a 
unit  input  voltage. 


In  a  similar  manner,  a  receiving  transfer  function  HR  can  also  be 
defined.  Recall  that  the  output  voltage  across  a  load  impedance 
connected  to  the  antenna  terminals  can  be  written  as 


VL(<o) 


( C-30 ) 


where  E s  ( co )  =  E s  ( r  *  0,  oj)  is  the  scattered  field  at  the  antenna 
terminals.  Without  loss  of  generality,  the  antenna  terminals  are 
assumed  to  be  located  at  the  origin  of  coordinates  (F  =  0).  The 
receiving  transfer  function  is  defined  such  that 


VL(w)  *  HR(f,  w)  •  Es(w) 


(C-31  ) 


Thus,  Hr  is  the  response  to  a  unit  amplitude  scattered  field  polarized 
in  the  same  sense  as  the  receiving  antenna.  Explicitly, 


(C-32) 


and  h  is  the  antenna  polarization  unit  vector. 


A  fundamental  result 
is  real,  namely 


can  be  derived  for  the  special  case  when  ZL  = 


id 

4ttZ^ 


to. 


(C-33) 


*Both  Zin  and  Z L, 
transimpedance  Zt  defined 


Eq.  (C-30), 
in  Section  C. 1. 


include 


the 


transmi tter 


2p 


where  u  is  the  permeability  of  the  medium  and  the  superscript  denotes 
the  scalar  component.  Since  in  practice  u  is  nondi spersi ve 

(Independent  of  u>),  Eq.  (C-33)  is  a  statement  of  the  well-known  fact 
that  the  transmitting  impulse  response  of  an  antenna  is  the 
time-derivative  of  the  receiving  impulse  response.  The  fact  is  also 
evident  from  Eq.  (C-27). 

A  final  comment  regarding  the  transfer  function  model  of  Eqs. 

(C-29)  and  ( C - 3 1 ) .  In  the  analysis  of  Section  3,  an  assumption 

regarding  the  spatial  and  frequency  dependence  of  the  equivalent  current 

distributions  J  and  J  of  the  form 
e  m 


W'  “>  '  Je>n,(r)a(“) 


( C-34) 


was  made.  It  was  argued  that  this  assumption  is  valid  for  horn-like  or 
reflector-like  antennas.  While  this  may  be  true  over  narrow  or  even 
moderate  bandwidths,  it  is  in  general  not  applicable  for  wideband 
systems  of  large  relative  bandwidth  such  as  those  being  considered  in 
this  research  program.  It  can  remain  valid  for  array  antennas,  however, 
when  the  individual  elements  are  "point-like"  radiators  and  interactions 
between  elements  are  negligible.  This  was  exactly  the  antenna  type 
considered  in  Section  C.l.  However,  satisfaction  of  Eq.  (C-34)  may  not 
be  necessary  for  the  results  of  Section  3  to  remain  valid. 

C.2.2.  SCATTERING  TRANSFER  RELATIONS 

To  complete  the  model  of  a  wideband  radar  system,  it  remains  to 
define  the  connection  between  the  incident  and  scattered  fields  E1  and 
E"s ,  respectively.  A  rigorous,  explicit  relationship  requires  the 
solution  of  a  boundary  value  problem,  generally  in  the  form  of  an 
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Integral  equation,  over  the  surface  of  the  scatterer.  This  is  a 
difficult  problem  except  in  few  situations,  particularly  if  the  effects 
of  a  clutter  background  are  included.  An  implicit  relationship  that  is 
often  used,  as  it  was  in  Section  C.l,  is  to  define  a  (dyadic)  scattering 
matrix  S,  whereby 

i  kr  _ 

E5^)  *  S[u)-I\r,  o>).  ( C-35 ) 

As  with  the  antenna-based  functions  HT  and  H  ,  S  is  a  transfer  function 

-s  T  r 

which  relates  the  output  E  of  a  linear  system  (the  scatterer)  to  the 
system  input  E1 . 

Now  the  scattered  field  Es  satifies  a  relationship  similar  to  Eq. 
(C-23) ,  i.e. 


i  kr  _  _ 

■ikZ  - - (r  x  r  x  N.  +  r  x  M.  )  , 

w  ,  l*  L 

4  irr 


( C-36 ) 


where  N^,  are  the  electric  and  magnetic  current  moments  induced  on 
the  target  by  the  incident  field.  If  Jtfi  and  Tt  are  the  equivalent 
electric  and  magnetic  currents  defined  on  the  target  surface  S,  then, 
analogous  to  Eqs.  ( C- 24 )  and  (C-25), 


v/ 


Jte(V  aj)  9 


ikr-r 


o  .2 


d  r 


V  '  Z0fJtJr 0 


_  ikr-r 

(r  ,  a,)  e  0  j2 


d  rQ, 


( C-37  ) 


( C-38 ) 
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where  rQ  Is  a  local  position  vector  relative  to  the  target  "center." 
Equations  (C-36)  through  (C-38)  have  been  evaluated  In  the  backscatter 
direction,  which  Is  equivalent  to  replacing  r  by  -r  In  Eqs.  (C-24)  and 
(C-25).  It  Is  evident  from  Eqs.  (C-35)  though  (C-38)  that  the  explicit 
form  of  1>  Is  tied  to  the  relationship  between  the  current  densities 
Jtm  and  the  direction  of  Incidence  r,  and  can  be  determined  only  by  the 
solving  the  above-mentioned  Integral  equations  (which  Is  in  general  not 
possible  analytically). 

Consider  the  possibility  of  reconciling  the  scattering  model  of  Eq. 
(C-35)  with  the  "reflectivity"  model  used  in  Section  3.  First,  it  is 
important  to  note  that  Eq.  (C-35)  assumes  the  incident  field  to  be 
locally  plane  in  the  vicinity  of  the  scatterer  (which  is  consistent  with 
the  far-zone  approximation). 

The  model  used  in  Section  3  for  the  received  signal  from  the  target 
is  of  the  form  (in  the  present  notation) 

vHU)  “  2^  J  Vw)  e"1Wt  dt,  ( C- 39 ) 


where  is  the  spectrum  of  received  voltage  (analogous  to  VL  in  the 
present  model ) .  Explicitly, 


V“> 


a(u))Vin 


i2kr 

e 


w) 


/■ 


W 


i  2kr  •  r 
.  o 

a)  j  c 


(C-40) 
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Here,  a(uj)  »  ay(u)ap(oj)  is  an  antenna-related  spectral  transfer  function 
and  E_  is  a  far-zone,  field-like  distribution.  Of  course,  gt  Is  the 
target  reflectivity  density.  In  order  to  express  these  quantities  In 
terms  of  the  parameters  used  in  this  appendix,  it  is  first  noted  that 
from  the  far-zone  assumption,  and  the  requirement  of  a  localized  target, 
Eq.  (C-22)  implies 

t  (r  +  rQ,  w)  *  F(r,  u>)e  ,  (C-41) 

that  is,  the  Incident  field  is  locally  plane,  as  assumed.  Thus,  from 
Eqs.  ( C- 27 )  and  (C-28) 


i kr  ^  i kr • r 

E1  (r  +  rQ,  ou)  =  -ikVin(ai)  aT(a>)  h(r,  «)  e  0 


(C-42) 


a  i<r-r 
=  EQ(r,  u)h  e 


where  the  unit  vector  h  is  defined  in  Eq.  (C-32). 

Since  the  scattering  process  is  (assumed)  linear,  it  follows  that 
Jtp  and  Jtm  are  proportional  to  the  amplitude  of  the  incident  E1.  By 

i  nt  r.trtiir  i  nn  t  hp  nn  rma  1  i  7o/i  i  nHurori  rn  rron  f  none  itior  i  w  a  a 

j  -te.m 


ikr-r 


te  ,m 


=  EQ( r ,  ^)e 


o 


^te,n/  ro’ 


( C  -  4  3 ) 


and  the  corresponding  co-poiarized  "reflectivity  density," 
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( C-44 ) 


MV"  )  ■  b.r  x  (r  x  Ju  ♦y-Jtm)  • 

0 

Eq.  (C-30)  for  the  output  voltage,  combined  with  Eqs.  (C-36)  through 
(C-3S)  and  (C-41)  through  (C-45),  becomes 


VLU)  *  aT(u)aR(w)Vin(u) 


|H(r,  ) !  2 


/-12kr-r 
9t(r0.  «)  e 


lC-«5) 


where 

ZL  (co ) 

aR(w)  =  ya.)  "+  ZL(U)  • 


(C-46) 


The  similarity  to  tht  result  of  Section  3,  Eq.  (C-40),  is  evident. 
However,  it  must  be  emphasized  that  g.  is  also  a  function  of  the  aspect 

A  * 

direction  r,  and  of  the  form  of  the  incident  field  (here  assumed  to  be 
planar). 


Whether  or  not  it  is  necessary  to  introduce  the  concept  of 
reflectivity  into  the  analysis  of  Section  3  is  not  immediately  evident. 
In  fact,  a  quantity 

~  p  -  f  _  i2kr-r  0 

Eo  (r,  j)  J  9t(r0,  co)  e  0  d2rQ,  ^C-47) 

S 
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Is  defined  therein  such  that 


VL(w)  *  a(id)Vin(oo)e12kr  ut(r,w)  . 


( C- 48 ) 


From  Eqs.  (C-27)  though  (C-30)  and  (C-35),  one  can  easily  deduce 


ut(r,  co) 


(C-49) 


from  which  a  relationship  between  the  reflectivity  and  scattering  matrix 
can  be  obtained,  viz. 


i  s'' 

^•§-h 


A 

12kr* r 
co)  e  0 


(C-50) 


It  appears  that  u^  is  the  only  important  parameter  for  the  approach 
used  in  Section  3,  it  is  therefore  possible  that  Eq.  (C-49)  can  be  used 
in  place  of  Eq.  (C-47). 
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